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A  comparative  study  is  made  of  various  methods  for  computing  the  free 
vibration  modes  and  natural  frequencies  of  thin  plates  with  clamped  and  rota¬ 
tional  supports  and  cylindrical  curvature.  The  methods  include  closed  form 
analytical,  digital  computer,  nomographic,  and  graphical  computations.  Based 
on  the  results,  preferred  methods  of  computation  are  recommended.  These 
methods— Option  2— are  of  particular  value  in  extending  previously  formulated 
digital  computer  programs  for  obtaining  the  vibroacoustic  response  to  turbu¬ 
lence  excitation  of  a  plate.  Computer  results  for  a  particular  case  provide  a 
comparison  of  the  effect  of  clamped-clamped  and  simply  supported  boundaries 
on  the  vibratory  response  of  a  plate  subject  10  turbulence  excitation. 

ADMINISTRATIVE  INFORMATION 

This  study  was  conducted  at  the  Naval  Ship  Research  and  Development  Center  (NSRDC) 
and  supported  by  the  Naval  Ships  Systems  Command  (NAVSHIPS)  Code  0311.  Funding  was 
provided  by  NAVSHIPS  0311  under  Subprojects  S-F1453  21  06  and  R  003 03,  Task  15326. 

INTRODUCTION 

Reference  1*  documents  four  available  computer  programs  for  determining  the  vibratory 
response  and  associated  acoustic  radiation  of  a  finite  rectangular  plate  to  fully  developed 
turbulence  excitation.  Reference  2  treats  a  modification  of  these  computations  to  include  the 
effects  of  pressure  pickup  dimensions  and  boundary  layer  thickness  (Option  1).  These  pro¬ 
grams  include  the  response  of  simple  and  clamped  plates  in  air  and  in  water.  Several  compu¬ 
tational  frameworks  are  provided  which  can  be  modified  and  extended  through  additional  re¬ 
search  to  furnish  more  accurate  programs  capable  of  meeting  naval  needs  in  an  increasingly 
realistic  manner.  The  chief  objective  of  the  original  study  was  to  furnish  a  base  for  future 
development. 

Reference  1  contains  vibroacoustic  solutions  for  all  programs  using  simply  supported 
plate  boundaries  and  for  the  following  programs  using  clamped  plate  boundaries: 

1.  Boeing  Program  I  (Maestrello) 

2.  Boeing  Program  II  —  Finite  Element  (Jacobs  and  Lagerquist) 

3.  Electric  Boat  Program  (Izzo  et  al.) 

Boeing  Program  I  uses  the  Warburton  method  for  computing  the  modes  and  natural  fre¬ 
quencies;  it  may  not  be  adequately  accurate  for  square  plates  or  preferable  with  respect  to 
accuracy,  computer  running  time,  computer  cost,  and  ease  of  computation  etc.  compared  to 


♦References  are  listed  on  page  149. 


other  methods  of  computation.  The  finite  element  method  of  Boeing  Program  II  yields  results 
whose  accuracy  decreases  with  mode  number.  Finally,  the  particular  aspect  of  the  Electric 
Boat  Program  which  deals  with  the  normal  modes  and  frequencies  of  clamped  plates  is  con¬ 
sidered  proprietary  by  General  Dynamics  Corporation;  hence  although  their  numerical  results 
for  a  particular  clamped  plate  computation  are  accessible,  the  associated  program  is  not 
available  to  NSRDC.  Nor  are  other  programs  for  obtaining  the  response  of  clamped-clamped 
plates  presently  available  at  NSRDC.  Thus,  there  is  a  need  for  evaluating  methods  for  ob¬ 
taining  the  normal  modes  and  natural  frequencies  of  clamped  plates  in  order  (1)  to  select  a 
method  or  methods  which  are  relatively  accurate,  simple  to  apply,  and  inexpensive  to  run  on 
a  computer  (if  necessary)  and  (2)  to  extend  the  applicability  of  those  programs  in  Reference  1 
which  are  presently  limited  to  the  case  of  simply  supported  boundaries. 

Accordingly,  the  present  report  presents  a  modification  (Option  2)  of  any  of  the  pro¬ 
grams  of  Reference  1  for  continuous  thin  plates.  The  modification  is  an  attempt  to  incorpo¬ 
rate  into  the  programs  accurate  methods  for  computing  the  normal  modes  and  natural  frequen¬ 
cies  of  plates  with  clamped  and  rotational  supports.  A  method  is  also  presented  for  including 
the  effects  of  clamped  thin  plates  with  cylindrical  curvature  in  the  modified  programs.  The 
selected  methods  for  the  clamped-clamped  finite  rectangular  plate  are  based  on  a  comparison 
of  experimental  results  to  results  of  closed  form  analytical,  digital  computer,  nomographic, 
and  graphical  computations. 

The  following  titles  identify  the  methods  treated  in  the  comparative  study  and  their 
location  in  the  report;  notations  relevant  to  each  method  are  also  included  in  the  Appendixes. 

Appendix  A  -  Warburton  Method 

Appendix  B  —  Young  Method 

Appendix  C  -  Ballentine-Flumblee  Method 

Appendix  D  —  Greenspon  Method 

Appendix  E  -  White  Method 

Appendix  F  —  Crocker  Method 

Appendix  G  —  Sun  Method 

Appendix  H  —  Claassen-Thorne  Method 

The  corresponding  computer  programs  and  flow  charts  are  given  in  Appendix  I. 

Fcr  the  convenience  of  the  reader,  the  Appendixes  include  an  adequate  amount  of 
mathematical  development  underlying  these  methods.  An  understanding  of  the  development 
will  assist  the  reader  to  appreciate  the  merits  and  shortcomings  of  a  particular  method  and  to 
compare  and  apply  the  various  methods.  Relevant  figures  and  tables  are  adapted  from  the 
basic  references. 

In  addition  to  the  references,  a  bibliography  of  other  pertinent  published  papers  is 
given  for  background  information. 
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DISCUSSION 


All  of  the  computer  programs  in  Reference  1  include  a  treatment  for  determining  the 
vibroacoustic  response  for  simply  supported  plates  subject  to  turbulence  excitation.  How¬ 
ever,  both  theory  and  experiment  suggest  that  when  properly  interpreted,  these  programs  can 
also  be  used  directly  to  obtain  the  response  for  clamped  plates.  The  interpretation  is  based 
on  the  following  considerations. 

As  discussed  in  Appendix  C  of  Reference  1,  Izzo  compared  the  computed  sound  pres¬ 
sure  level  for  a  clamped-clamped  plate  with  that  of  a  simply  supported  plate.  The  compari¬ 
son  suggests  that  a  simplified  and  realistic  approach  to  the  investigation  of  plates  with 
nonsimple  supports  would  be  to  calculate  the  modal  frequencies  considering  the  true 
(clamped-clamped)  end  conditions  but  to  use  the  mode  shapes  considering  the  end  conditions 
to  be  simple  supports.  This  approach  requires  much  less  computation  and  its  results  are  in 
very  good  agreement  with  those  of  the  exact  approach  (clamped-clamped  frequencies  and  mode 
shapes). 

Snowdon3  lends  further  theoretical  confirmation  to  these  findings.  He  discusses  the 
first  few  modes  of  a  clamped-clamped  beam*  harmonically  driven  at  its  miapoint.  When  this 
beam  vibrates  in  its  first  four  resonant  and  first  four  antiresonant  modes,  its  displacement 
curves  are  closely  similar  in  appearance  to  those  of  a  simply  supported  beam.  At  the  ends 
of  the  clamped-clamped  beam,  however,  the  slope  as  well  as  the  displacement  of  the  beam  is 
constrained  to  zero.  The  results  for  the  simply  supported  and  clamped-clamped  beams  differ 
principally  in  the  frequencies  at  which  the  resonant  and  antiresonant  modes  of  beam  vibration 
occur. 

Other  investigators  have  found  that  nodal  lines  on  plates  may  be  equivalent  to  simple 
supports,  i.e.,  a  plate  with  any  boundary  conditions  oscillating  in  one  of  its  higher  modes 
thus  behaves  virtually  like  a  slightly  smaller  plate  on  simple  supports.  Moreover,  the  effect 
of  boundary  conditions  on  the  natural  frequencies  of  a  plate  diminishes  with  increasing  fre¬ 
quency  (or  node  number);  see  Figure  1. 

Recent  measurements  made  by  Smith  et  al.4  or.  the  fundamental  and  higher  modes  of 
vibration  of  clamped  stiffened  plates  show  that  the  different  clamp  arrangements  used  did  not 
erfect  the  mode  shapes  but  did  affect  the  frequencies. 

Thus  to  obtain  a  reasonable  approximation  to  the  vibroacoustic  response  for  a  clamped 
plate,  we  need  merely  determine  the  frequencies  for  the  freely  vibrating  c’amped  plate  and  in¬ 
sert  these  predetermined  eigenvalues  as  input  dato  to  the  appropriate  programs  of  Reference  1. 

In  view  of  the  above,  we  seek  to  devise  optional  methods  (including  programs)  for  de¬ 
termining  the  frequencies  of  freely  vibrating  clamped  plates.  The  establishment  of  accurate 
methods  of  calculation  of  the  frequencies  for  all  modes  requires  comparing  the  theoretical 
frequencies  as  computed  by  various  methods  to  the  experimental  frequencies  and  using  the 


*The  modes  for  a  plate  are  usually  treated  in  terms  of  products  of  the  modes  for  a  beam  (see  Appendixes  A— G). 
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SIMPLY  SUPPORTED 


Figure  1  -  Examples  of  Mode  Shapes 

NOTE.  The  analysis  in  Reference  5  suggests  that  a  clamped  edge  panel  has  approximately  the  same  transverse 
vibrational  behavior  as  a  simply  supported  panel  whose  orthogonal  dimensions  and  bending  wavelengths  are  smaller 
1.0S  1.05 

by  the  ratios  £m  =  j  ^  q  g —  and  £n  -  g-g-  respectively;,  m,  n  are  node  numbers  (number  of  half  wavelengths  in 

the  plate  in  the  *•  and  y-coordinate  directions).  Here,  m-n.  Thus,  £m  and  <£,  can  be  termed,  “bending  wavelength 

a 

equivalency  factors.’’  The  physical  significance  of  these  ratios  is  clear  from  the  figure  where  —  =  £ :n . 
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results  of  this  comparison  to  select  the  best  methods.  The  modes  which  are  intrinsically 
associated  with  the  frequencies  can  also  be  computed  using  the  methods  or  programs  recom- 
meade d;  the  codes  may  be  of  value  to  users  interested  in  making  modal  comparisons  and  in 
applying  the  results  presented  here  to  other  problems. 


CALCULATION  AND  RESULTS 


Table  1  compares  computed  and  experimental  results  obtained  for  the  natural  frequen¬ 
cies  of  a  clamped-clasped  steel  plate.  The  methods  and  programs  used  in  the  computations 
are  respectively  described  in  Appendixes  A— H  and  Appendix  I. 

Toe  frequencies  versus  mode  numbers  gives  in  Table  la  for  each  method  are  plotted 
as  Figure  Sa.  The  frequencies  versus  method  given  in  Table  lb  for  each  mode  number  are 
plotted  as  Figure  2b.  Experimental  results  cited  by  Izzo  are  also  included  in  Table  la. 

Figure  3  compares  the  effect  of  clamped-clamped  and  simply  supported  boundaries  on 
the  vibratocy  response  of  a  plate  subject  to  turbulence  excitation.  The  results  were  obtained 
fay  osieg  the  Warburton  metfecd  for  computing  the  natural  frequencies  of  clamped-clamped 
plates  (see  Appendixes  A  and  I)  and  the  average  of  the  natural  frequencies  obtained  from 


£be  simple  frequency  expression  oj_c 


and  from  Warburtons 


method  fee  simply  supported  plates  in  the  Maestrello  program  for  vibratory  response.  Note 
that  the  computer  program  for  the  Warburton  method  given  in  Appendix  I,  yields  results  for 
hodfe  the  clamped-clamped  and  She  simply  supported  plates  (see  pages  97  and  103). 

Table  2  summarizes  key  features  associated  with  the  basic  references.  Some  of  these 
features  exceed  those  investigated  ia  this  paper.  They  may,  however,  be  of  interest  to  users 
ard  investigators  who  wish  to  extend  the  work  of  the  present  study. 


EVALUATION 

A  comparison  cf  the  computed  natural  frequencies  obtained  by  several  methods  (see 
Table  1  and  Figures  2a  and  2b)  shows  that  all  of  these  methods  yield  frequency  results  which 
are  in  good  agreement  with  each  other.  Hence  on  purely  theoretical  grounds,  any  method  can 
be  used  if  the  percentage  deviation  (obtained  from  the  results  of  Table  1)  between  the  mini¬ 
mum  (or  maximum)*  frequency  value  and  the  value  computed  by  the  specific  method  is  accept¬ 
able  for  a  particular  mode. 

However,  a  comparison  of  the  computed  and  experimental  natural  frequencies  given  in 
Table  la  and  Figures  2a  and  2b  as  well  as  an  appreciation  of  the  significant  features  involved 
in  carrying  out  a  computation  lead  to  a  preference  for  the  Warburton  method.  Using  Izzo’s  ex¬ 
perimental  results  as  a  standard  the  data  in  the  table  and  figures  show  that  for  the  modes 
treated,  the  maximum  error  attributable  to  the  Warburton  method  is  less  than  3.0  percent  for 


•The  deviation  from  the  minimum  or  maximum  is  taken  according  to  which  one  produces  the  greater  deviation 
for  a  particular  modal  frequency. 
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TABLE  1 

Comparison  of  Natural  Frequencies  Computed  by  Various  Methods  for  a  Clampod-CIampod  Steel  Plate 
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Table  la  —  Computed  Natural  Frequencies  for  Plate  1  (Izzo-Electric  Boat)  with 
Dimensions  2.0  x  2.33  x  0.0313  Feet  (see  Appendix  I) 
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986.4  983.6  982.9  986.6  986.9  1002  985.8 


D 

Wilily* 

(E<pcria«fltol) 

H earn on* 

Worburton 

Young** 

Bollentine* 
Plumb! ee 

Green spoe 

lliile^ 

Crocker 

W 

Ctoei»»-tt 

Thorne 

u 

541 

586 

577.9 

- 

581.0 

577.4 

581.1 

598.6 

- 

577.0 

1.2 

1307 

1439 

1402 

- 

1394 

1402 

1395 

1433 

- 

1398 

1,3 

2498 

2726 

2647 

- 

2636 

2647 

2646 

2484 

- 

2638 

2, 1 

833 

904 

9128 

- 

907.2 

9123 

9021 

941.0 

- 

923.3 

2,2 

1567 

1730 

1714 

- 

1703 

1714 

1741 

1758 

- 

1717 

2.3 

2747 

3010 

2954 

- 

2937 

2954 

2962 

3009 

- 

2948 

3.1 

1351 

1443 

1474 

- 

1465 

1473 

1503 

1502 

- 

1499 

3,2 

200S 

2241 

- 

2229 

2240 

2276 

2287 

- 

2259 

3,3 

- 

3488 

3461 

- 

3449 

3460 

3462 

3525 

- 

- 

4.1 

2047 

2186 

2247 

- 

2237 

2245 

- 

2273 

- 

2290 

4,2 

2646 

2939 

2986 

- 

2969 

2985 

- 

3030 

- 

3022 

*  Results  obtained  fro*  Reference  11.  Wilby’*  experimental  results  were  found  to  lie  between  tbe  simply  supported  and  fully  fixed 
edge  conditions  ;n  this  reference.  Hence,  com  per  Ison  between  .r.eory  and  experiment  is  of  limited  value. 

**Not  computed  for  tfiis  plate  but  computed  for  plate  in  Toble  la. 

*See  third  footnote  to  Table  la. 

**See  lost  footnote  to  Toble  la  (Izxo-  Wilby). 

Table  lb  —  Computed  Natural  Frequencies  for  Plate  2  (Wilby)  with 
Dimensions  4.C  x  2.75  X  0.015  Inches  (see  Appendix  I) 


n,  n 

Wilby* 

(Experimental) 

Heormon* 

Worburton 

Young** 

Bdlentine* 

Plunblee 

Green  spon 

White^ 

Docker 

$»«** 

Cloossen-^ 

Thorne 

1. 1 

1058 

925 

935.1 

- 

935.2 

934.6 

935.8 

954.9 

- 

935 

1.2 

2495 

2409 

2433 

- 

2439 

2432 

2435 

2464 

- 

2434 

1265 

1215 

1214 

- 

1211 

1214 

1236 

1249 

- 

1211 

2742 

2689 

2708 

- 

2706 

2709 

2781 

2756 

- 

2704 

1723 

172? 

1711 

- 

1704 

1711 

1731 

1751 

- 

1703 

3140 

3165 

3174 

- 

3173 

3175 

3332 

3230 

- 

3168 

4,  1 

2403 

2456 

2423 

- 

2411 

2423 

- 

2465 

- 

2409 

5,  1 

3321 

3392 

3341 

- 

3322 

3341 

- 

3382 

- 

- 

*See  first  footnote  to  Toble  lb. 

|  **Not  computed  for  this  plate  but  computed  for  plate  in  Table  la. 

^See  third  footnote  to  Toble  la. 

**$ee  last  footnote  to  Toble  la  (Iz20-*  Wilby) 

Table  lc  —  Computed  Natural  Frequencies  for  Plate  3  (Wilby)  with 
Dimensions  4.0  X  2.0  x  0.015  Inches  (see  Appendix  I) 
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Comparison  of  Theoretical  ami  Experimental  Natural  Frequencies 


CROCKER  CREENSPON 


IK  SQUARE  INCHES 


a -10  FT 
i- 0.54  Mi  FT 


Figure  3  —  Modal  Mean  Square  Plate  Displacement  for  Clamped-Clamped 
and  Simply  Supported  Aluminum  Plate 
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all  inodes.  Thus  it  is  acceptably  accurate  for  many  (probably  roost)  applications.  In  addition, 
the  Warburton  program  is  relatively  easy  to  run  on  a  computer  and  requires  little  running  time 
per  mode  (1.1  minutes  for  50  modal  frequencies  on  the  IBM  7090);  this  makes  for  a  relatively 
inexpensive  computation  for  each  frequency. 

The  error  of  3  percent  may  be  exceeded  for  square  plates  (see  Appendix  A),  and  hence 
an  alternative  method  of  computation  may  be  desirable  for  this  case. 

If  a  computer  is  not  available,  calculation  of  the  natural  frequencies  for  a  finite  rec¬ 
tangular  clamped-clamped  plate  can  be  performed  manually  by  any  of  several  methods  present¬ 
ed,  using  closed  form  analytical  or  nomographic  or  graphical  computations  (see  Appendixes 
A— F,  Appendix  H,  and  Table  2). 

The  frequencies  of  clamped-clamped  thin  plates  with  cylindrical  curvature  can  be  ob¬ 
tained  by  use  of  the  Ballentine-Plumblee  method. 

The  frequencies  of  thin  plates  with  clamped  and  rotational  supports  can  be  obtained 
by  use  of  the  White  method  (Appendix  E)  or  by  an  extension  of  the  Greenspon  method  (Appen¬ 
dix  D)  given  in  Reference  12. 

Figure  3  shows  that  at  the  convection  velocities  considered,  the  value  of  the  modal 
mean  square  displacement  for  any  mode  of  clamped  plates  subject  to  turbulence  excitation  is 
less  than  the' corresponding  value  for  simply  supported  plates.  The  difference  in  the  plate 
response  corresponding  to  the  two  boundary  conditions  increases  with  convection  velocity  for 
any  mode,  but  the  difference  is  relatively  constant  at  higher  convection  velocities  in  the 
region  of  maximum  response. 

The  nature  of  the  curves  in  Figure  3  suggests  that  at  low  convection  velocities  (1) c  - 
300  ft/sec),  the  difference  between  the  response  of  a  clamped-clamped  and  a  simply  supported 
plate  is  significantly  greater  for  the  lower  mode  (m,  n  =  5,1)  than  for  the  higher  mode  (m,  w  = 

7,  1).  It  appears  from  this  result  that  the  statement  previously  made,  namely,  that  the  effect 
of  the  boundary'  conditions  on  the  natural  frequencies  of  a  plate  diminishes  with  increasing 
frequency  (or  mode  number),  can  be  extended  to  include  a  diminishing  influence  of  boundaries 
on  the  higher  mode  response  to  turbulence  at  low  convection  velocities.  For  very  low  convec¬ 
tion  velocities,  the  trend  of  the  curves  suggests  that  the  concept  is  also  applicable  to  the 
lowest  modes. 

The  magnitude  of  the  curves  indicates  that  the  contribution  of  the  higher  mode  to  the 
total  response  is  not  negligible  for  either  boundary  condition,  i.e.,  the  contribution  of  the 
(7, 1)  mode  to  the  total  response  is  of  the  same  order  of  magnitude  as  that  of  the  (5, 1)  mode 
for  a  given  boundary  condition.  Thus,  determination  of  the  total  response  requires  that  the 
computations  include  the  contribution  of  the  several  modes  of  vibration  deemed  to  be 
significant. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  and  recommendations  are  based  on  the  results  of  the 
present  investigation. 

L  For  computing  the  vibroacoustic  response1  of  thin  clamped-clamped  rectangular 
plates,  the  modes  and  natural  frequencies  are  adequately  represented  when  the  modal  frequen¬ 
cies  are  calculated  by  considering  the  true  (clamped-clamped)  end  conditions  but  using  the 
mode  shapes  considering  the  end  conditions  to  be  simple  supports. 

2.  For  a  thin,  finite,  rectangular  clamped-clamped  plate,  the  Warburton  method  of  compu¬ 
tation  (including  computer  program)  of  the  natural  frequencies  is  acceptably  accurate.  For 
this  reason  as  well  as  for  its  relative  simplicity,  short  running  time,  and  inexpensiveness  in 
computer  application,  it  is  preferred  to  the  other  computer  methods. 

3.  If  a  computer  is  unavailable,  any  of  the  manual  methods  of  computation  presented  in 
Appendixes  A— F  and  H  can  be  used.  The  results  shown  in  Table  la  indicate  the  degree  of 
accuracy  to  be  expected  from  a  particular  method.  Moreover,  as  shown  in  the  tables  and  dis¬ 
cussed  in  the  Appendixes,  because  of  the  limited  data  available,  certain  methods  are  appli¬ 
cable  for  only  a  limited  range  of  mode  numbers. 

4.  For  clamped  thin  plates  with  cylindrical  curvature,  the  Ballentine-Plumblee  method 
(Appendix  C)  should  be  used  to  obtain  the  natural  frequencies. 

5.  For  thin  rectangular  plates  with  clamped  and  rotational  supports,  the  White  method 
(Appendix  E)  or  the  extension  of  the  Greenspon  method  (Appendix  D)  given  in  Reference  12 
should  be  used  to  obtain  the  natural  frequencies. 

6.  The  effect  of  the  boundary  conditions  on  the  natural  frequencies  of  a  plate  and  on  the 
response  of  a  plate  subject  to  turbulence  excitation  at  low  convection  velocities  diminishes 
with  increasing  frequency  (or  mode  number). 
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APPENDIX  A 

THE  WARBURTON  METHOD 


NOTATION 

A 

a,  h 

c,  k 
E 

fmn 

*«.  «*  K 

<V  Hy%  Jy 
9 
h 

m,  n 
T 
t 
U 
W 


w 


x,  y 

y,  € 


e,  $ 

A 

P 

a 

o 


Amplitude 

Length  and  width  of  sides  of  rectangular  plate  along 
#-  and  y-directions  respectively 

Ratios  in  expression  for  displacement 

Young’s  modulus 

Frequency,  modal  frequency 

Functions  of  m  in  frequency  expression 

F unctions  of  n  in  frequency  expression 

Acceleration  due  to  gravity 

Thickness  of  plate 

Mode  numbers  in  x-  and  y-directions,  respectively 

Kinetic  energy 

Time 

Potential  or  strain  energy 

Waveform  defined  by  Equation  (A2)  or  amplitude  of 
displacement  w,  i.e.,  w  =  W  sin  wt 

Transverse  displacement  of  a  point  on  the  plate 

Coordinate  distances  in  plane  of  plate 

Factors  in  amplitude  expression  defining  modal  pattern 

Functions  of  x  and  y,  respectively,  defining  waveform 

Nondimensional  frequency  factor  defined  by  Equation  (A8) 

Weight  per  unit  volume  of  plate 

Poisson’s  ratio 

Circular  frequency,  equal  to  2  vf 
lo 


-4 

■■  >1 


j 

$ 

i 
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■l 


♦ 


■i 
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DESCRIPTION 


Using  this  piste  theory,  Wsrbnrton13  derived  an  approximate  frequency  formulation  for 
all  modes  of  vibration  by  applying  the  Rayleigh  method  and  by  assuming  that  the  waveforms 
of  transversely  vibrating  rectangular  plates  and  beams  are  similar.  For  a  fully  clamped 
plate,  the  waveform  is  assumed  to  be  the  product  of  the  characteristic  functions  (discussed 
below)  foe  two  beams  with  fixed  ends.  The  plates  are  assumed  to  be  isotropic,  elastic,  free 
fcoar  applied  loads,  and  with  a  thickness  that  is  both  uniform  and  small  compared  to  the 
wavelength.  The  frequency-  is  expressed  in  terms  of  boundary-  conditions,  the  modal  pattern, 
she  dimensions  of  the  plate,  and  the  constants  of  the  material.  Because  of  the  imposition  of 
additional  constraints  oa  the  system  required  by  the  Rayleigh  method,  the  resulting  frequen¬ 
cies  are  higher  than  those  given  by  an  exact  analysis.  To  use  this  method,  the  modal  pat¬ 
terns  must  consist  of  lines  approximately  parallel  to  the  sides  of  the  plate.  This  require¬ 
ment  Is  satisfied  for  clamped  rectangular  plates,  and  the  errors  are  small.  The  exceptions 
and  their  effect  on  frequency  associated  with  some  modes  of  square  plates  are  discussed  in 
Reference  13. 

DERIVATION 

The  homogeneous  equation  for  a  freely  vibrating  thin  plate  is14 


0*t: 

d*tc 

-  a 

d*u?  !2p(l  —  a2) 

d2ic 

dx* 

dx2dy2 

By*  Egh 2 

dt 2 

(Al) 


The  solution  of  Equation  (Al)  is  assumed  to  have  the  form  of  a  product  of  separable 
solutions. 


tr(x,y,  Q  sin  at  =  A  0(x)  6(y)  sin  at  (A2) 

(The  nrotion  in  each  mode  is  vaa(x,y,  f)  =  B'ran  sin  amrt  =  Amn  0m(x)  6n{y)  sin  amnt  where 

the  actual  A  may  be  obtained  from  measurements.)  Here  0(z),  6(y),  the  characteristic 

beam  functions  or  mode  shapes  which  satisfy  the  boundary  conditions  for  plates  with  fixed 
dis  dw 

edges  (ir  =  —  =  0  at  x  =  0,  a  and  i c  =  — -  =  0  at  y  =  0,  b)  are  assumed  as  follows  ( m  and  n 

dx  ay 

are  node  numbers  and  correspond  to  m-  1  and  n-1  modes  respectively;  see  footnote  at  end 
of  ibis  Appendix). 


«>(*>.  cos  -  !)- 
6(i)  .  sin  r'fe  -  j)* 


k  cosh  y 


£'sinh  y 


(A3a) 


(A3b) 
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.  y 

sin  — 

2  y  y 

where*  k  =  -  and  tan  —  +  tanh  —  =  0  in  Equation  (A3a) 

sinh  — 

2 

y' 

sin  — 

2  y'  y' 

and  k'  = - -  and  tan - tar.h  —  =  0  in  Equation  (A3b). 

v  2  2 

sinh  — 

2 

The  corresponding  expressions  for  <f>(y)  are  obtained  by  substituting  y ,  b,  e,  and  c  for  x,  a, 
y,  and  k,  respectively. 

For  a  rectangular  plate,  the  potential  and  kinetic  energies  are  respectively  given  by15 


«t.f  f*  1 

J0  J0  12(1  -a2)  LW'  \dy2/  dx2  dy2  \dxdy  J  J 


(A3c) 


and  the  maximum  values  of  these  quantities  are 


”-r  r(‘[©‘*©' 


d2W  d2W 

dx2  dy2 


(d2W  \21 

[i~a)  W)  _p‘ 


i  phco2  r°  rfc 

=  2"  ~r  J  J 

•'0  Jo 


W2  dx  dy 


,  y  y  y  y 

*The  equation  tan  +  tanh  -g*  =  0  is  transcendental  and  may  be  solved  by  plotting  —  tanh  -rj-  and  tan  -tg 

and  looking  for  the  series  of  intersections.  Then  m  =  1  corresponds  to  the  value  of  y  for  the  first  intersection, 
m  =  3  for  the  second,  etc. 


I 


Equating  T  and  U  a?  required  by  the  Rayleigh  method,  we  have 


U_ 


,y2  = 


|ij  J 

J o  •'o 


(A7) 


By  the  Rayleigh  principle,  if  a  suitable  waveform  W  =  A  6(x)  6(y)  is  assumed  and 
approximately  satisfies  the  boundary  conditions,  the  resulting  frequency  value  is  slightly 
higher  than  the  true  value  because  the  assumption  of  an  incorrect  waveform  is  equivalent  to 
the  introduction  of  constraints  in  the  system. 

Substituting  the  expressions  for  the  characteristic  beam  functions  &x  and  tj>y  given  by 
Equations  (A3a)  and  (A3L)  which  satisfy  the  boundary'  conditions  for  the  clamped  plate,  into 
Equations  (A2)  and  (A7),  the  following  expression  for  the  approximate  frequency  is  obtained 


f  = 


F 


'Eh2, 


V 


4  n2Pa4  12(1  -  c2) 


(AS) 


where 


A2  =  ff  4  +  e4 

x  7 


2  a2 

V 


(A9) 


Here  coefficients  G%,  Gy,  Ux,  Hy,  and  Jy  depend  on  the  modal  pattern  and  boundary 
conditions.*  Values  of  these  coefficients  are 


G  J1-056 

*  \m  -  1/2 


#*  =  •4 


Hy-Jy 


for  m  =  1 

for  m  =  2, 3,  4  . . . 

for  n  =  1 

for  n  =  2, 3, 4  . . . 

for  m  =  1 


for  m  =  2,  3, 4,  . . . 
for  n  =1 

for  n  =  2, 3, 4,  . . . 


*In  Reference  13,  m  refers  to  the  number  of  nodes  along  the  plate  length  and  hence  to  m  -  1  modes.  In  the 
present  paper,  however,  m  refers  to  the  mcde  number.  The  letter  notation  is  more  common  and  is  consistent  with 
the  notation  used  by  Maestrello  and  other  investigators.  This  definition  for  •n  is  now  reflected  in  the  numerical 
values  of  m  used  in  computing  the  coefficients  G^,  H x,  ]  whereas  the  values  for  m  used  previcusiy  (Equations 

(A3a)  and  (A3b))  correspond  to  the  Warburton  definition  in  Reference  13.  A  similar  situation  holds  for  n. 
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Hence  for  a  given  m,  n  mode  and  —  ratio,  we  obtain  the  appropriate  value  of  the  co- 

b 

efficients  for  use  in  determining  A2  from  Equation  (A9).  For  a  given  ratio  a/b,  the  corre¬ 
sponding  approximate  frequency  is  found  from  Equation  (A8)  to  be 


Utt  r  Eg  "I 

a2  L48p(l-a2)J 


(A10) 


For  mode  numbers  mn,  As  and  / sL,  and  a  =  The  corresponding  mode 

5  mn  1  1 mn  mn  1 mn  *  ° 

shape  is  then  Wmn  =  Amn  6m  ( » )  <j>n  (y). 
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APPENDIX  B 


THE  YOUNG  METHOD 


Coefficient  used  in  series  representation  of  deflection 

Length  and  width  of  plate  along  x -  and  y-directions 
respectively  ’ 

Coefficients 

Bending  stiffness  of  a  plate  equal  to  Eh3/12(l  -  p2) 
Modulus  of  elasticity 

Definite  integrals 

Frequency 
Poisson’s  ratio 
Thickness  of  plate 

Positive  integers 
Length  of  beam 

Elastic  strain  energy  of  bending  of  a  plate 

Lateral  deflection  of  plate 

Function  of  x  alone 

Rectangular  coordinates 

Function  of  y  alone 

Parameter  in  expressions  for  <f>T 

Kronecker  delta 

Parameter  in  expressions  for  <f>T 
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Characteristic  value  equal  to 


o2 pha3 b 


D 


Poisson’s  ratio 

Mass  density  of  plate  material 

Characteristic  function  of  a  vibrating  beam 

Angular  frequency  equal  to  2s  f 
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DESCRIPTION 

Young16  uses  the  Ritz  method  to  obtain  approximate  solutions  for  the  freqnencies  and 
modes  of  vibration  of  thin,  homogeneous  plates  of  uniform  thickness;  the  freqnencies  calcu¬ 
lated  by  the  Ritz  procedure  are  always  higher  than  the  exact  values.  To  represent  the  plate 
deflection.  Young  treats  combinations  of  the  characteristic  functions  which  define  the  normal 
modes  of  vibration  for  a  uniform  beam.  He  computes  and  tabulates  values  of  these  r unctions 
as  well  as  associated  integrals  and  derivatives  of  the  functions.  W»th  me  aid  of  these  tables, 
the  user  can  set  up  and  solve  the  necessary  equations  with  reasonable  effort.  A  simple  iter¬ 
ation  procedure  is  used  to  solve  the  equations. 


DERIVATION 

The  maximum  potential  and  kinetic  energies  for  a  harmonically  vibrating  uniform  plate 
are,  respectively  (see  Appendix  A), 


V  = 


d2« 

-  +  2(1  -  ft) 

dy2 


dxdy 


(Bla) 


T  = 


(Bib) 


Equating  these  expressions,  we  obtain 


a 


2 


2_  V 
// w2  dxdy 


(B2) 


The  Ritz  method  consists  of  assuming  the  deflection  w(x,  y)  as  a  linear  series  of 
“admissible”  functions  and  adjusting  the  coefficients  in  the  series  so  as  to  minimize  Equa¬ 
tion  (B2).  For  rectangular  plates  with  edges  parallel  to  the  z-  and  y-axes,  Young  represents 
w  by  the  following  approximate  series: 


»(*,,).  I  1  Am,  *<*>  r<’> 


m  =  1  n  =  1 


(B3) 


Each  function  Xm  Yn  must  be  admissible,  i.e.,  it  must  satisfy  the  so-called  artificial  bound¬ 
ary  conditions  which  are  the  prescribed  values  for  the  deflection  and  -for  the  slope.  It  need 
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not  satisfy  any  natural  boundary  conditions  which  require  that  second  or  third  derivatives  or 
combinations  thereof  vanish  at  the  boundary.  Satisfaction  of  these  latter  conditions,  if  pos¬ 
sible,  is  desirable  however  in  accordance  with  practical  consideration  of  the  rate  of 
convergence. 

Substituting  for  »o  (x,  y)  in  Equation  (B2)  using  Equation  (B3)  and  minimizing  the 
right-hand  side  by  taking  the  partial  derivative  with  respect  to  each  coefficient  Amm  and 
equating  to  zero,  we  obtain  a  set  of  linear  homogeneous  equations  in  the  unknown  Amn  each 
of  which  Isas  the  form 


cP’  ph  d 


ffv>2  dxdy  *0 


where  Aik  is  any  one  of  the  coefficients  Amn.  The  natural  frequencies  &v  o2  are  determined 
from  the  condition  that  the  determinant  of  the  system  must  vanish. 

For  a  clamped-clam-ptd  beam ,  the  infinite  set  of  characteristic  functions  is  given  by 


d,  =  cosh - cos 


— - ar  ^sinh  -  -sin  —  J - '-1,8,! 


0  -*  -£ 


(The  method  for  determining  the  set  of  characteristic  functions  which  define  the  normal  modes 
is  gjven  in  References  15  and  17.) 

The  numerical  values  of  «r  and  er  for  each  set  of  functions  is  given  in  Table  3.  Ref¬ 
erence  8  tabulates  values  of  these  functions  to  five  decimal  places  at  intervals  of  the 
x 

argument  —  =  0.02. 

The  function  $r  given  by  Equation  (B5)  satisfies  both  the  boundary  (i.e.,  end)  condi- 

d$T 

tions  for  the  clamped-clamped  beam  6r  =  — —  =  0  at  x  -  0,  l  and  the  differential  equation  for 

a 

d**r  <r*r 

the  beam  -  =  -  .  Also  any  set  of  functions  4>r  and  <3  are  orthogonal  for 

dx*  E4 

0  $x$t,  i.e.. 


J  * 

—  n 


4>  dx  =  1  (for  t  =  s') 


=  0  (for 


'-  *)  1 
r  +  a)  j 


The  second  derivatives  of  the  functions  of  the  set  are  also  orthogonal  and  satisfy  the 
relations 


i 


(for  T  =  «) 


r 


}.  • 

i; 


i 


i 
>.  * 

l- 

i 


I 


f  ^2<5r  d2of  < 


dx2  dz2 


=  0 


(for  t  ±  s) 


(B7) 


Numerical  values  of  are  given  in  Table  3.  In  addition  to  the  integrals  defined  by  Equa¬ 
tions  (B6)  and  (B7),  tie  Ritz  method  also  requires  evaluation  of  tie  integrals 


t 


d26 


Cl  do,  do. 


ds2 


■  s  I  -s  -  s 

—  dx  and  I  -  -  dx 


J  dx  dx 

Jo 


Table  4  gives  the  values  of  these  integrals  computed  by  Young. 

The  characteristic  functions  are  those  that  are  used  for  XB  and  Yc  in  Equation  (B3). 
Consider  a  rectangular  plate  bounded  by  the  lines  z  =  Q,  x  =  a,y  =  0,  y=b.  When  the  func¬ 
tion  is  used  for  A=,  we  take  £  =  a;  if  used  for  Yr,  we  take  1  =  b  and  replace  x  by  y.  Appro¬ 
priate  changes  of  the  subscripts  r  and  s  to  either  m  and  i  or  ton  and  k  are  to  be  made  in  the 
set  of  functions. 

It  is  convenient  Co  introduce  the  following  notation: 


=  a 


j 

•'A 


d2  X. 


dx 2 


dx. 


~f 

•*n 


X. 


rb  d2  Y  rb 

Fkn=>  \  —  dy,  Fak  =  b  Yn 

J  n  dy  Jn 


dx2 


d2  Y, 


dx 


dy 2 


dy 


r*  dxi  dx .  r 

•--‘J  -3T-*r'fe’  '-•*1 

J0  Jc 


b  dJk  d\ 

dy  dy 


dy 


(BS) 


(89) 


(BIO) 


Since  the  appropriate  ^-functions  are  to  be  used  for  Xm  and  Yn,  the  numerical  value  of  these 
integrals  can  be  obtained  directly  from  the  data  given  in  Table  4. 

From  Equations  (Bla)  and  (B3)  and  the  orthogonality  relations  (Equations  (B6)  and 
(B7)),  the  set  of  Equations  (B4)  can  be  reduced  to  the  form 


£  1  lC<£*>  -  A 5m  J  =  0 

_ ,  rnn  n  rnn 

m  =  1  n  =  1 


ttiti  mn 


(Bll) 
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Type  of 
Beam 


Clamped 

Clamped 


TABLE  3 
Values  of  or.  and  f 


TABLE  4 

Integrals  of  Characteristic  Functions  of  Clamped-Clamped  Beam 
r?  do.  d6_ 


C  d*a 

Values  of  l  1  — - —  — —  dx 


V 


1 

2 

3 

1230262 

0 

-  9.73079 

0 

46.05012 

0 

-  9.73079 

0 

98.90480 

0 

-17.12892 

0 

-  7.61544 

0 

-24.34987 

0 

-15.19457 

0 

-  17.12892 

0 

171.58566 

0 

-  31.27645 


f  <P6a  r  d<f>r  d4>a 

E:  |  <f>r  -  dx  -  -  I  -  -  dx 

L  dx 2  d*  d* 


5 

6 

-  7.61544 

0 

0 

-  15.19457 

-  24.34987 

0 

0 

-  31.27645 

263.99798 

0 

0 

376.15008 

where 


A  =  - 


!  ph  a3  b 

~D 


(B12) 


and 


S  =1  for  mn  =  ik 
mn 

=  0  for  mn  £  ik 


*  2(1  -c)  f  (BIS) 


which  is  valid  for  mn  ^  For  mn  =  z'£,  the  coefficient  is 


7  ,4  +  ^  */  +  *!■  f  VH  +  S(l-rt-j.*„jrw  (B14) 


In  Equation  (B14),  is  to  be  taken  from  the  data  in  Table  3  corresponding  to  the  ^-function 
that  represents  -Am,  whereas  is  to  be  taken  from  data  for  the  ^-function  that  represents 

y„. 

There  will  be  one  equation  of  the  type  (Bll)  for  each  of  the  q  combinations  of  ik. 

In  general,*  an  iterative  procedure18  is  used  to  find  the  characteristic  values  of  A  from  the 
condition  that  the  determinant  of  this  system  of  equations  must  vanish.  Results  for  a 
clamped  square  plate  are  given  in  Reference  16. 


*A  manual  computation  can  be  performed  for  systems  with  no  more  than  three  or  four  equations. 
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NOTATION 


APPENDIX  C 


THE  BALLENTINE-PLUMBLEE  METHOD 


Simple  panel  aspect  ratio;  ratio  of  arc  length  to  straight 
edge  length 

Midplane  radius  of  simple  panel 
Panel  arc  length 

Young’s  modulus  for  isotropic  material 

Simple  panel  thickness 

Panel  length  (for  simple  and  sandwich  panel) 


Generalized  coordinate 


Kinetic  energy 

Length  to  thickness  ratio  for  simple  panel 

Strain  energy 

Generalized  coordinate 

Strain  energy  density 

Midplane  displacement  in  ^-direction 


Generalized  coordinate 


*.(*> 


Midplane  displacement  in  y-direction 
Midplane  displacement  in  radial,  a-direction 
Mode  shape  for  ^-coordinate 
Shell  midplane  coordinate 
Mode  shape  for  y-coordinate 
Snell  midplane  coordinate,  y  =  a<j> 

Shell  midplane  coordinate  through  thickness 
Constant  appearing  in  clamped  mode  function 


Preceding  page  blank 


L 


0= 

7- 


€- 


l  1 
1  I 
l  1 

m 


Constant  appearing  in  mode  function 
Constant  appearing  in  mode  function 
Strain 

Constant  appearing  in  clamped  mode  function 
Nondimensicnal  frequency 
Poisson’s  ratio  for  isotropic  material 
Mass  density 
Stress 

Angle  which  defines  cylindrical  coordinate  y 
(generalized  coordinate) 

Circular  frequency 

Row  matrix 

Column  matrix 

Rectangular  matrix 

Diagonal  mrtrix 
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DESCRIPTION 

Ballentine19  uses  the  Rayleigh-Ritz  energy  method  for  finding  the  frequencies  and 
normal  modes  of  a  cylindrically  curved  panel  with  clamped  edge  conditions*;  the  results  in¬ 
clude  those  for  the  flat  plate.  For  clamped  edges,  inexact  mode  functions  which  satisfy  only 
the  geometric  boundary  but  not  the  differential  equations  are  used.  The  analysis  assumes 
that  the  material  is  linearly  elastic  and  orthotropic  and  that  the  panel  thickness  is  much  less 
than  the  major  panel  dimensions,  i.e.,  the  elasticity  theory  of  thin  shells  is  applicable.  Only 
the  main  analytical  steps  and  chief  results  are  discussed  here.  The  reader  interested  in 
studying  the  associated  details  of  matrix  manipulation  is  referred  to  Reference  19. 


DERIVATION 


The  total  strain  energy  U  of  the  curved  plate  (Figure  4)  obtained  by  integrating  the 
strain  energy  density  U0  over  the  volume  of  the  plate 
is 


where 


(Cl) 


(C2) 


a.  is  expressed  in  terms  of  strain  ei  and  ;hen  the  strain  in  terms  of  displacements  which  are 
represented  by 


*2  7T  °mnWW 

Pm 

*-  Ssf  VmnXm(x)Yn'(y) 
Yn 

»  "  *.<*)  r.M 


CCS) 


V 


‘Results  for  simply  supported  conditions  are  also  presented  in  this  reference. 
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which  can  be  expressed  in  matrix  form.  The  boundary  conditions  for  a  curved  plate  with 
clamped  edges  are 

«(0,  y)  =  to(f, y)  =  v(x, 0)  =  w(x,  b)  =  0  ^ 

(0,  y)  =  wx  (£,  y)  =  wy  (a-,  0)  =  wy  {x,  b)  =  0  I 

v  (0,  y)  =  v(t,  y)  =  v(x,  0)  =  v(x,  b)  =  0  I 

u(0,  y)  =  u(£.  y)  =  u( x,  0)  =  u(x,  b)  =  0  ' 


The  assumed  mode  shapes  for  a  plate  with  clamped  edges  are 

Xm(x)  =  Cosh  fimx  -  Cos  Pmx  ~<*m  (Sinh  fim  x  -  sin  pmx) 
Yn(?/)  =  Cosh  yny  -Cos  yny  -dn  (Sinhyny  -  sin  yny) 


where 


Cosh  0m£  -  cos 
*“  ~  Sinh  pm  £  -  sin  pj 

Cosh  y„b  -  cos  yn  b 

6  =  - - 

"  Sinh  ynb  -sin  yn  b 


and  Pm  and  yn  are  determined  from 


Cosh  Pn.  £  cos  Pm  £  »  1 
Cosh  vn  b  cos  va  b  =  1 


Tne  kinetic  energy  of  the  vibrating  plate  obtained  by  integrating  the  product  of  mass 
and  one-half  velocity  squared  over  the  volume  of  the  plate  is 


T=  - 


b  1 

J  J 


(k2  +  v2  +  w2)  dz  dy  dx 


where  w,  v,  w  can  be  expressed  in  matrix  form  using  Equation  (C3). 

V  and  T  are  now  substituted  in  the  Lagrange  equation  of  motion  to  obtain  an  equation 
for  the  natural  modes  of  vibration  which  can  be  written  in  the  form 


[m  ~  <y2  ph  [</]]  \qr  1  =  0 


where  the  terms  in  the  [K]  and  [7]  matrices  are  given  in  Reference  19.  This  equation  can  be 
solved  for  the  modal  frequencies. 
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Reference  19  indicates  that  inasmuch  as  the  integrals  of  X'  X^,  X'' Xm  and  Xp  X" 
(which  were  used  in  deriving  the  terms  in  [/if]  [«/])  for  clamped  edge  conditions  are  nonzero 
when  p  4  m  then  the  analysis  does  not  display  the  desired  orthogonality  between  the  modes. 
However,  a  numerical  analysis  for  one  of  the  test  panels  used  in  the  reference  program 
showed  insignificant  differences  when  compared  to  a  numerical  analysis  which  assumed 
orthogonality.  A  complete  investigation  of  the  effects  of  including  this  non  orthogonality 
relationship  has  not  been  evaluated  because  of  computer  time  requirements.  Finally  a  simpli¬ 
fication  of  considerable  interest  to  the  orthotropic  curved  panel  frequency  analysis  occurs, 
provided  the  modal  integrations  are  taken  to  be  orthogonal  and  the  material  is  isotropic.  In 
this  case  the  modes  are  uncoupled,  and  assuming  that 


h2 

—  <  <  1 


(C9) 


we  find  that  the  determinant  of  the  coefficients  is 


and 


[[<;]  -  a2  [t]  |  ■  o 

r  ,  Eh 3 

[tf]=. -  [Q] 

l2(  l-v2) 

[J]  =tb[L] 

x2  m  f£i£ z*  *2 

Eh2 


(CIO) 


where  the  terms  in  [6]  and  [L\  are  given  in  Reference  19.  Equation  (CIO)  can  be  solved  for 
the  modal  frequencies. 

b 

If  a  =  oo  (flat  plate,  <f>  =  —  =0),  then  the  3x3  matrix  is  reduced  to  a  2x2  matrix  and  one 

a 

equation  in  terms  of  A2  in  the  3, 3  position.  The  equation  resulting  from  the  3,  3  element 
yields  the  flat  plate  flexural  modes,  whereas  the  2x2  matrix  gives  the  in-plane  or  longitu¬ 
dinal  vibration  modes. 

Some  important  simplifications  can  be  made  in  the  frequency  theory  if  the  angle  which 
the  panel  subtends  is  small.  For  angles  4>  less  than  0.2  radians,  the  frequency  of  flexural 
vibration  can  be  approximated  by  the  following  equation  when  all  edges  are  clamped: 


-  a  25.2  41.7 

A2  =  41.7  A  +  -  +  -  + 

A  A3 


t2  <£2 


;  (fj  <  0.2  radians 


(CU) 


where  A  = 


b 


l  ’ 


and  t  =  —  . 
h 
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It  follows  from  the  foregoing  equations  that  the  ratio  of  the  curved  panel  frequency  to  that  of 
the  infinite  panel  for  the  1, 1  mode  of  vibration  is 


/anc  \2  =i+  C(At6)2 

\°ii»  '  A*  +  0.61  A2  +  1 


where  the  theoretical  value  of  C  is  0.024  for  clamped  edges. 

The  frequency  analysis  for  isotropic  curved  panels  with  no  coupled  modes ,  Equation 
(CIO),  has  been  programmed  in  Fortran  language  for  solution  on  the  IBM  360/91  at  the 
Applied  Physics  Laboratory  of  Johns  Hopkins  University.  The  equations  are  nondimension- 
alized  in  terms  of  three  independent  variables  A,  t  and  the  dependent  variable  which  is 
nondimensional  frequency.  Calculation  of  the  frequency  for  clamped  plates  was  made  over 
the  following  range  of  variables: 


0  ^ 


b  < 

-  =  <f>  ^3.14 
a 


20  ^ 
0.5  ^ 


e_ 

h 

6 

e 


=  t  ^  1000 

=  A  -  2.0 


For  particular  values  of  aspect  ratio  A,  nondimensional  frequency  is  plotted  for  six 
modes  and  six  values  of  length-to-thickness  ratio.  Figures  5  to  9  give  clamped  edge  frequen¬ 
cies.*  Once  nondimensional  frequency  is  found,  the  actual  frequency  can  be  determined  from 
the  nomogram  shown  in  Figure  10. 

As  an  example,  the  natural  frequencies  of  a  clamped,  curved  panel  calculated  in  Refer¬ 
ence  19  are  presented.  The  panel  dimensions  are 

Radius  a  =  100  in. 

Arc  length  b  =  10  in. 

Length,  £  =  20  in. 

Thickness  h  =  0.05  in. 

The  nondimensional  ratios  are: 

A  =  0.5 
<f>  =0.1 
t  =  400 


‘Similar  results  are  presented  in  Reference  19  for  simply  supported  edges. 
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Code 

Mode  Number 

Ode 

Subtended 

Angle,  4 

Straight 
Edge,  m 

Curved 
Edge,  n 

c 

1 

1 

mmm 

0 

d 

1 

2 

fBEX 

e 

1 

3 

K3 

mBMmm 

f 

2 

1 

■3 

g 

2 

2 

'  Q 

1.0 

h 

3 

1 

m 

3.14 

1.0 

3.14 


Mode  Number 


Code 

Straight 
Edge, m 

Curved 
Edge,  n 

Code 

Subtended 
Angie,  j/ 

c 

■ 

■n 

mta 

0 

d 

KS 

0.05 

e 

vQ 

0.1 

f 

a 

'-El 

0.2 

g 

1.0 

h 

1 

H 

3.14 

Table  5  shows  values  of  A  for  the  different  combinations  of  mode  number.  These  values  were 
taken  from  Figure  5  for  A  =  0.50.  The  frequencies  converted  through  the  use  of  the  nomogram 
are  also  displayed  in  Table  5. 


TABLE  5 

Natural  Frequencies  for  Sample  Problem 


APPENDIX  D 


THE  GREENSPON  METHOD 


Area  of  plate 
Width  of  plate 
Length  of  plate 
Aspect  ratio 


Plate  modulus  = 


Eh 3 


12(1  -  v2) 
Differential  element  of  area 
Modulus  of  elasticity  of  plate  material 
Thickness  of  plate 


Distance  in  direction  normal  to  boundary  of  a  flat  plate  of 
arbitrary  shape  (has  dimensions  of  length);  n  lies  in  plane 
of  plate 


Circular  frequency  of  rth  mode  of  vibration 
Circular  frequency  of  ij  th  mode  of  vibration 


A  function  of  time  such  that  w  =  wm  qm  satisfies  the  homo- 


a2 


w 


geneous  plate  equation  D\'*w  ~  ph  -  =  0 

di2 


Distance  in  Jirection  of  boundary  of  a  flat  plate  of  arbitrary 
shape  (has  dimensions  of  length) 


Time  variable 
Lateral  deflection 

Deflection  function  in  rth  mode  of  vibration 
Normal  mode  functions  for  the  modes  of  vibration  of  a  beam 
Factors  defining  modes  of  vibration  cf  a  beam 
Frequency  numbers  of  the  modes  of  vibration  of  a  beam 
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Poisson’s  ratio 

Mass  per  unit  volume  of  plate  material 


Differential  operator  ( •——  +  2  d4 

dx 2  (?V2 

rectangular  coordinates 
Slope  of  plate  boundary 


L 
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DESCRIPTION 

Using  the  general  theory  of  small  vibrations  of  plates,  Greenspon7,  12,  20  presents  a 
method  for  calculating  the  frequency  and  deflection  response  of  a  clamped  rectangular  plate.* 
The  calculation  is  based  on  a  knowledge  of  the  normal  modes  of  vibrations  which  are  approxi¬ 
mated  by  the  product  of  two  beam  functions  (or  characteristic  functions)  identical  to  that 
used  by  Young  (see  Appendix  B). 


DERIVATION 

The  homogeneous  equation  for  a  freely  vibrating  thin  plate  is  given  by7,  12,  20 


For  a  clamped  boundary 


£>V4to  +  ph 


w  =0  along  s 
dw 

— —  =  0  along  s 
dn 


(Dl) 


(D2) 


The  deflection  of  the  plate  is  taken  to  be  the  sum  of  the  normal  modes. 

w(x,y,t)=  wr(x,y)  q{t)  (D3) 

Substitution  of  Equation  (D2)  into  Equation  (Dl)  yields 
D  A  oo  co  d2(lr 

—  V4  2  war~  2  to  -  =0  (D4) 

ph  r=l  r=  1  dt2 

Integration  of  the  product  of  Equation  (A3)  and  one  of  the  normal  mode  functions  wm 
over  the  plate  area  A  gives 


(Do) 


i 


•Results  for  isotropic  plates  are  gives  ia  Reference  7  and  for  otho  tropin  (i.e.,  stiffened)  plates  in  References 
12  and  20. 
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As  shown  in  Reference  12,  the  first  term  in  this  equation  which  contains  integrals  of 
the  form  / wm  wr  dA  is  zero  if  r  dm  and  the  second  term  in  this  equation  which  contains 
integrals  of  the  form  /  w  wrdA  is  also  zero  if  r  dm  and  the  plate  is  clamped.  Thus  if  the 
AP 

plate  is  vibrating  freely  in  one  of  its  modes  w=wf  sin  pTt ,  Equation  (Do)  can  be  written 


7T  /  ** w’dA " p'  S 


to  wr  dA 

rn  r 


(D6) 


and  since  the  integrals  have  a  value  only  for  r=m,  the  circular  frequency  of  the  mth  mode  of 
vibration  is 


P»  = 


If  Wm  '^4  Wn,  dA 

/  w2  dA 
. A  m 
P 


(D7) 


To  calculate  the  frequency  and  deflection  response,  the  normal  modes  of  the  clamped 
plate  are  approximated  by  the  product  of  two  beam  (or  characteristic)  functions,  i.e., 
wm  =  X{  Yj,  which  depend  on  the  boundary  conditions  of  the  plate. T  (For  the  first  mode 
i  =  1,  j  =  l;  for  the  second  mode  i  =  1,  /  =  2,  etc.)  (For  the  clamped  plate,  the  values  of  Xi 
and  Y-  used  by  Greenspon  are  identical  to  those  used  by  Young;  see  Appendix  B.) 


Pix 

X.  =  cosh -  -  cos  - 

a  a 


Yj  =  cosh 


/  Pi*  Pi*\ 

-  «-  (sinh  -  -sin  - / 

\  a  a  / 

Pjy  PjV  (  ftp  PfA 

— —  -cos  — —  -  «.  Isinh  — —  -sin  — — ) 
0  b  >  \  o  b  / 


(D8) 


Substituting  the  value  of  wm  =  Yj  into  Equation  (D7)  using  Equation  (D8),  we  find 
(see  page  30  of  Reference  12  for  details). 


Pir 


jrJw4  ,  2fo °£x‘xrYiYi 

V  Ph  7  a4  b 4  f°  fb  X?  Y-2  dx 

o  o  1  1 


"dx  dy 

dy 


(09) 


d2  X.  d2  Yj 

where  X"  =  - —  and  Yf'  - - - 


dx2 


dy2 


"The  product  of  the  bean  functions  is  not  an  exact  expression  for  the  nodes  of  a  damped  plate  because  it 
generally  does  not  satisfy  the  plate  equation. 
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The  values  of  /J  and  a  as  well  as  the  integrals  /“  X^-'dx,  f°  X f  dx  and  the  values 

of  X{  and  X"  which  are  contained  in  References  8,  9,  and  16  were  used  by  Greenspon7  to 
compute  Table  6. 

For  purposes  of  the  present  report,  the  final  expression  for  the  deflection  response 
derived  in  References  7,  12,  and  20  is  omitted  here. 

Following  jf  similar  procedure,  Reference  12  presents  a  frequency  equation  for  a 
fluid-loaded ,  cross-stiffened  plate,  i.e.,  orthotropic  plate.  It  also  gives  the  procedure  for 
determining  the  orthotropic  constants  and  other  data.  The  beam  functions  A';  Yj  are  written 
for  a  beam  with  rotational  constraint  which  includes  simply  supported  and  clamped 
constraints.  Thus  Equation  (D9)  is  a  special  case  of  the  more  general  frequency  equation 
given  in  this  reference. 


TABLE  6 

Function  Values  for  a  Clamped-Clamped  Beam 

(Here  a  or  b  is  the  length  of  the  beam,  and  the  origin  x  =  0  is  located  at  one  end. 
The  tabulations  will  remain  valid  if  X-  is  replaced  by  Y-.) 
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1 

P,P, 
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lb  Yfdy 

Value  of 

Point  at  Which 
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.V(  Occurs 
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a 

a 
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4.7300 
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i 

1.5882 

2 

n 

0.8X9 

2 

1.0008 
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0 

3 
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10.9956 
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NOTATION 


APPENDIX  E 


THE  WHITE  METHOD 


a  ,  a 
mn 5  rs 


Beam  width 

Coefficients  used  in  series  representation  of  deflection 

A  constant  which  determines  the  amplitude  of  response 
for  the  mth  and  nth  modes  respectively  of  a  beam;  beam 
nondimensional  frequency  parameters 

Beam  length 

Rotational  spring  stiffness  per  unit  length  along  the  ith 
edge 

Quantity  defined  by  Equation  (El6) 


Plate  bending  stiffness  equal  to  El  = 


12(1  -v2) 


Young’s  modulus  of  elasticity 
Gravity  acceleration 
Plate  thickness 


m ,  n  and 
r,s 


Moment  of  inertia  of  cross  section  of  the  beam  about 
the  neutral  axis 

Mass  moment  of  inertia  per  unit  length  along  the  zth 
edge 

Plate  mass 

Mode  numbers,  i.e.,  number  of  elastic  half-waves  parallel 
to  the  x-  and  y-axes,  respectively 

Edge  mass  per  unit  of  length  along  the  i th  edge 

Kinetic  energy 


Equal  to  -  :  defined  by  Equation  (E7) 


Potential  energy 


Preceding  page  blank 


1 


anO’«nL 


0<r)  0(.y) 

n  >  s 

A,  A 

5  mrc 


ft 

p 

*„(?)>  *,(sr) 


<y,  o>mF> 
7  m  jj 


Equal  to 


2P53 


;  defined  by  Equation  (Ell) 


Da 

Plate  deflection 

Rectangular  coordinates 

Beam  nondimensional  frequency  parameters 

Plate  nondimensional  frequency  parameters 

Nondimensional  frequency  parameters  for  the  rath  mode 
of  a  symmetrically  constrained  beam  which  has  springs 
of  stiffness  CQ  and  C respectively,  at  both  ends  of 
the  beam 

Defined  by  Equation  (E17) 

Beam  mode  shapes  (functions  of  y  only) 

Nondimensional  plate  frequency  parameters  defined  by 
Equations  (E13)  and  (E19),  respectively 

Plate  mass  per  unit  of  area 

Poisson’s  ratio 

Nondimensional  rotational  stiffness  parameter 
Mass  density 

Beam  mode  shapes  (functions  of  x  or  y  only) 

Plate  mode  shape,  approximately  equal  to  <£m(z)  <?n(y) 

Beam  functions  defined  by  Equation  (E19) 

Circular  frequency  and  circular  resonance  frequency 
of  plate,  respectively 

Designates  a  nondimensional  integral 
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DESCRIPTION 


i 

i 


Using  the  Rayleigh-Ritz  technique,  White21  derives  a  set  of  simultaneous  algebraic 
equations  for  computing  the  resonance  frequencies  and  modes  of  a  rectangular  flat  plate  hav¬ 
ing  a  uniform  distribution  of  elastic  and  inertial  edge  fixities.  These  fixities  are  equivalent 
to  a  uniform  attribution  of  independent  masses,  translational  springs,  and  rotational  springs 
along  each  edge  of  the  plate;  the  various  edges  of  the  plate  can  have  equal  or  different  elas¬ 
tic  constraints  and  inertial  loadings.  The  only  coupling  between  the  individual  masses  along 
an  edge  is  the  coupling  provided  by  the  deflection  of  the  plate.  Certain  integrals  of  products 
of  beam  mode  shapes  and  derivatives  of  these  mode  shapes  are  expanded  in  terms  of  modal 
displacements  and  derivatives  of  these  displacements  at  the  ends  of  the  beam.  These  inte¬ 
grals  are  used  to  develop  expressions  for  plate  frequencies.  All  effects  of  rotary  inertia  and 
shear  deformation  of  the  beam  are  neglected. 

Once  the  masses  and  springs  along  the  four  edges  of  the  plate  are  known,  the  frequen¬ 
cies  and  modes  can  be  numerically  evaluated.  Solutions  of  the  simultaneous  set  of  algebraic 
equations  can  be  obtained  by  iteration  using  standard  digital  computer  techniques. 

Reference  21  treats  the  special  case  in  which  the  edges  of  the  plate  are  translation- 
ally  fixed,  elastically  constrained  in  rotation  by  a  uniform  distribution  of  rotational  springs, 
and  not  loaded  by  edge  masses.  In  this  special  case,  each  edge  of  the  plate  can  have  a  fix¬ 
ity  arbitrarily  between  a  pinned  and  clamped  support  and  the  four  edges  can  have  different 
elastic  constraints.  The  special  case  is  further  specialized  in  the  present  report  to  treat 
only  the  completely  clamped  case.  Although  exact  solutions  of  the  corresponding  set  of  si¬ 
multaneous  frequency  equations  require  an  iteration  of  the  Ritz  type,  it  was  found  that  rea¬ 
sonably  accurate  estimates  of  the  plate  resonance  frequencies  c«n  be  obtained  by  using  a 
single  term  from  the  appropriate  equation  in  the  set.  The  resulting  approximate  frequency 
equation  is  given  as  well  as  nomographs  for  quick  computation  of  these  frequencies.1'  The 
White  method  as  applied  to  the  completely  clamped  plate  follows. 

DERIVATION 

The  partial  differential  equation  which  defines  the  undamped  resonant  vibration  of  a 
thin,  uniform  rectangular  plate  is 

d4  d4  d4  ,  phi 

-  +  2  -  +  -  <u  —  W(a:,y)  =  0  (El) 

dx 4  dx~  dy 2  dy4  ^  J 

Using  the  Rayleigh-Ritz  technique,  the  approximate  solution  W{x,y)  of  Equation  (El)  is  ex¬ 
pressed  as  a  doubly  infinite  series  of  products  of  normalized  uniform  beam  modes. 


*The  nomographs  yield  results  for  the  special  case  cited  above  which  includes  the  clamped  plate. 
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*(*.?)-  2  2  ^(y)  (E2) 

m  =  1  n  =  1 

where  the  mode  shapes  and  0n(y)  are  associated  with  the  mode  shapes  of  uniform 

beams  having  end  fixities  which  are  the  same  as  the  corresponding  edges  of  the  plate;  the 
particular  form  of  these  beam  modes  for  particular  boundary  conditions  can  be  obtained  from 
Reference  21.  These  forms  are  not  required  for  the  present  analysis. 

The  kinetic  energy  T  of  ihe  clamped  plate  is* 

T=  —  phfafbW2(x,y)dzdy  (E3) 

2  0  0 


Substituting  Equation  (E2)  into  Equation  (E3),  we  obtain 


T  =  —  a >2  2  a  a  M  <f>  <f>  0C 

2  mnrs  mn  TS  P  m  r  n  5 


(E4) 


From  the  condition  of  orthogonality  of  beam  modes 


0n  8S  =  0  if  n  4  s 

<t>m4r  =°  if®^f 

writing 

T  =-  <o2  M  T 
2  P 

we  have 


(E5) 


(E6) 


T  =  2  a  a  d>  d>  6  6 

‘  ,  mn  rs  ^ ^ r  n  $ 
m ,  n  j  t  j  s  y 


(E7) 


The  integral  expression  for  the  potential  energy  V  of  a  flat  rectangular  clamped  plate 


is 


‘Assuming  no  edge  masses,  all  M-  =  0  in  Reference  21.  With  no  mass  moments  of  inertia  at  the  boundaries, 
all  J  ■  =  0  in  Reference  21. 

“For  the  clamped  plate,  we  assume  infinite  stiffness  in  the  translational  and  rotational  springs  along  the 
edges  of  the  plate  so  that  no  potential  energy  is  associated  with  these  springs.  The  spring  energies  are,  how¬ 
ever,  included  in  the  potential  energy  term  in  Reference  21. 
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a  b 

*  - 1  (  J  w,a,  *  »&  "V, ♦  2d  -  ->  ** 

0  0 

or  (E8) 

a  b  a  b 

r»|  J  j  »U +  J  {  Ky 

0  0  0  0 


Substituting  Equation  (E2)  in  Equation  (E8),  we  get 


Da  f/5  \ 4 - -  - 

2  $3  mnrs  mn  rsL\o/  9m  9r  n  s  9m  9r  n  5 

+  (I)2  IO,  KK* ] 


<E9) 


0 


+  —(!_,/)  s  (2  a  [6' 6'  9'  6'  -  6" 6  6  d "] 

mn  rs  ™m  “r  n  s  “ m  .  r  n  s  -1 


This  equation  can  be  simplified  by  use  of  the  integral  relationships  between  <£m  <f>r, 
<f>"  </>'',  <j>'m  and  6n  0S,  0"  6",  0'  0'  given  by  Equation  (42)  of  Reference  21.  The  steps 
involve  a  lengthy  integration  by  parts.  The  resulting  expression  for  the  potential  energy 
becomes. 

D  a  - 

V  =  -  -  V  (E10) 

2  £3 


where 


2  a  a 
mn  rS 

ui  n  r  $ 


+  £ 
mnrs 


+  2 
mnrs 


<£<£  0  0  +  a  4  0  0  6  6 

m  \a  /  ns  n  n  s  ^  ^  r 

«„  (7)  W 

<■„  12(1  -  (j)  «*;  «“  <»„ «;( 


-(<£'<£)  0  0"  -  (0  0')  6"  6  ] 

'“m  ■  rt.  n  s  '  n  j'.  r'm  vr 


(Ell) 
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AH  the  expressions  necessary  to  evaluate  the  derivatives  and  integrals  of  mode  shape 
appearing  in  Equations  (E7),  (Ell),  (El6),  and  (E17)  have  been  developed  in  Reference  21 
and  are  also  used  in  Appendix  F.  Hence  the  quantities  Cf,sn  and  8™D  can  Le  numerically 
evaluated  for  a  clamped  plate.  Solution  of  the  set  of  Equations  (E15)  can  be  obtained  by 
iteration  using  standard  digital  techniques.  These  methods  are  briefly  discussed  in  Refer¬ 
ences  16,  21,  and  22  for  certain  special  cases. 

In  Reference  21  numerical  evaluation  of  Equation  (El 5)  showed  that  accurate  estimates 
of  the  plate  frequencies  can  be  obtained  by  using  a  single  term  from  the  appropriate  equation 
out  of  the  set  of  Equations  (Elo).  To  obtain  the  approximate  frequency  equation,  set  rs  =  mn 
in  Equation  (E15)  and  equate  to  zero  the  coefficient  of  amn  giving 

1  -L 

where 

kmn  =  “m  +  «.  +  2(J/fl>2  * n  > 

(E19) 

^m  =  KK  /  *m 

J,  =  9"  9  /  9 2 

r  ti  n  n  In  j 

Actually  Equation  (E19)  and  the  quantity  tpm  (or  tpn)  was  numerically  evaluated  for  the 
beam  having  translationally  fixed  ends  and  rotational  spring  ends.  Thus  Equation  (E19)  is 
the  approximate  solution  to  an  equation  somewhat  more  comprehensive  than  Equation  (El5), 
given  by  Equations  (66)  in  Reference  21.  For  a  clamped  plate,  the  rotational  spring  has  infi¬ 
nite  stiffness.  The  results  are  presented  in  Figures  11—13  for  the  first  three  beam  modes. 
Thus  approximate  frequencies  can  be  obtained  for  the  first  nine  modes  of  the  plate  for  any 
aspect  ratio  b/a  by  using  the  above  equation  and  the  data  presented  in  Figures  11—13  for 
\pm  (and  ipn)  and  Figures  14-16  for  am  (and  an).  For  symmetric  edge  fixity  in  which  oppo¬ 
site  edges  are  equally  constrained,  the  numerical  results  obtained  agree  within  2  or  3  percent 
with  those  computed  in  Reference  22  using  a  36-term  series.  The  approximation  is  increas¬ 
ingly  more  accurate  the  smaller  the  plate  aspect  ratio  and  has  the  greatest  error  for  the  square 
plate,  particularly  in  the  fourth  and  fifth  modes  when  equally  constrained  on  all  four  edges. 
Approximate  mode  shapes  <f>mn{x,  y)  ~  <f>m(x )  &n(y)’  locati°ns  of  peak  deflections,  locations 
of  node  lines,  etc,  can  be  obtained  from  the  data  presented  in  Figures  19-53  of  Reference  21. 
A  nomograph  constructed  by  White  is  presented  in  the  present  report  to  aid  in  evaluating  the 
approximate  resonance  frequencies  of  the  plate,  Equation  (E19),  corresponding  to  the  first 
nine  modes  for  any  aspect  ratio  b/a.  The  opposite  edges  can  have  equal  or  different  elastic 
constraints.  Note  that  graphical  techniques  can  account. for  only  the  most  significant  term  or 
terms  in  a  mathematical  solution  which  may  involve  a  large  number  of  terms. 
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.  3.2  3.4  3,6  3.8  4.0  4.2  4.4  4.6 

BEAM  FREQUENCY  PARAMETER  a, L 

Figure  11  -  Parameter  r!>x  versus  a  10  andctj^.  First  Mode 


Figure  12  -  Parameter  \p2  versus  a  2(J  and  a^[  .  Second  Mode 
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FREQUENCY  PARAMETER  a 


Figure  17  presents  nomographs  developed  by  Dr.  White  for  nine  modes  of  a  rectangular 
plate.  These  permit  the  graphical  computation  of  resonance  frequencies  of  a  plate  of  arbi¬ 
trary  aspect  ratio  when  the  four  edges  of  the  plate  are  translationally  fixed  but  elastically 
restrained  against  rotation.  The  compliances  of  the  rotational  supports  are  assumed  to  be 
uniform  along  each  edge,  but  the  compliances  may  be  different  for  all  four  edges.  The 
clamped  plate  is  represented  by  rotational  springs  of  infinite  stiffness  along  all  edges.  Each 
nomograph  contains  a  sample  calculation  which  is  indicated  by  arrows  and  which  is  tabulated 
on  the  nomograph.  Note  that  it  is  necessary  to  transfer  numerical  values  from  certain  scales 
to  other  scales;  these  transfers  are  indicated  by  arrows  at  the  bottom  of  each  nomograph.  If 
opposite  edges  of  the  plate  have  different  rotational  elastic  constraints,  the  tp1  and  «*j  scales 
should  be  used  instead  of  the  £  scales.  Values  of  ax  are  obtained  from  Figure  14  for  unsym- 
metric  edge  fixities.  In  the  nomographs  \/Amn  is  replaced  by  amn.  Symbols  used  in  the 
nomographs  correspond  to  those  used  in  Reference  21. 


> 
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Nondimensiona!  Rotational  Stiffness  Parameters, 


Nondimensional  Rotational  Stiffness  Parameters, 


Nondlmenslonal  Rotational  Stiffness  Parameters 


Ratio. 


u  faT1 


e  -  _ 2_ 

'4  El 


*3¥-TT 

*  10 


Exocplt =  >,  =  20 

W® 

b/o  =  .5 
=21.4 


tO  V?  H 

\  Y-~\  »4 


L_ 


Nondlmenslonal  Rotational  Stiffness  Parameters, 


Figure  17e  —  Nomograph  for  Plate  Nondimensional  Frequency  Parameter  oc-,2 
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Figure  17f  —  Nomograph  for  Plate  Nondimensional  Frequency  Parameter**^ 


Piute  Nonrilmenslonal  Frequency  Parameter 


Nondimenslonal  Rotational  Stiffness  Parameters, 


Nondimensional  Rotational  Stiffness  Paiameteis, 


Nondimensionat  Rotational  Stiffness  Parameters, 


Nondimensional  Fiequency  Parameter, 


Beam  Nondlmenslonal  Frequency  Parameter, 


APPENDIX  F 
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Density  of  material 
Poisson’s  ratio 


a 

4>  Normalized  mode  shape 

ip  Resonant  freque~<iy  parameter 

c o  Circular  frequency 

Subscripts 

m ,  n  Refer  to  m  th  and  n  th  modes,  respectively 

n  Refers  to  direction  normal  to  certain  direction 

r  Refers  to  rth  mode 

x  Refers  to  z-direction 

y  Refers  to  y-direction 
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DESCRIPTION 


Crocker23  presents  an  analysis  for  computing  the  normal  modes  and  frequencies  of  a 
uniform  flat  panel  with  fully  fixed  edge  conditions.  The  method  involves  an  approximate 
solution  of  the  frcqi  vncy  equations. 

DERIVATION 

The  mode  shapes  of  a  clamped-clamped  panel  are  approximately 


xmV)  xjy)  1  r 

fjx,  y)  =  - - -  - -  , 

"  [xmmixam  \xj\xn\  L 


cosh  a„  —  +  sinh  or  _  — 


m  m 


+  C  cos  <*  —  +  Z?  sin  a„ 

i7i  m  ^  m  ci 


•  [*.  cost  «.  l.B,  sinh  l 
cos«„  y_tD'  sin  on  X] 


where  the  quantities  in  brackets  or  Xa,  Xn  represent  the  mode  shapes  of  vibrating  uniform 
beams  lying  along  the  x-  and  y-axes,  respectively,  and  \Xm\  and  !A'n|  are  their  respective 

values.  Applying  the  boundary  conditions  for  a  clamped-clamped  plate,  i.e.,  for  either 

v  y  y  ^X  \x= 0,  y=0 

Xmo,X„X.—  .0^  y=J 

Then 

-C 
B=  -D 

and  (F2) 

0  =  A  cosh  a  +  B  sinh  a  +  C  cos  a  +  D  sin  a 

0  =  A  sinh  at  +  B  cosh  a  -  C  sin  a  +  D  cos  a 

Equations  (F2)  may  be  solved  in  order  to  obtain  the  frequency  equations  for  a  clamped- 
clamped  plate: 


cosh  ai7,  cos  am  -  1 
cosh  an  cos  an  =  1 
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Solution  of  Frequency  Equations 

The  solution  of  Equations  (r'3)  may  be  shown  to  be  of  the  form: 

am  ~  (2m  1)  *  ~  +  A:  m  =  1,2,3...  ,« 

where  A-»0  as  Now 

cosh  am  =  [cosh(2m*l).|]coshA*  [sinh  (2n  +  1)  .  ij  sinh  A 

=  [sinh  (2m  -s-  1)  -  — J  [cosh  A  +  sinh  A] 
and  from  Equation  (F4) 

cos  <*ra  =  -  [sin  (2m  -s-  1)  -  -  J  sin  A,  [jsince  cos  (2m  ~  1)  .  1  =  oj 
=  -  (-  l)m  sin  A 

Thus  from  Equations  (F3),  (Fo),  and  (F6): 

(cosh  A  +  uinh  A)  sin  A  =  - - - ^ _ _ 


(F4) 


(F5) 


(F6) 


sinh  (2m  +  1)  — 
2 


(FT) 


But  A  =s  0.  Thus  for  small  values  of  A, 


cosh  A  =  — 
2 


sinh  A  =  — 


[i  +  At  5!l 

2  2  L  2  •  2  J 

sin  A  =  A 

Thus  substituting  Equations  (F8),  (F9),  and  (FID)  into  Equation  (F7) 


(F8) 
A  (F9) 
(F10) 


gives: 


-2(-D-e"9<"  +  ,>T 


sinh  (2m  +  1)  -  r- 


(Fll) 
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and  neglecting  terms  of  order  greater  than  A,  then: 


A  =  2(-l) 


=  rl 


-  e 


~(2a  -r  1)  - 


(F12) 


lis,»g  Equations  (F4)  and  (Fl-2),  values  of  o ,  to  «„  aero  calculated  in  Reference  23  and 
are  presented  ,n  Table  7.  It  was  found  that  for  the  higher  frequency  parameters,  the  value  of 
4  became  negltgtble  and  Equation  (Frj  was  sufficiently  accurate.  For  example.  A.  ,  -  i  43t 
x  10  and  was  thus  negligible.  Equations  (F3)  may  also  be  solved  by  assuming  a  solutio, 

such  as  Equation  <F4>  with  4 . 0  and  using  the  Newton  method  to  refine  the  original  apprexi- 
mate  solution.  ' 

Determination  of  the  Modal  Constants 

Arbitrarily  putting  one  of  the  modal  constants  l >=  .  1,  the  other  mode!  constants  may 
be  determined  from  Equations  (F2). 

Thus  B  =  -  D  =  -  1  and  Aa  sinh  orc  -  cosh  «c  -  Aa  sin  -  cos  a  =  0. 


cosh  cr  -  cos  a 

C3  £3 

a  sinh  a  ~  sin  a 


(FI3) 


But  using  Equations  (F5)  and  (F9), 


Thus: 


cosh  a 

SI 


(2m  ~  1) 
e 

2 


2 

=  sinh  a 

SI 


COS  Of 

1  ^ 

sinh  -cos  am  ~  sinh  a n 

sinh  a  ~  sin  a  sin  a 

1  ^ 

sinh  a 

SI 


(  COSam\ 


/ 


sm  a_  \ 


TABLE  7 


Parameters  for  a  Clamped- Clamped  Mode  Shape 


mom 

Frequency  Parameter 
a=oraa 

Resonant 

Frequency  Parameter 

3***  3s 

Haaaua 

Disploceaent 

A.  or  A’ 

jb  c 

Model  Coefficient 
ABorAn 

4.73004 

11302 

1.6162B 

1.017804 

7.85320 

44050 

1.50605 

0.999224 

W.99560 

98.905 

L51259 

1.000034 

1413720 

171.590 

1.5122B 

C.999998552 

17.27880 

2641376 

1.5125 

1.0000000627 

20.420352 

3741092 

1.5125 

0.99999999729  j 

21561945 

5048633 

L5125 

1.0000000001175 

8 

26.703537 

659.4048 

1.5125 

0.99999999999491 

mm 

29.845130 

830.7431 

1.5125 

1.000000000000220 

H 

31986722 

- 

1.5125 

0.99999999999999046 

So!c  The  cod<]  coef£dc3ts  C  =  —  A.  aai  B  —  — 

n  as  an 

given  ahere  t hey  are  repaired  for  accurate  calculations. 

=  —  3-  Store  significant  Semes  are 
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Bat  from  Equation  (F6) 


cos  a  _  =  -  (-  1)~  sin  A 


and  from  Equation  (F4) 

=  £sin  (2*  *!)•  —  J  cos  A;  jsince  cos  (2ej  -r  1)  -  —  =  oj  =  (-  1)“  cos  A 


sin  a. 


Thus 


--1=1- 

Si 


I (_  i)®1  .  cos  A  -  (-  1)“  •  sin  A] 


sin . 


(-  1)~  •  {cos  A  -  sin  A] 

Aa  =  1  -  - — i - 

“  sinh  a_ 


1)°  *  !-l  -  A]  -  2e  (2“*I>  2  (Fl4) 


Thus  using  Equation  (Fl4)T  values  of  Aa  and  CB  were  calculated  for  as  =  1  through  10;  see 
Table  7 

Determination  of  Resonant  Frequency  Parameters 

Eh3 

From  Equations  (E19)  of  Appendix  E  with  Aac-»/?_c  and  D=El==  -  ,  the 

12(1  —o2) 

undamped  resonant  circular  frequency  of  the  cm  th  mode  of  the  plate  is: 

t’=n='/^7  -J — 

b2  F  12p(l  - o2) 

where 

fb\*  .  .  fb\2 

R  =  I  —  I  .  o4ia4  +  *2  —  I  -ib  •>, 

mn  yaJ  n  n  \a  )  ■  n 

a  a  was  derived  above  in  this  appendix  and  values  are  given  in  Table  7.  Also  the  following 
relations  were  derived  in  Reference  21  and  used  in  Appendix  E. 


(F15) 


(F16) 


where 


-  Try  [(J“*B=)sinh  *  7  *  (e  ““  *-c°s«a*y)*sin«a-yJ 

=  P7j-  taa-  B=>  00511  «»  *  J~  *  Aa  (“«  f  *  Sln  **«=  *  f)  ^  COS'“»  '  J  ] 

=  jyy  [(4:*  S=)  sinb  “a  *  7  *  («  *“  *  *  cos  «e  *  y)  -s«°  «„  *  7J 

^="  =  jjy  coshaa  -  7  *4,  (-e  “  F*sin  «B  -  y)  -cos«=  -  lj 


(FIS) 


and  where 


=  M  '-ib-4 


*I*J2 


ib^-a^c^dB 

C3  J S  C5  CS 


(FI9) 


—  3 

o  ~  =  — 


zero  zero 

-\,A  ,-L-kd  *_L 

J0  4-^1/  J0  2|Ac 


U2-B2-Ci  +  £?2] 


£3  £3  S3  £3 


(F-20) 


The  terms  shown  zero  in  Equations  (F19)  and  (F20)  are  zero  due  to  the  boundary  conditions 
d>m  -  =  0  at  x  =  0  and  x  =  t  for  a  clamped-ciamped  mode. 

Substituting  Equations  (F19)  and  (F20)  into  Equation  (FIT)  and  utilizing  Equations 
(F18)  gives,  after  simplification. 


^m=  {-7  *  {(Aa  -  BJ  cosh  a  m  +  A(-e  m  -  sin  aj  -  cos  aj  [  (Am  +  Bm)  sinh  a, 
+  A  (e  m  +  cos  a  )  -  sin  a  ]  *  2.4  (1-Z?  )  1  -a2  [B  2-j4  2  +  d?2  +  Z?2]l  / 

cx  ci'  in  si  '  m*  in  Z7Z  m  ci  ml/ 


(4  2  -  B  2  +  C  2  +  0  2] 

m  m  nj  in 


(F21) 


Equation  (F21)  was  evaluated  for  m  =  1  through  9;  the  values  are  presented  in  Table  7. 
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Value  of  Position  of  Maximum  Displacement  for  Each  Mode 

In  order  to  simplify  the  computer  program  for  the  response  of  a  clamped-clamped  panel, 
it  was  necessary  to  determine  the  maximum  value  Xa  ,  denoted  |A'n|,  for  each  mode.  In  fact, 
the  simplest  and  most  accurate  method  found  was  to  calculate  the  mode  shape: 

[x  x  x  i  "| 

Aa  cosh  <*_  *  -  *  Ba  sinh  «n  *  -  +  CB  cos  a=  -  -  -  Da  sin  an  -  -  \ 

(F22) 

by  means  of  a  computer.  The  computer  program  written  by  Crocker  is  given  in  Figure  18. 

Both  numerical  values  and  computer  plots  were  obtained  for  m  =  1  through  10,  and  the  computer 
plots  are  given  in  Figures  19—23.  In  this  manner,  both  values  of  [Xe|  and  A’  for  x  =  a/2 
were  obtained.  Since  the  whole  mode  shape  was  calculated,  the  response  of  any  point  of  the 
panel  could  be  computed  by  using  the  appropriate  values  of  A‘=  (x)  and  A’n(y).  It  is  interesting 
to  notice  that  Figures  19—23  indicate  that  the  maximum  displacement  |A’c|  does  not  occur  at 
the  center  of  the  span  except  for  the  first  mode,  but  two  maxima  lA'^j  occur  for  the  higher 
modes,  one  nearest  to  each  support.  The  other  maxima  are  found  to  be  slightly  smaller,  to  be 
of  approximately  constant  value  for  the  higher  modes,  and  to  lie  between  positions  of  the 
maxima  |A’  |. 

An  approximate  method  is  given  below  for  determining  the  value  and  position  of  the 
maximum  displacement  |A'=1  for  the  higher  modes.  Although  approximate,  |A’mj  calculated  by 
this  method  is  seen  to  be  only  0.66  percent  smaller  than  w-hen  calculated  by  the  more  exact 
computer  program. 

The  mode  shape  as  given  by  Equation  (F22)  may  be  rewritten: 


A'_  =  (A„  ~  B  )  sinh  a  ■  —  +  A„ 
=  ~  a1  m  a  02 


(  -anx  a  \  aX 

I  e  -  cos  -  1-5-  sin - 

\  a  J  a 


(F23) 


but  for  a  maximum  or  minimum  value  of  X  : 

El 


f dX*  \ 

-rhr )-° - 


(4o  +  BJ  COsh 


(-<*  X 

or  x  \  ax 

El  1  m 

-e  ~  sin  -  I  +  cos  - 

a  /  a 


(F24) 


Since  for  the  higher  modes-. 


(F25) 
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- aip«h  - 
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- Hi  ;■*!,*  Jt-3S - 
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- moiRTWESW - 
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M«l— 5.95CE-I2 
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HtiOI*-9.5*E-iS 

Wl  1  ~  . 

10  t*3. 

- AC  TO  1S031S3I  nS^TC^t  1  i - 

C  «T  *Cw  ORICIn 

— 5CC — UU  PLCTtU.-|L77"3I - 

c  0*AU  X-AXIS 

- urn.  pl'ju  wauCTTi - 

c  tRaW  TIC  HAVAS  9*  X-AXIS 

- mm - 

»2*.I2S 

- TTTTZ - 

00  SC 2  1*1*10 

- LAM.  FLaUXTIIfOi - 

CALL  FL0TCX.X2»2I 

-c - im  cmo  os  i»»m - 

CALL  STK*L*1X«--2S».I4»4CC< I l»0.*2J 

—502 — X*T*T  V - 

C  *RITE  *•  AA3  ITS  VALUE 

- LALL  >fH5LXtl,5SI.T5l.l41JWWI*.li  J  ■ 

CALL  K7i*)£4(7.tS>l.7St.l«Hij.,}Ml)) 

t - 3tW  T-ATI3 - 

CALL  FL0T(0.«2.»3> 

- CALL  PLOT?  5. 1*2. *2  7 - 


0RAW  TIC  HAVAS  On  T-AXIS 


XI  *0* 

X2-.I2S 

-T*=7Z - 


“557 


00  S03  I*l » 9 
Call  Piouiifffjj 

CALL  FLQTCX2»T»2I 
t*t**5 
*S»|TC  LASELS  cv  T_*X|S 
’  CALL  SlSsSt— J - 


*5»3.v.K*Ih  J.fl.B,.*] 

CAit  Sri5.«C-.«5.  . . AH  USm.iO 

CaCl  5»«BL«| -.*-»»  I. g*.T*»*H  |.a»0«»A) 
CALL  SfH8U(-.H»  O.S».l«»«w  a.S»u*»«> 
call  5iHsX«f-.*9»o.». > 
CALL  st**9L4(  —  *4S«-0*S** I a«ah-o*S*o« *4  1 

mi  STnL4l-*l9>-I.B.*U.lu-l.j»J..I » 

CALL  STHf L4C  7 

CALL  5IOSUl-.»J»-2.n..,U»«4-J.U.Uf«  ) 

WRITE  0  ANO  X  Ov  CR10 

CICC  SI WLAI‘-*LI,J. I* IX l|MLT5*J*<  1 1  J - 


HflHXiO.il  I 


CALL  STH5L«CS.O*-.7S* 
tftIUUR  IO~OUIUR  ■"  1  1  1 

CALL  PLOT(o.fO**31 
_  lOHT (nut 

C*B(  n  )•(  EXPft  ALPHA!  n  )•  t  i-E  tPi-'C  -alpha!  n 

'U»AI  W  TWCTPn‘-ACFWA|  r.  I»X  J - 

E-AC m )*COSF( ALPHA! H J*X  ) 

TXJTXfTTCWUl  K  1111  1 

G*C*0-E*F 

CO~TO  C50XISUS  IISSRILRri I  J - 

CONTINUE 

-plot  PJINI5 - 

XP*IO.*X 


- LALL  PLOrtAPILUl  ‘  ”  - 

CO  TO  S06 

503 — CONTINUE - 

WRIT£CCI«2)H»X'CtDft*F*6 

- 2 '  '  rURWXU  miJlbCLIA.3  U  '  '  '  - 

SOS  CONTINUE 

- *-A*.W^i3  —  ■  ■  ■ 

iFC X-| .  ) I • I » 3 

- 3 - CONTINUE - - - - 

I  Ft  N-IO)<fS*S 

- « - It»trri - 

TO  TO  10 

—5 - STOP  — —  ’  -  ~  - — . . . 

END 

Figure  18  -  Program  to  Calculate  and  Plot  Clamped-Clamped  Mode  Shapes 


i 


This  program  is  not  the  one  used  at  NSRDC  to  obtain  the  'Vequencies. 
The  NSRDC  program  is  given  in  Appendix  I. 


MOPE  SHAPE  A.  MODS 


I 


it  is  seen  by  inspection  of  Equation  (F24)  that  the  first  maximum  will  occur  at: 


x  3 

a_  •  —  =  —  77+8 

m  a  4 


where  8  is  a  small  number. 

Thus  making  the  approximations  then  cosh  5=1  and  cos  5  =  i, 


x  Sir  Sit 

cosh  a  •  —  =  cosh  —  +5  sinh  — 


X_  -3  TT 

e~*m  °  =  (1  -8)e  4 


x  -1 

sin  a  „  •  —  =  — ^=r  (1-5) 


(F26) 


(F27) 


cos”'"'  f '  7f 


Then  substituting  Equations  (F27)  into  Equation  (F24): 


—3  TT 


(^  +  SJCosh  Y+  5(4m  +  8m)  Sinh  —  ~(l-S)Ame  4  +-^~(1-S)-  -jf=-  (1  +  5)  =  C 


and  thus 


3  TT 


-3  TT 


5  = 


-( Am  *  Bm >  c0sh  —+Ame  4  +  (*  “  AmV^ 


(F28) 


3  TT 


~3  TT 


(• Am  *  Bm)  sinh  —  +  Am  e  4  ~  (1  r  AmVJ2 


Using  the  approximations  in  Equations  (F25),  Equation  (F28)  reduces  to: 


0-3  rr/4 


5  = 


0.0948 


e— 3  7t/ 4  _  J2  0-0948  -  1.4142 


„  0.0948 

5=  -  =  0.0719 

1.3194 
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N 


Again  using  the  approximations  of  Equations  (F25)  and  (F27),  Equation  (F23)  reduces  to: 


(F29) 


(F30) 


The  value  of  \Xm\  obtained  by  the  above  approximate  method  and  presented  in  Equation  (F29) 
compares  well  with  the  computed  values  (presented  in  Figures  19-23)  and,  in  fact,  is  only 
about  0.66  percent  smaller.  The  position  of  the  maximum  displacement  as  given  by  Equation 
(F30)  is  also  in  good  agreement. 
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APPENDIX  G 


THE  SUN  METHOD 


NOTATION 

[-4] 

.4. 

I 

Amn 

a 

IB] 
b 

IC ] 

D 

F 

G\Gl 

9 

h 

L,  L ' 
m,n 


P 


P,Pi 


R 

T 

t 


Symmetric  square  matrix  or  order  n  whose  elements  are 
defined  by  Equation  (Gl4b) 

Coefficient  in  equation  for  displacement  surface  function 
Coefficient  in  equation  for  displacement  surface  function 
Length  of  rectangular  plate 
Symmetric  real  matrix  defined  by  Equation  (Bio) 

Width  of  rectangular  plate 

Symmetric  square  matrices  of  order  n  whose  elements  are 
defined  by  Equation  (Gl4a) 

Eh3 

Flexural  rigidity  of  plate  equal  to  - 

12(1  -  a2) 

Function  satisfying  the  boundary'  condition  for  clamped  plate 
Polynomial  in  equation  for  displacement  surface  function 
Acceleration  due  to  gravity 
Plate  thickness 

Defined  by  Equations  (Gloa)  and  (Gl8),  respectively 
Mode  numbers 


Equal  to  R~@  = 


b\~P 


Circular  natural  frequency;  pi  = 
i  =  1,  2 . n 


b 


Equal  to  — 
a 


where 


Kinetic  energy 


Time 
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Preceding  page  blank 


Potential  energy 

Surface  displacement  function  of  plate  in  direction  perpen¬ 
dicular  to  plate;  subscript  t  indicates  a  time  derivative 

x  y 

Equal  to  —  and  —  ,  respectively 
a  a 

Column  matrix  containing  elements  of  A'  where  X  =  L'v 

Variables  in  cartesian  coordinate  system 

Exponent 

Exponent 

Plate  mass  per  unit  of  surface  area  where  y  is  the  weight 
per  unit  volume  of  plate 

a2  a2 

Equal  to  -  +  - 

dx~  dy2 

Kronecker  delta 

Equal  to  — 
a2 

Poisson’s  ratio 

Transverse  displacement  of  plate  in  free  vibration; 
subscript  t  signified  a  time  derivative 

Column  matrix  of  A.,  A2 . Ai . An  defining 

the  eigenvector  of  the  specific  natural  mode  concerned,  i.e., 
nodal  pattern  of  i  th  vibration  mode 

Eigenvalue  defined  by  a  =  p 
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DESCRIPTION 


Son24  presents  a  method  for  computing  the  normal  modes  end  frequencies  for  a  clamped 
thin  rectangular  plate  undergoing  transverse  vibrations.  Vertical  shear  and  rotary  inertia 
effects  are  ignored.  The  method  uses  the  Rayleigb-Ritz  procedure,  but  the  deflection  of  the 
plate  is  represented  by  a  series  of  polynomials  rather  than  the  product  of  beam  norma]  mode 
functions. 

DERIVATION 

The  transverse  displacement  for  a  freely  vibrating  thin  plate  is 

=  »"(r,y)  cos  pi  (Gl) 


The  potential  energy  of  the  plate  is 


V  =  JJjdF  =  JJW*  *  *yr  *  ~G  -(1  “  c)  dxd9  (G2) 


The  kinetic  energy  of  the  plate  is 


,  d.  ff«2 

2?  Jj  ' 


dx  dy 


Substituting  Equation  (Gl)  into  (G2)  and  (G3)  and  setting  cosine  and  sine  values  equal 
to  1  in  Equations  (G2)  and  (G3),  respectively,  the  maximum  potential  and  kinetic  energies  are 

,rnai  -  |  jjW  R  )2  -  2(1  -  a)  [Hrx  Ryr  -  B  jp  I  dxdy  (G4) 

raax=  ^  p2j\w2dzdy  (Go) 

Equating  Equations  (G4)  and  (G5)  as  required  by  the  Rayleigh  principle 


Zg  n  ax 

p  =  -T  - : — 


vh  //II'2  dx  dy 

Now  there  is  a  class  of  plate  geometries  governed  by  the  equation 


7  "  f  - 


S7 


1 


Equation  (G7)  includes  the  approximated  rectangle.  Dividing  through  Equation  (67)  by  c  and 
letting  X  =  z/c,  Y  =  y /a,  R  =  b/a,  P  -  B~@,  the  resultant  normalized  equation  replacing 
Equation  (67)  is 

Xa  +  PYP=  1  (68) 


Then  to  determine  the  natural  frequency  p  of  the  clamped  rectangular  plate  in  terms  of 
a,  /5,  and  ?,  let  the  displacement  surface  function  be  expressed  as 

W  (X,  T,P,a,fi)  =  F  (X,  Y,  P,  a  ,0)  1  1  X*  7* 

a  =0  as  =0 

=  F(X,T,Pta,P)(Aoo+AJ0X*A01T*AllXT*..  (69) 


=  F  2 

i=  1 


where  for  a  clamped  plate 


F  =  (1  -  Xa  -  PY&f 


(GlO) 


aw  aw 

satisfies  the  requirement  — —  =  — —  =  ,  B  =0  along  the  boundaries. 

dX  dY 

Following  the  Rayleigh-Kitz  procedure,  the  A- s  in  Equation  (GS)  have  values  obtained 
from  a  minimization  of  Equation  (64). 

—  [jj{(v* *02  -  2  (i  -  o)  iws=wn  -  b;2}J  -  fr-2}  dAV/y  (Gn) 


z=  1, 2, ...  n 

For  the  clamped  plate,  satisfaction  of  the  natural  boundary  conditions23*  (also  see  Appendix 
B)  reduces  Equation  (611)  to  the  simpler  form 


l  [III 


(y2  fjf2  J  =  o 

r  =  1, 2, ...  n 


(G12) 


*There  are  no  natural  boundary  conditions  for  the  clamped  plate  and  therefore  they  need  not  be  satisfied. 
However,  as  discussed  in  Appendix  B,  practical  consideration  of  the  rate  of  convergence  makes  such  satisfaction 
desirable. 


/ 


Substituting  Equation  (G9)  with  F  =  (1  -  A'**  -  PY&)  into  tfee  above  equation,  a 
matrix  equation  results  as 


W)  lr>|  -  o2  f.4]  \6\  =  0 


(Gl3) 


where  1-4]  and  ]C]  are  square  matrices  of  order  n  whose  elements  are  respectively  defined  as 


l  /fa-Y*>£ 

C{LJ)  =  /  /  v2  (FGJ)  V'2  (FGJ)  dXdY 

o  o 


(Gl4a) 


,  l  fi{ \-Xa$  ,  . 

-■*  (/,  J)  =  I  f  ( FG ')  ( FGJ )  dXdY 


(Gl4b) 


where  F  =  (l  -  A'a  -  PY&)2. 

Matrices  fCl  and  1.4]  are  therefore  symmetric  square  matrices  with  all  real  number 
elements. 

The  column  matrix  \M  of  .4j,  .4, . .4t- . .4c  defines  the  eigenvector  of  the 

specific  natural  mode  concerned  and,  in  turn,  yields  the  modal  patterns  of  the  corresponding 
vibration  mode. 

The  eigenvalues  of  Equation  (GI3)  are  o2  =  p~  ( yh/gD )  where  p  is  the  natural 
frequency. 

In  order  to  reduce  Equation  (Gl3)  to  standard  matrix  pencil,26  let  C  =  LL'.  A  =  l/&>2. 
and  A'  =  L" &.  Equation  (G33)  then  becomes 


L -1  A{L  T1  A  =AA' 
or  IB]  1A|  =  A|A| 


(Gloa) 

(Glob) 


where  [B]  is  symmetric  and  real  and  thus  }A }  is  orthogonal  with  respect  to  each  natural  mode, 
that  is27 


a;  A'.  =  5-. 


(G16) 


where  8~  is  Kronecker  delta.  Tne  natural  frequencies  can  then  be  expressed  as 


(G17) 


and  the  corresponding  eigenvectors  !^{-l  can  then  be  obtained  through  the  following  trans¬ 
formation: 


=  T1  IX  il 

The  modal  pattern  of  the  i th  vibration  mode  is  given  by  ]$.}. 


(G18) 
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To  achieve  a  good  approximation  to  the  fundamental  and  higher  mode  frequencies.  Sun 
used  an  xy  (or  XY)  polynomial  consisting  of  21  terms.  The  computational  methods  include 
both  a  beta  function  evaluation  and  a  Gaussian  quadrature  integration  technique.*  The  latter 
has  no  restriction  as  to  the  values  of  a  and  J3  but  requires  approximately  twice  the  computa¬ 
tional  time  of  the  former.  The  method  of  reduction  (i.e.,  iteration)  is  used  to  find  the 
eigenvalues  and  the  corresponding  eigenvectors  are  obtained  from  Equation  (Gl5b). 
Polynomial  expressions  for  the  fundamental  and  higher  modes  as  well  as  other  details 
relevant  to  the  computational  methods  are  given  in  Reference  24.  The  reference  also 
includes  computed  results  which  were  carried  out  on  an  IBM  7094. 


•When  Of  and  (3  values  are  less  than  or  equal  to  1.5,  the  beta  function  is  not  properly  defined.  Hence,  a 
numerical  integration  using  the  Gaussian  quadrature  r  lie  of  order  64  was  used  in  the  range  below  a  =  /3  =  1.6. 
A  Gaussian  quadrature  double  integration  formula  is  given  in  Appendix  B  of  Reference  24. 
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APPENDIX  H 


THE  CLAASSEN-THORNE  METHOD 


NOTATION 

a 


4> 

d  , 

h  , 

ft* 

C  , 

0  . 

171 5 

-  n' 

€  , 

i 

n 

E 

1 

h 

K 

K' 


tC 


k' 
m,  n 

t 

H '(X,Y) 

X,Y 


x,y 


Plate  length  lying  along  2>axis 

Coefficient  of  doubly-infinite  Fourier  series  defined  by 
Equation  (H6) 

Plate  width  lying  along  y-axis 

Coefficients  of  Fourier  series  defined  by  Equation  (Ho) 


Young’s  modulus 
Frequency 
Half- thickness 
a2 

Equal  to  —  K , 

r,2 


Equal  to 


K_ 

k2 


Equal  to 


i 


Sp(l-u2)  (2 nfj1 


Eh 2 


Equal  to  — 
b 


Equal  to  — 


Harmonic  order  for  sine  waves  along  x  and  y,  respectively; 
see  Equation  (H5) 

Time 

Amplitude 

Rectangular  coordinates 
jt  n 

Equal  to  —  X  and  —  Y,  respectively 
a  b 


Poisson’s  ratio 


i/,  a 

P 


Mass  density  of  plate 
Phase  angle 


'  >^<5  fy  ,  '  J*^10***  {*t**f;  **  i 


DESCRIPTION 

Claassen-Tfcorne10  present  a  Fourier  series  method  for  computing  the  frequencies 
and  modes  of  free  transverse  vibrations  of  thin,  rectangular,  isotropic,  fully  clamped  plates.* 
Corves  axe  given  for  determining  the  first  ten  frequencies  and  their  modal  patterns  as  a 
function  of  the  aspect  ratio. 

DERIVATION 

The  governing  differential  equation  for  sinusoidal  free  vibrations  of  a  thin  rectangular 
isotropic  plate  is 

d*us  d4tc  94w  3p(  1-v2)  d2w 

- *  o  „ -  - - — - -  -  (HI) 

a.v4  a  A'2  ay2  ay4  Eh2  at2 

For  sinusoidal  vibrations,  tc( X,  ?,  t)  =  H’ (A,  Y )  sin  (2 rtf  l  +  6)  Equation  (HI)  becomes 


v  v  if  y  n  _ 

-  *2  -  +  -  =  K2  W 

dX4  3X2dY2  dX23Y2 


where  K2  = 


3p(l  -  v2)  (2  =f)2 


For  a  clamped  plate  the  boundary  conditions  are 


W  (A,y)  =  0 


where  the  subscript  n  denotes  the  normal  derivative. 

The  origin  of  the  rectangular  coordinate  system  is  taken  at  one  comer  of  the  plate, 
with  one  side  of  length  a  iying  along  the  A^-axis  and  the  other  of  width  b  along  the  T-axis. 
Thus.  Equation  (HI)  is  valid  for  0  <  X  <  a  and  0  <  Y  <  b. 

77  it  a 

It  is  useful  to  transform  the  coordinate  system.  Let  x  =  —  X,  y  =  —  Y,  k  =  —  ,  and 

abb 

a 2 

K  =  - K j.  Then  Equation  (Hi)  becomes 

7T2 

a4w  ,  d4w  d4w  ,  o  <  *  < 

-  +2k2  -  + -  =K2W,  (H4) 

dx4  dx2dy2  dy4  0  <y  <n 

*The  frequencies  and  modes  are  also  computed  for  plates  with  two  edges  clamped  and  two  edges  free. 


I 

/ 


A  solution  for  W  is  assumed  to  be  in  the  form  of  a  doubly  infinite  Fourier  series 


0  <  x  <  n 

W  ( x,y )  =  2  2  a  sin  nx  sm  my,  (Ho) 

m  n  ‘  0  <  y  <  tt 


where  2  denotes  2  and  2  denotes  2  . 

Ft  m  =  1  n  n  -  1 

Further  Fourier  series  that  are  assumed  to  exist  for0<£<!7or0<y<;r  (i.e.,  the 
boundary  conditions)  are: 


£ 

K 

sin 

my 

lK(0,y)  = 

2 

Cm 

sin 

my 

m 

m 

W(x, 

tt)  = 

2 

fn 

sin 

nx 

W{x,0)  = 

2 

9n 

sin 

nx 

n 

n 

y)  = 

2 

m 

dm 

sin 

my 

2 

m 

sin 

my 

tt)  = 

2 

rt 

K 

sin 

nx 

Wyyix,  0)  = 

2 

n 

*„ 

sin 

nx 

(H6) 


d2W 

where  =  -  ,  etc. 

dx2 

The  authors  apply  an  available  technique  to  Equations  (H5)  and  (H6)  to  obtain  formulas 
for  the  higher  derivatives  and  cross  derivatives  of  the  Fourier  series.  These  results  are  then 
used  to  obtain  a  solution  for  each  amn  of  Equation  (H5)  in  terms  of  the  coefficients  in  Equa¬ 
tion  (H6).  Higher  derivatives  and  cross  derivatives  required  by  Equation  (H4)  are  then  ob¬ 
tained  from  Equation  (H5)  using  the  solution  obtained  for  each  amn.  Moreover,  since  the  de¬ 
flection  on  all  edges  and  corners  is  zero  for  the  case  of  a  clamped-clamped  plate, 

^m~cm~^n  =  9n~^ ®)  =  ^  (n> ®)  =  ^ (05  n)  =  ^ n)  =  0-  Also  the  normal  derivatives  at  all 
four  edges  are  zero  so  that  Wy  (x,  0)  =  Wy  ( x ,  tt)  =  Wx  (0,  y)  =  Wx  (tt,  y)  =  0.  Finally,  applying  to 
Equation  (H4)  these  boundary  conditions  as  well  as  the  higher  derivatives  and  cross  deriva¬ 
tives  previously  obtained,  an  infinite  set  of  homogeneous  equations  is  obtained.  The  authors 
then  present  a  method  for  the  approximate  determination  of  K  satisfying  these  equations. 

For  the  completely  clamped  plate,  K's  are  graphed  only  for  0  <  h  <  1.  Setting  k'  =  — 

k 


and  K'  = 


—  ,  a  value  of  K  can  be  found  for  k  >  1  by  locating  the  value  of  K  for  —  and  multi 
k2  * 


plying  by  k2.  Appendix  I  gives  the  method  for  determining  the  frequency  from  these  quantities 
as  well  as  a  sample  computation. 

The  frequency  and  mode  data  computed  in  Reference  10  are  presented  there  in  both  tab¬ 
ular  and  graphical  form.  Interpretation  of  the  results  are  given  as  well  as  computer  times  in¬ 
volved  ;n  obtaining  the  results.  A  copy  of  this  reference  is  available  in  the  computer  files 
associated  with  this  investigation  at  the  Computation  and  Mathematics  Department. 
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APPENDIX  I 


COMPUTER  PROGRAMS 


Appendixes  A— H  have  presented  several  methods  for  computing  the  natural  frequencies 
of  vibration  of  clamped-clamped  plates.  The  corresponding  computer  programs  including  flow¬ 
charts  are  given  here;  computer  program  decks  are  now  available  at  the  Computation  and 
Mathematics  Department  of  NSRDC.  Table  1  gives  the  results  of  these  programs  for  particular 
plate  input  data  representing  the  plate  geometry  and  mass-elastic  properties.  Figures  2  and  3 
are  plots  of  the  data  in  Table  la  only.  Thus,  the  first  set  of  results  shown  in  Table  la  con¬ 
tains  the  computed  frequencies  for  a  plate  with  geometry  and  properties  identical  to  those 
used  by  Izzo  (Electric  Boat)1;  experimental  results  cited  by  Izzo  are  also  included.  The 
second  and  third  sets  of  results  shown  in  Tables  lb  and  lc,  respectively,  are  the  computed 
and  experimentally*  obtained  frequencies  for  two  plates  used  by  Wilby.11  The  corresponding 
input  data  for  the  three  sets  of  results  are: 


Data 

Plate  1 
(Izzo-Electric 

Boat) 

Plate  2 
(Wilby) 

Plate  3 
(Wilby) 

Dimension  in  ^-direction 

2.0 

ft 

4.0 

in. 

4.0 

in. 

Dimension  in  y-direction 

2.33 

ft 

2.75 

in. 

2.0 

in. 

Plate  thickness  h 

0.0313 

ft 

0.015  in. 

0.015  in. 

Young’s  modulus  E 

4.175  x  109 

lb/ft2 

33.7 

x  106  lh/in.2 

31.0 

x  106  Ib/in.2 

Poisson's  ratio  a 

0.33 

0.3 

0.3 

Weight  density  pw 

466.56 

lb/ft3 

0.27 

Ib/in.3 

0.27 

Ib/in.3 

Gravitational  constant  g 

32.2 

it/  sec2 

386.4 

in./  sec2 

386.4 

in./  sec2 

Five  sets  of  computer  programs  and  one  manual  method  of  computation  are  presented. 
Their  designations  and  the  computers  used  in  making  the  calculation  are: 

1.  WCGFRE  on  the  IBM  7090  of  NSRDC:  This  program  includes  the  methods  of  Warburton 
(Appendix  A),  Crocker  (Appendix  F),  and  Greenspon  (Appendix  D).  Figure  24  presents  a  flow 
chart  of  this  program. 

2.  WHITE  on  the  IBM  7090:  This  program  treats  the  conversion  of  the  White  nomographic 
values  (Appendix  E)  to  dimensional  frequencies.  Figure  25  presents  a  flow  chart  of  this 
program. 


*The  measured  frequencies  were  obtained  by  Wilby  in  Reference  11. 
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3.  PLFREQ  on  the  IBM  360/9 1  of  the  Applied  Physics  Laboratory,  Johns  Hopkins 
University:  This  program  treats  the  Ballentine-Plumblee  method  (Appendix  C).  Figure  26 
presents  a  flow  chart  of  the  program. 

4.  SUNFRE  on  tho  IBM  360/91:  This  program  treats  the  Sun  method  (Appendix  G). 
Figure  27  presents  a  flow  chart  of  this  program. 

5.  YNGFRE  on  the  IBM  360/91:  This  program  treats  the  Young  method  (Appendix  B). 
Figure  29  presents  a  flow  chart  of  this  program. 

6.  Claassen-Thorne  manual  method  of  computation. 

In  all  computations,  the  frequency  /  (in  hertz)  is  obtained  as  the  product  of  the  fre¬ 
quency  parameter  Am  n  (or  am  n)  and  a  factor.  For  particular  computations,  the  factors  are: 


Warburton: 


r48pm(l-o2) 


Crocker: 


Greenspon: 


Plumblee: 


2 nb2  Vl2pm(l-£r2) 

Vl2pm(l-^2) 

^  Pm?3j(*-Cr2)  ( 


= 

=  Pm) 


Young:  -  y 

12 pm  b2a(\~a2) 

h 

.r~e 

>*  nite. 

2  na2 

Y  12pm(l-cr2) 

h 

1  E  ... 

bun:  - 

2  na2 

V  12y(l-*2) 

k2kn 

J  E 

Claassen-Thorne: 

2a2 

V3Pm(l-o2) 

NOTE:  The  user  submits  weight  density  pw  which  is  converted  by  the  program  to  mass 

Pw 

density  pm  where  pm  = - . 


t 


l 

§ 

K 

? 

fc- 

r 

c 

V 


f. 


WCGFHE  (see  Table  8  and  Figure  24) 

This  combined  program  yields  separate  solutions  corresponding  to  the  Warburton, 
Crocker,  and  Greenspon  methods.  The  program  card  IOPT  contains  data  input  to  the  program 
which  permit  the  user  to  compute  the  natural  frequency  for  either  one  or  all  of  these  methods, 
i.e.,  IOPT  =  1  -►  Warburton  method,  IOPT  =  2  Crocker  method,  IOPT  =  3  -»  Greenspon  method, 
IOPT  =  4  -►  all  of  these  methods. 

Warburton13  treats  the  frequency  parameter  subscripts  m,n  as  the  number  of  nodal 
points  along  the  plate  length  and  width,  respectively;  see  Appendix  A.  However,  most  other 
authors  treat  m,n  as  the  mode  numbers  along  these  dimensions  (or  define  it  for  the  opposite 
dimensions).  Thus  (m  =  2,  n  =  3)ffarbtlrtoa  means  the  1,2  mode  containing  2  nodes  along  x  and 
3  along  y  whereas  (m  =  2,  n  =  3)0thers  means  the  2,3  mode  containing  either  3  nodes  along  x 
and  4  along  y  or  4  nodes  along  x  and  3  along  y  depending  on  the  definition  of  m,n  with  re¬ 
spect  to  the  x,  y  coordinates.  To  avoid  confusion  and  for  compatibility  with  most  investi¬ 
gators,  the  program  assigns  the  modal  (not  nodal)  meaning  to  m,n  for  all  computations. 

WCGFRE  Restrictions 

For  IOPT  =  3,  M  -  5,  A'  -  5.  That  is,  the  Greenspon  option  computes  the  frequencies 
for  M  -  5  and  N  -  5.  However,  for  this  option,  if  the  user  requires  higher  modes  he  may  change 
the  Greenspon  subroutine  to  read  in  the  values  of  the  integrals  discussed  in  Appendix  D. 

The  integrals  are  given  in  References  7,  8,  and  9. 

The  simply-supported  frequencies  may  be  computed  by  the  Warburton  method.  In 
this  case,  the  value  of  SPEC  must  be  1.0.  Clamped  frequencies  are  computed  with  any 
value  of  SPEC  not  equal  to  1.0. 

Units 

All  length  units  are  shown  in  feet.  However,  if  all  length  data  are  converted  to  inches, 
this  is  acceptable  to  the  program,  and  is  actually  preferable  in  the  case  of  a  very  small  plate 
because  of  simpler  handling  and  greater  accuracy. 
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TABLE  8 

Program  Listing  for  WCGFRE  Computer  Program 

COMMENT  ****  PROGRAM  IfCGFRF  ********************* 

COMMON  M.N.A.8.H.E.SIGMA.RHO.PI.G 

************************************************************************* 

M  -  MODES  IN  X  DIRECTION 
N  -  MODFS  IN  Y  DIRECTION 
A  -  LENGTH  IN  X  DIRFCTION 
8  -  LENGTH  IN  Y  DIRECTION 
H  -  PLATE  THICKNESS 
E  -  YOUNGS  MODULUS 
SIGMA  -  POISSONS  RATIO 
RHO  -  PLATE  DENSITY 
G  “  ACCELERATION  DUE  TO  GRAVITY 

***************************^********-»#********#^##*sr  *##**#**-HHf**#%-*** 

PI=3. 14159?7 

REAP (5 .2)  IOPT»  NC ASF 

DO  500  L=1»NCASE 

READ(5.2>  M  eN 

READ (5*3)  A.B.H 

READI5.4)  E.SIGMA.RHO  .6 

2  FORMAT (215) 

3  FORMAT (3F12»6 I 

4  FORMAT (E16.8.3F12.6) 

RHO=RHO/G 

GO  TO  ( 10*20*30  .  10)  «  IOPT 
10  CALL  HARB 

IF(IOPT.LE.l)  GO  TO  500 
20  CALL  CROCK 

IFU0PT.LE.2)  GO  TO  500 
30  CALL  GRFEN 
500  CONTINUE 
STOP 
FND 

SIBFTC  WARBER 

SUBROUTINE  WARB 

REAL  LAMBDA. JX.JY.K.KP 

DIMENSION  OMEGA  120*10) 

DIMENSION  FREQ( 25. 10) »  GX( 100 ) .HX < 100 ) »JX t 100 ) .GY( 100) .HY(IOO) * 

1  JY(IOO) 

COMMON  M.M,A.B.H.E.SIGMA.RHO»PI»G 
READ (5. 9979)  SPEC 
9979  FORMAT (F10.0) 

A2=A*A 

R2=B*B 

A4=A2*A 2 

B4=B2*B2 

MPl=M+l 

NP1=N+1 

IF( SPEC.EQ.  1.0)  GO  TO  510 
GX ( 1 ) =1 « 

HX< 1 )=1. 

JXIU-1. 

GY( 1 )=1 « 

HY  ( 1 J  =»1« 

JY(1)=1. 

GX( 2  )  =  1.506 
HX ( 2 ) =1 • 248 
JX( 2  )  =  1.248 
GY ( 2 ) =1.506 
HY ( 2  )  =  1  *  248 
JY( 2  )  =  1 .248 


TABLE  8  (Continued) 


00  100  M1=3«KP1 
GXCMl}=FL0AT(Ml)-*5 

HX(H1)=( lFL0AT(Ml)“.5}**2)*<l.-2«/t  f FL0AT(Ml)-.5)*Pn ) 

JX(  M1 3  =HX(M1 ) 

100  CONTINUE 

DO  150  Nl=3tNPl 
GY ( Nl > =FLOAT C  Nl ) “• 5 

HY { N1 )  =  <  (FLOAT{ftH“«5)»*21*(l.“2./(  ( FLOAT (Nl  )-.5}*PIJ  1 
JY(N1 )=HYCNl J 
150  CONTINUE 
GO  TO  590 

510  DO  500  Ml  =  1*WPJ 

GX(M1)  s  FLOAT (Ml)  -1*0 
HX(Ml)  =  GXCM1I  **2 
500  JX(M1I  =  HX(M1J 

DO  550  N2  =  1*MP1 
GY(Nl)  =  FLOAT  t  fill -1.0 
HY(Nl)  =  GY(N1>**2 
550  JY(Nl)  =  HYJN1 ) 

590  WRITE(6*20‘A*B»H»E»SIGMA»RH0 

20  FORMAT C 1H1»3H  A=*F7.2t3H  B=*F7*2-3H  H=*F7.4*3H  E=»F11.4»7H  SIGMA=, 
1  F7.2.5H  RH0=»E11»4) 

WRITE(6»19J 

19  FORMAT t //23X»  22H  WARPURTON  FRFOuENCIFS) 

IW  =  1 

DO  400  N2=2»NP1 

N21=N2-1 

WR I TE ( 6  »  2 1 >  N2 1 

21  FORMAT (3H  N=*I2> 

WRITE(6»22J 

22  FORMAT ( 9X » 1 HM ♦ 1 5X  #6HLAMRPA 1 1 6X  *  5H  FREQ) 

DO  300  M2=2fMPl 

M21-M2-1 

XL AMSQ=GX (M2 ) *GX ( M2 ) *GX ( M2 ) *GX  CM2 ) + 1 «Y ( N2  J *GY ( N2 ) »GY( N2 )  *GY ( N2 ) 

1  *A4 ) /B4+ ( 2 • *A2 /B2 ) * ( S I GMA*HX C  M2 ) *HY ( N2  )  +  < 1. “SIGMA )  * JX ( M2 ) * JY ( N2  )  > 
LAMBDA* SORT ( XLAMSO) 

FRE0CM2 »N2 ) = ( ( LAMRDA*H*P I ) /A2 ) *SORT ( F  / (48«*RH0* ( 1 .-SIGMA**2 ) > > 

WRITE(6*23)M21»LAMBDA»FRFO(M2*N7J 

23  FORMAT ( 5X » 15 » 5X »E15. 8 »5X* El  5. 81 
OMEGA(M2»N2)  =  2«  *  PI  *  FRFO(M2»N2) 

WRITE(6»30)  OMFGA(M2»N2 ) *  IW 

WRITE(8*30|  OMEGA (M2  *N2 ) *  IW 

30  FORMAT! F10#4»65X» 15} 

IW  *  IW  +  1 
300  CONTINUE 
400  CONTINUE 
RETURN 
END 


TABLE  8  (Continued) 


SIBFTC  CRCKER 

SUBROUTINE  CROCK 
DIMENSION  FREQ 120 *10) 

COMMON  M»N*A»B*H»E»SIGMA*RHO»PI *G 
REAL  LAMBDA 

VRITE(6*4)A»B»H*E»SIGMA»RH0 

4  F0RMAT(1H1,3H  A=,F7.2*3H  B=.F7.2»3H  H*.F7.4.3H  E=.Ell.4*7H  SIGMA 
1  F7.2*5H  RH0*»F7.2) 

KRITE(6#19) 

19  FORMAT C//23X«20H  CROCKER  FREQUENCIES) 

DO  40  J*1»N 

GAMN=(2.»Fl0ATCJ)-s-1.)*PI/2» 

AN= ( COSH  C  GAMN  > -COS ( G AMN ) ) / { S I NH ( GAMN ) +S I N ( G AMN ) J 
WRITE(6»13)J 

13  FORMAT ( 3H  N=*I2) 

WRITE(6,14) 

14  FORMAT (9X*1HM»15X»6HLAMBDA»16X»5H  FREQ) 

Z I  N=  ( GAMN/2** C  t ( AN-1 • ) *C0SH t GAMN ) +AN* t -EXP ( -GAMN ) -SI N ( GAMN ) ) 

1  -COS ( GAMN } ) *  C ( AN-1 • 1 *S I NH ( GAMN  >  +AN* ( EXP l -GAMN ) +COS  C  GAMN ) ) - 

2  S I N ( GAMN ) ) +4. *AN ) -2. *GAMN**2 ) / 2. *AN*AN 
DO  30  1=1 »M 

GAMM= ( 2  «*FLOAT  1 1 )  +l«)*PI/2. 

AM=(COSH(GAKMJ-COS(GAMM) )/(SINH(GAMM)-SIN(GAMM) ) 

Z I M= ( GAMM/2 • * ( ( ( AM-1 • ) *COSH ( GAMM ) +AM*  t -EXP  C -GAMM ) -S I N  C  GAMM  > ) 

1  -COS  t  GAMM ) ) * ( ( AM-1 . ) *S I NH{ GAMM  >  +AM»  t  EXP  t -GAMM ) +COS ( GAMM ) ) - 

2  SINCGAMM) :+4.*AM)-2«aGAMM**2?/2**AM*AM 
LAMBDA* ( B*GAMM/A ) **4+GAMN**4+2 . *Z I M*Z I N* i 6/A } **2 

FREQ ( I • J ) =SQRT ( LAMBDA*E/ ( 12.*RHO» ( 1 .-SI6MA**2  > ) ) *H/B**2 
FREQ(I*J)=FREQ(I.J)/(2.*PI> 

WRlTEt  6»7) I »LAMBDA»FREQ( I » J) 

7  FORMAT (5X» I5»5X »E15»8»5X*E15»8) 

30  CONTINUE 
40  CONTINUE 
50  CONTINUE 
RETURN 
END 

SIBFTC  GRNSP 

SUBROUTINE  GREEN 

DIMENSION  FREQ(5»5)»P(5)»X<5) »YI5)*XSQ(5)»YSQ{5) 

COMMON  M»N»A«B»H»E»SlGMA»RHO»P I *G 

P(l)=4.73 

P ( 2 ) =7*8532 

P(3 ) =10*9956 

P(4)=14«1372 

P(  5 ) =17*2788 

X(l)=-12. 3026/A 

X<2)=-46« 0501/A 

X(3)=-98. 9048/A 
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TABLE  8  (Continued) 


X(4) =-171. 2560/A 
X(  51 =-263. 9980/A 
Y(l)=-12.3026/B 
Y( 2) =-46.0501/8 
YC3)=-98.9048/B 
Yl4) =-171.2560/8 
Yf  5) =— 263.9980/B 
DO  1  1=1*5 
XSQ(I)=A 
YSQ( I)=B 
1  CONTINUE 
A4=A**4 
84=B**4 
H3=H**3 

WRITE(6»8)A»B»H,E»SIGMA.RH0 

8  FORMAT ( 1H1 »3H  A=.F7.2»3H  B=*F7.2»3H  H=»F7.4.3H  E=.E11.4»7H  SI6MA 
1  F7.2»5H  RHO=»F7.2) 

D=E*H3/( 12.* ( 1 .-SIGMA**2 ) ) 

F=SQRT(D/CRHO*H) 1 
IF  (M  .GT.  5)  M=5 
IF(N  .GT.  5)  N=5 
WRITE! 6*19) 

19  FORMAT (//23X»22H  GREENSPON  FREQUENCIES) 

DO  20  J=1.N 

WRITE(6»4)  J 

4  FORMAT (///3H  N=»I2> 

WRITE(6»5) 

5  FORMAT (9X»1HM» 15X.5H  FREQ) 

DO  10  1=1 *M 

FREQ(I,J)=F*SQRf((P(I)**4/A4)+(PCJ)**4/B4)+(2.*X(I)*Y( J) )/ 

1  (XSG( I5*YSQCJJ) ) 

FREQ(I*J)=FREQ( I*Ji/!2-*PIl 
WRITE(6»6)  I»FREQ{I»J) 

6  FORMAT (5X*I5*5X*E15.8) 

10  CONTINUE 

20  CONTINUE 
30  CONTINUE 

RETURN 

END 


START 


Figure  24  -  Flow  Chart  for  WCGFRE.  Computer  Program  for  Computing  Natural  Frequencies 
of  a  Plate  by  Warburton,  Crocker,  and  Greenspon  Methods 

The  printed  output  of  the  program  contains  FREQ  (V,/V).  However  the  value  FREQ  (A/,  Ar)  x  2n  may  be  used 
as  the  input  OMEGA  (A/, /V)  to  Subprogram  A  in  Appendix  B  of  Reference  1. 
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Input  Description 

.  The  input  description  is  as  follows. 


Card 

No. 

Program 

Symbol 

Theory 

Symbol 

Description 

Units 

Format 

i 

IOPT 

OPTION  for  methods: 

1  —  Warburton  only; 

2  —  Cocker  only; 

3  —  Greenspon  only; 

4  —  all  methods 

15 

i 

NCASE 

Number  of  plates  to  compute 
frequencies  for 

15 

2 

M 

m 

Number  of  modes  in  x- 
direction 

15 

2 

N 

n 

Number  of  modes  in  y- 
direction 

15 

0 

< 

A 

a 

Plate  dimensions,  in¬ 
direction 

ft 

F12.6 

3 

B 

b 

Plate  dimensions,  y- 
direction 

ft 

F12.6 

3 

H 

h 

Plate  thickness 

ft 

F12.6 

4 

E 

E 

Young’s  modulus 

lb/ft2 

E16.8 

4 

SIGMA 

o  or  v 

Poisson’s  ratio 

F12.6 

4 

RHO 

Pw 

Weight  density  of  plate 

lb/ft3 

FJ2.6 

4 

G 

9 

Gravitational  constant 

ft/sec2 

F12.6 

5 

SPEC 

OPTION  for  W'arburton 
simply-supported  frequencies. 
Used  :  if  IOPT  =  1  or  =  4; 

SPEC  -  1.0  moans  simply- 
supported  case. 

FlO.O 

Cards  2—4  axe  repeated  NCASE  number  of  times. 


Output  Description 

The  input  data  and  results  are  labelled  and  printed  out  for  each  plate  (or  each  value  of 
NCASE).  The  first  printout  is  Waiburton,  followed  by  Crocker,  and  finally  Greenspon.  The 
mode  numbers  (oi,n),  nondimensional  frequency  A,  and  final  frequency  f  (in  hertz)  are  given. 

A  sample  problem  using  all  subroutines  to  compute  25  modes  each  for  two  plates  took  a 
total  of  1.1  minutes  on  the  7090. 


WHITE  (see  Table  9  and  Figure  25) 

White  has  provided  a  set  of  nomographs  that  permit  manual  computation  of  the  frequency 
parameters  am  n  =  \fXm  ~  for  the  first  nine  modes.  A  short  subroutine  handles  the  conversion 
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TABLE  9 

Program  Listing  for  WHITE  Computer  Program 


DIMENSION  FREO(20#7)*ALPHA<20*7) 

PI*3*14159?7 

WRITE<6«1) 

1  FORMAT (1H1,18H  WHITE  FREQUENCIES) 

RFAO<5*2>  NCASF 

2  FORMAT! 15) 

4  FORMAT (2 15) 

5  FORMAT !4Fl? .6) 

$  FORMATt£16.8i<2Fl2*6) 

T  FORMAT ( //3H  A*#F8*3»3H  R**F8»3#3H  H««F8*3»3H  E*»E11«4#7H  SIGMA»* 
1  F7*2#5H  RHO«#F8«3) 

9  FORMAT ! 9X  *1 HM# 1 5X *6HALPHA  #16X*5H  FRFO) 

8  FORMAT  (3H  N«*-l2) 

10  FORMAT (5X»I5*5XtE15«8»5X«E15*8) 

M  *  3 
N  *  3 

DO  40  L*1*NCASE 

RFAO(5»3)  (( ALPHA <I*J) *I«1*3) *J«1*3) 

3  FORMAT! 3F12.6) 

RFAO(5*5)  A*P»H  «G 
RFAD!5*6)  FtSIGMAiRHO 
WRITE! 6»7)  At.R*H*E»SIGMA »RHO 
A4  *  A**4 

B4»R**4 

H3«H*#3 

D=F*H3 /! 12«#( 1«-SI GMA*»2 ) > 

F=SORT! !D*G)/(RH0*H*A4) ) 

DO  30  N2«1»N 
WRITF(6*8)  M2 
WRITF(6»9) 

DO  20  M2«1»M 

FRFQ(M2»N2)»ALPHAfM2»N2>*F/t2«*PI) 

WRITE!6»10)  M2 *ALPHA  »FRFQ!M2»N2) 

20  CONTINUE 

30  CONTINUE 
40  CONTINUF 

STOP 

END 
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Figure  25  —  Flo«\  Chan  for  WHITE.  Computer  Program  for  Converting  Nomograph 
Frequency  Parameters  a  _  to  Frequencies  / 


The  printed  output  includes  FR£Q(V,A).  However,  the  value  2  Jr  >  FRHQ  (I/,  V)  may  be  used 
as  the  input  OMEGA  to  St.bp.'ogras  A  in  Appendix  B  of  Reference  !. 
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of  these  frequency  parameters  to  hertz  using  a  formula  given  by  White  (Appendix  E).  The 

b 

nomographs  are  read  for  various  aspect  ratios  —  <  1.  Thus  the  user  must  make  adjustments 
b  a 

for  the  case  —  >  1,  e.g.,  interchanging  m  and  n.  The  nomographs  are  applicable  to  nine 
a 

combinations  of  m  =  1,  2,  3  and  n  =  1,  2,  3. 


input  Description 

The  input  description  is  as  follows. 


Card 

No. 

Program 

Symbol 

Theory 

Symbol 

Description 

Units 

Format 

1 

NCASE 

Number  of  plates 

15 

(There  are  NCASE  sets  of  remaining  cards.) 

2-4 

(ALPHA)  (I,J), 
(1=1,3), 

(J=l,  3) 

Modal  frequency 
parameter,  found 
from  nomographs 

3F12.6 

5 

A 

a 

Dimension,  ^-direction 

ft 

F12.6 

0 

B 

b 

Dimension,  y-direction 

ft 

F12.6 

5 

H 

h 

Plate  thickness 

ft 

F12.6 

5 

G 

9 

Gravitational  constant 

F12.6 

6 

E 

E 

Young’s  modulus 

mmm 

E16.8 

6 

SIGMA 

■ 

Poisson’s  ratio 

F12.6 

6 

RHO 

Plate  weight  density 

lb/ft3 

F12.6 

Output  Description 

Both  ALPHA  and  FREQ  ( fm  n)  are  given  according  to  mode.  The  7090  computer  time 
is  about  30  seconds. 

PLFREQ  (see  Table  10  end  Figure  26) 

PLFREQ  is  a  computer  program  developed  by  Plumblee28  and  Ballentine19  to  yield 
the  natural  frequencies  of  vibration  of  either  a  simply  supported  or  clamped  thin  plate,  flat  or 
curved.  The  original  program  was  in  nondimensional  form.  However,  for  the  comparison  pur¬ 
poses  of  this  report,  the  program  was  modified  so  that  additional  input,  in  units  permitted  the 
frequency  to  also  be  computed  in  hertz. 

The  mathematical  subroutines  needed  from  the  IBM  SHARE  library'  are  EIGEN,  LOC, 
and  MINV.  The  sample  problems  for  36  modes  were  run  on  the  IBM  360/91  and  took  18  sec¬ 
onds  per  plate. 
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TABLE  10 

Program  Listing  for  PLFREQ  Computer  Program 


3k 

S 


T 

■? 


REAL  BETAL ( 20 ) ,M1 (20, 20 ) ,M2 ( 20, 20 ) , R ( 27 ) , ML,NU 
DOUBLE  PRECISION  L ( 378 ) , VECTOR ( 729 ) , VEC( 27 5 , XX 

DOUBLE  PRECISION  FR(5)  , 

INTEGER  LR(45),LM(45),P,Q,PP»QQ,QQQ,S,T,PQ,QI  i 

READ( 5,415)  RH0,AL,B,G,E  S 

415  FORMAT (4F12.6,E16.8) 

3  READ( 5,1 )TH£TA ,TL , A ,NU  : 

READ ( 5 , 2  ) MM , NN ,MV,LL, L BOUND 

1  FORMAT (4E1 0.4) 

2  FORMAT (5 12)  ' 

WRITE(6,15)THETA,TL,A,NU 

15  FORMAT (4X, •THETA=‘ , F10.4, 'TL=',F10.4, ' A=* , F10.4, * NU= * , F10. 4 ) 

W  R I T  £  { 6 , 1 6 )  MM ,  N N ,  M  V ,  L  L  ,  L  B  0 U  N  D 

16  FORMAT (4X, *MM= *, 12, ' NN  =  * , 12, »MV=' ,12, 1 LL=*, 12, r LB0UND= ' , 12) 

R(1)=LB0UND 

CALL  B ET A ( MM , NN , R , B ETAL , M2 , M 1 ) 

IF(MM-NN)  41,41,42 

41  KK=2-NN 
GO  TO  43 

42  KK=2*MM 

43  WRITE(6,46) 

DO  44  1=1, KK 

WRITE (6,48)  (M1(I,J),J=1,KK) 

44  CONTINUE 
WRITE(6,47) 

DO  45  1=1, KK 

WRITE(6,48)  (M2 ( 1 , J ) , J=1 ,KK) 

45  CONTINUE 

46  FORMAT( 1H1,4X, 'MATRIX  Ml(I,J)',//> 

47  FORMAT! 1H1 ,4X, 'MATRIX  M2(I,J)',//> 

48  FORMAT ( 5X,9E12.5 ) 

MN=MM*NN 

MN5S3SMN 

P=1 

GO  TO  11 


107 


•M*'1 


TABLE  10  (Continued) 


10  P=P+1 

11  CALL  SUBSCP(P,MN,NN,LL,PP,S,T) 

Q=P 

GO  TO  IB 

12  0=Q+1 

IB  CALL  SUBSC  P ( Q  »  MN » NN , L  L » QQ » M , N ) 

CALL  LOC(P,Q,PQ,MN5,MN5,l) 

GO  TO  ( 101 , 102, If  3 ) 5  P° 

101  GO  TO  ( 1011 ,1012,1013) ,QQ 

1011  L(PQ)=A*BE7AL(M)**3*MliS,M)*Ml(T,N)/BETAL(S)+( l.-NU)*M2(S,M) 
1*M2(T,N)/(6ETAL(M)*BETAL(S)*2.*A) 

I F { P— Q )  12,10111,12 

10111  R(P)=M2(S,M}-M1(T,N)/(8ETAL($)*BETAL(M>> 

GO  TO  12 

1012  L(PQ)  =  (1.+NU)=!=H2(S5M)*H2(T,N)/(BETAL(S)*BETAL(N)*2.0) 

GO  TO  12 

1013  L  {  PQ  i  =  -NU*THETA*M2(S,M)=?M1<T,N)/BETAL{S) 

I F ( 3=S=MM#NN-Q ) 10 , 10 , 1 2 

102  QQQ=QQ-1 

GO  TO  (1022, 1023), QQQ 

1022  L(PQ)  =  M1($,K)*K1(T,N)-~BETAL(N)**3*(  1 . +THETA**2/ ( 12. * < TL*A) **2 )  ) 

1  /( A*BETAL(T) ) -Ml .-NU)*A*M2 (S,M)*K2 (T,N )*(!.  +  ( (THETA/A/TL)**2/3.0) 

2  )/(2.*BETAL(T )*BETAL(N)  ) 

IF(P-Q)12, 10222, 12 

10222  R ( P ) =H1 ( S ,M )*M2 (T ,N )/ ( BET AL ( T)*BETAL ( N) ) 

GO  TO  12 

1023  L(PQ)=THETA*M1(S,M)*M2(T,N)/(A*BETAL(T)  )+THETA*M2(S,M)*M2( T,N) 
l*f\'U/(  12.*AvTL=!:TL:5:BETAL(N))+THETA*)';i(S,H):!:Ml(T,N)*BETAL(N)*=!:4/12. 
2/(TL*TL-A**3=f=SETAL(T  ) )  +( 1  «-NU  )*TKETA*M2(  S,M)  *M2(  T,  N)  /  (  6.*A*TL*TL) 
3/SETA  L( T ) 

IF(MN5-Q)10, 10,12 

103  L(PQ)=THETA=i=*2:!:Ml(S,M)«Ml(T,N)/A+A-~HllS,M)*BETAL(M):5:-4*Ml(T,N) 
l/(12.*TL*TL)+Kl(S,K)*MlCT,f  -BETAL(N)**4/(12.*TL-TL*A=!=-3)  + 
2M2(5,M)*M2(T,N)/(6.*TL*TL*A) 

1F(P-Q)1033, 10333, 1033 
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TABLE  10  (Continued) 


10333  R(P)=M1(S,M)*M1(T,N) 

1033  1 F { MN5-Q )  1034,1034,12 

1034  IF(MN5-P)100,100,10 
100  00  110  1=1, MN5 

110  R( 1 )=SQRT ( R ( I ) ) 

00  120  1=1, MN5 
DO  120  J=I,MN5 
CALL  L0C(I,J,IJ,MN5»MN5,1) 
L{IJ)=L(IJ)/(RU  KRIJD 
120  CONTINUE 
0FACT=1. 

140  DN0RM=1. 

00  150  1=1, MN5 
CALL  L0C(I»I,II,MN5,MN5,1) 
DN0RM=DN0RM*L{II I/DFACT 
IF (0N0RM-1.D+70 1145,155, 155 
145  1F(0NGRM-1. 0-70)160, 160, 150 
150  CONTINUE 
GO  TO  165 

155  DFACT=10.*DFACT 
GO  TO  140 

160  OFACT  =0.1*DFACT 
GO  TO  140 

1 65  DNORM=  ( ABS  ( DNORM ) ) ** ( 1  /MN5 ) 

DO  170  1=1 ,MN5 
DO  170  J= I ,MN5 
CALL  LOCH  ,J,IJ,HN5,MN5,1) 

170  L(IJ)=L(IJ)/ (DNORM^OFACT ) 

00  125  1=1, MN5 
DO  125  J=1,MN5 
CALL  LOCI  I ,J» 1K,MN5,MN5,0) 

CALL  LOC(I,J,IJ,MN5,MN5,l) 

125  VECTOR! IK)=L{IJ) 

MN52=MN5*MN5 

CALL  MINV ( VECTOR, HN5 ,XX ,LH ,LR ) 
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TABLE  10  (Continued) 


WCITE ! 6, 130 )  XX 

130  FORMAT {  «0»,'THE  DETERMINANT  IS',  E12.53 
DO  135  1=1, MN5 
DO  135  J= 1 ,MN5 
CALL  LOC ( I , J , I J ,MN5,MN5, 1 > 

CALL  LOC ( I , J, IK,MN5,MN5,0) 

135  L  (  1 J )= VECTOR { IK ) 

CALL  EIGENIL, VECTOR, MN5, MV) 

20  FORMAT !  '  1* ,8X, '0IM6NSIQNLESS  FREQUENCIES  ARE  NORMALIZED1, 

1  2X,« EIGENVECTORS') 

WRITE(6,20) 

21  FORMAT ( 33X , ' FOR  '  ) 

WR1TE!6»21) 

22  FORMAT ( 21X, 1 A  CYL INDRICALLY  CURVED  PANEL') 

HRITE(6,22) 

23  FORMAT { 32X, 'WITH* ) 

WRITE! 6,23) 

GO  TO  l 241 , 242 ) , L30UND 

241  WRITE! 6,24) 

24  FORMAT! 28X, 'CLAMPED  EDGES') 

GO  TO  251 

242  WRITE!6,245) 

245  FORMAT (23X, 'SIMPLY  SUPPORTED  EOGES*) 

251  WRITE! 6, 25) 

25  FORMAT! 'O' ,29X, •*#$#******• ) 

26  FORMAT! 'O' ,19X, 'NONDIMENSIONAL  INPUT  PARAMETERS') 

WRIT  E ( 6,26) 

27  FORMAT! 'O' , 'SUBTENDED  ANGLE= •  F7.4, 10X, 'ASPECT  RAT  1 0= » , F7„ 4 
WR1TE!6,27)THETA,A 

28  FORMAT! 'O' , 'LENGTH/SKIN  THICKNESS=' ,F7.2) 

WRITE (6, 28)  TL 

WRITE ! 6 , 29 )  NU 

29  FORMAT! 'O', 'POISSONS  RATIO=' ,F4.3) 

32  FORMAT!  *0«  , 'NUMBER  OF  SERIES  TERMS  ALONG  STRAIGHT  EDGE=',I). 
1», ALONG  CURVED  EDGE=',I1) 
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TABLE  10  (Continued) 


WRITE ( 6, 32 ) MM » NN 

33  FORMAT!  'O',  29X,  •******=!=***'//,  17X  *  COMPUTED  FREQUENCIES  AND', 
1  *  MODE  SHAPES'  ) 

WRITE(6»33) 

00  180  1=1, MN5 

CALL  L0C(I,I,II,MN5,MN5,1) 

IF(L( II) 1180,180,179 

179  L  ( I )  =0  . 1 59154=?SQRT  ( ONORM^OF  ACT  J/DSQRT  {  L<  II )  ) 

180  CONTINUE 
11  =  1 

60  TO  51 

50  11=11+1 

51  MI  =  5=S=(II-1)+1 
NI=5*II 

IF(MN-4) 520,520,523 

520  GO  TO  (521, 521, 522, 523), MN 

522  GO  TO  (521,531,521) , II 

523  IF(II-1)521,521,533 

532  FORMAT( *  1 «////////// ) 

533  WR ITE ( 6,532 ) 

GO  TO  521 

53  FORMAT ( '1'  ) 

531  WRI TE ( 6 , 53 ) 

52  FORMAT ( 'O' , 'FREQUENCY=',5(1X,E11.4)) 

521  WRITE (6,52  )  ( L ( I ) , I =M I ,N I ) 

J  =  1 

DO  5521  I  =  M I , N I 

FR(  J  )  =  L  ( I  )*  SQRT  (  (  E*G )  /  (  RHO^AL^BS  (  l.-NU*^ ) )  ) 

5521  J  =  J  +  1 

WRIT  E ( 6,5522  )  (FR(I),I  =  1,5) 

5522  FORMAT ( 10X,5(1X,E11.4)  ) 

54  FORMAT( 'O' , 'GEN  COORD *, 3X , 5 ( 2X ,' MODE  SHAPE')) 

WRI TE (6 , 54 ) 

Q=1 

GO  TO  61 
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TABLE  10  (Continued) 


i 


60  Q=Q+1 

61  CALL  SUBSCP(Q,MN,NN,LL,QQ,M,N) 

GO  TO  (7110, 7210, 7310), QQ 

7110  DO  711  I=MI,NI 

CALL  L0C(Q,I,QI,MN5,MN5,0) 

711  VEC(I)=VECTOR(QI) 

WRITE(6,71)M,N, ( VEC (I ) , I =MI ,NI ) 

GO  TO  60 

7210  DO  721  I=MI,NI 

CALL  L0C(Q,I,QI,MN5,MN5,O) 

721  VEC{ I ) = VECTOR ( QI ) 

WRITE ( 6, 72 )M,N» ( VEC ( I ) , I=MI »NI ) 

GO  TO  60 

7310  DO  731  I=M I ,N I 

CALL  LOC(Q, I ,QI,MN5,MN5,0) 

731  VEC( I )= VECTOR (QI ) 

WRIT E( 6,73 )M,N, (VEC (I ),I=MI,NI ) 

I F ( MN5-Q )76 ,76, 60 

76  IF(MN5-NI )77,77,50 

77  WRITE (6,53) 

80  CONTINUE 

IF ( LL-4)  3,74,74 

71  FORMAT (2X, 'U( * , II »  » , * , 1 1, • ) • ,4X,5 { IX, E11.4) ) 

72  FORMAT! 2X, *V(  Ml,',  *,I1»'  )  ‘  ,4X ,  5  ( IX,  E 11 .4)  I 

73  FORMAT (2X, • W ( * ,11, •, »,I1, • ) • ,4X,5( IX, Ell. 4) ) 

74  CONTINUE 

APL=SQRT (41 . 7*A+25 .2/A+41 TL*THETA) «*2/A) 
WRITE (6,78)  APL 

78  FORMAT { Ell »4) 

STOP 

END 

SUBROUTINE  BETA (M ,N , A , B ,G,H ) 

DIMENSION  A(1),B(1),G(20,20),H(20,20) 

IF(M-N) 1,1,20 
1  KK=2*N 
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TABLE  10  (Continued) 


GO  TO  2 
20  KK=2*M 

2  IF(A(1)-1. 519,9,10 
9  DO  5  1=5, KK 

IF( 1-5)4, 4, 3 

3  A ( I ) =1 .0 

B(I  )=(2*I+1)*1. 5707963 
GO  TO  5 

4  A(l)=. 9825022158 
A(2)=l. 000777311 
A(3)=. 9999664501 
A<4)=1. 00000145 
A( 5)=. 9999999373 
B{ 1 ) =4 . 7300408 
B(2)=7. 8532046 
B(3)=10. 9956078 
B(4)=14. 1371655 
B(5)=17. 2787596 

5  CONTINUE 

DO  8  1=1, KK 
DO  8  J= 1 , KK 
1F(I-J)7,6,7 

6  G  (  I  ,  J  )  =  A  { I  )  *  B  ( I  )  *  ( A  (  I  )  *  B  (I  )  -  2 . 0  ) 

H ( I , J ) =1 .0 

GO  TO  8 

7  G ( I , J ) =-4.*B ( I  )**2*B(J)**2*(A(I )*B<I ) -A( J )*B ( J ) ) * 
1  ( 1 .  +  ( -1 . ) ** ( I +J } ) / ( 3 ( I ) **4-B ( J ) **4  ) 

H ( I , J ) =0 . 0 

8  CONTINUE 
RETURN 

10  DO  11  1=1, KK 

B( I )=I*3. 1415927 
DO  11  J= 1 , KK 
I F ( I- J ) 12, 13, 12 
12  G( I ,J)=0.0 


■  " 
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TABLE  10  (Continued) 


H(I,J)=0.0 
GO  TO  11 
13  G(I,J)=e(I)**2 
H ( I , J}=1.0 
11  CONTINUE 
RETURN 
END 

SUBROUT INE  SUBSCP(NR»MN,NN,KK,NP, J,K) 

NP=( (MR-1 )/HN)+l 

1=NR-(NP-1 )«MN 

II  =  (  I-D/NN 

GO  TO  (1,2,3,4),KK 

1  J=2* I 1+1 
K=2*I-2#II*NN-1 
RETURN 

2  J=2*II+2 
K=2*I-2*II*NN-1 
RETURN 

3  J=2*II+1 
K=2*I~2*II*NN 
RETURN 

4  J=2*II+2 
K=2*I-2*II*NN 
RETURN 

END 


»***r%sf  -vi-  +**;?■*» ^ *WV» j^^r‘wywn^«nw:y,VF«tt»i?y*flfg^;yr'»"imm<?3^^ 


hi 

UmQ 

Z|-Q 

£S:9 

£  DC 

<UQ 

!  ••  COo 

Q-ra 

Uv>  — 

S«2  11 

l|<H 

i  CO  — I 


zz 

Q  IUILI 
Q>> 
OUJUi 

■  •  a 

z  z 

lUOUJ 

>o> 

UJ  o  til 

NM'tf 

II  II  II 


_I  —J  — J  — J 


-<ZN 

HlilJ- 

UiOjU 

>  lll£ 

-Z.CH 
OIL  = 


lilz  < 

3°0' 

Shi-! 

iq2i 

OZ^a' 

Uo°; 


llio  LU  oo 

>  II 

-I  _JO 
o>  <f- 
>  u 

..<x  Z  HI 
2  ^  Ui> 

LU  I  UZ 

<2  “Wo 

«url 


q  tr  < 

<?h 
UJ  Eb  < 
q;Zq 


.t  v> 

fiM  HI  I™1 

:n£9= 

<“bo 

;sp 

l£<=> 

z 


Q  «/i 
<  LU  UJ 
H  M  \- 

w  2*  J0  _i  :$ 

h^'*'<£ 

S  I—  '0£Q 
^  =}  J®  UJ  oc 

0.  ^  z  o 

Z  UJ  O 

—  o  u 


115 


Input  Description 


The  input  description  is  as  follows. 


Card 

No. 

Program 

Symbol 

Theory 

Symbol 

Description 

Units 

Format 

■ 

RHO 

Pw 

Plate  weight  density 

F12.6 

1 

AL 

a  at 

Panel  length 

ft 

F12.6 

1 

B 

b 

Panel  arc  length 

ft 

F12.6 

1 

G 

9 

Gravitational  constant 

F12.6 

n 

E 

E 

Young’s  modulus 

lb/ft2 

E16.8 

For  each  value  of  LL,  there  is  a  set  of  the  following  cards: 

2 

THETA 

6 

b 

Subtended  angle  —  (0  for  flat  plate) 

R 

E10.4 

2 

TL 

t 

h 

If  curved  panel,  R= panel  midplane 
radius,  ratio  of  panel  length  to 
thickness 

ElO.4 

2 

A 

b 

T 

Aspect  ratio 

E10.4 

2 

NU 

V 

Poisson’s  ratio 

E10.4 

3 

MM 

m 

Modes,  ^-direction 

12 

3 

NN 

n 

Modes,  y-direction 

12 

3 

MV 

0  eigenvalues  and  eigenvectors 

1  eigenvalue  only 

12 

3 

LL 

1  odd-odd  modes 

2  even-odd 

3  odd-even 

4  even-even 

12 

3 

LB013ND 

1  clamped  edges 

2  simply  supported  edges 

12 

Output  Description 

The  frequencies  are  printed  out  in  ascending  order  for  each  set  of  subscripts  (odd-odd, 
even-odd,  odd-even,  even-even).  The  nondimensional  frequency  is  given  first,  with  frequency 
in  hertz  on  the  next  line.  The  generalized  coordinates  and  mode  shapes  are  also  given  in  the 
same  column  as  the  frequencies  they  represent. 

SUNFRE  (see  Table  II  and  Figure  27) 

SUNFRE  is  a  computer  program  developed  by  Sun24  to  obtain  the  natural  frequencies 
of  vibration  of  a  class  of  thin  plates,  including  such  special  cases  as  the  circle,  square,  and 
rectangle. 
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TABLE  11 

Program  Listing  for  SUNFRE  Computer  Program 


C  FREQUENCIES  OF  GENERAL  PLATS  By  RITZ  METHOD 

DOUBLE  PRECISION  XHAI21.21 > .XHBI21 .21 ) ,XKC( 21 .21 ) »XI (43) »YI (48) 
DOUBLE  PRECISION  Hi (48) »HOR( 21) »VER(21> »AREA{ 462) »AREAU<462) 
DOUBLE  PRECISION  AREAV(462) »XWW(462) 

DIMENSION  XP ( 2 1 ) »  YP(21) 

DOUBLE  PRECISION  XUI21.21 ) »XMD(21 .21 ) ,A( 21) ,B(21 ) *C(21 J 
DOUBLE  PRECISION  VAL.blG.AdSD.P.COHV.AKPLTD.EIGENS 
DOUBLE  PRECISION  VXPI 21 ) *VYP ( 21) 

COMMON  XMA,  XKB,  XKC.  XI.  Yl.  Wl»  HOR»  VER,  AREA,  AREAu.  AREAV, 

1  XWW, P.B. ALPHA .BETA .RAT 10. NX .NROW.XP.YP, AMI »BK1» 

2  SWITCH,  VXP,  VYP 

READ  (5.  999  )  NX,  (XI(I).  1=  1,  NX  ),  (WI(I),  I  =  1,  NX  ) 

999  FORMAT (I 10  /  (4E20.10)) 

DO  2  I  =  1.  NX 
2  Yl ( I )  =  XI ( I ) 

SWITCH  =  0. 

10  READ  (5.  1000)  ALPHA,  BETA,  RATIO,  MODE,  NOIT,  NP,  LIMIT,  CONV 
1000  FORMAT  (  3F5.2.  415,  F10.7  ) 

LAST  *  0 

MODE  =1  X,  Y  TAXE  EVEN  POWER 
MODE  =2  X.  Y  TAXE  ODD  POWER 
MODE  =  3  X  TAXE  EVEN  POWER,  Y  TAXE  ODD  POWER 
MODE  =  4  Y  TAKE  EVEN  POWER,  X  TAXE  ODD  POWER 
NOIT  =  NUMBER  OF  EIGENVALUES  DESIRED 
NP  =  0  MO  POINTS  FOR  NODAL  LINES 
NP  *  20  20  POINTS  FOR  NODAL  LINES  PLOT 

LIMIT  =  800  (RECOMMENDED)  CYCLES  OF  ITERATION 
CONV  =  0*00001  IS  RECOMMENDED 
CALL  XPYP  ( XP. YP »NROW .MODE) 

WRITE  (6,1050)  ALPHA.  BETA.  RATIO,  NROW,  MODE 
1050  FORMAT ( /  2X,  7HALPHA  =,  F6.2.8H  BETA  =,  F6.2.9H  RATIO  *,  F6-2. 
1  AX,  25HNO.  OF  TERMS  IN  X  AND  Y  n  ,14,  8H  MODE  *  ,  1 3  ) 

WRITE  (6,1052)  <  ( XP  ( I )  ,  YP(I)J,  I  =  j,  ,\R0U  ) 

1052  FORMAT  (  7(2H  (,  F3.0,  F3.0,  2H)  !  ) 

P  =  1.  /  (RATIO  **  BETA  ) 

AMI  =  ALPHA  -  1. 

BM1  =  BETA  -  1. 

CALL  DUBINT 
ICCT  =  1 

DO  12  I  =  1,  NROW 
DO  12  J  =  I  ,  NROW 
XMCII.J)  =  AREA! ICCT) 

XMC(J.I)  =  AREA ( ICCT ) 

12  ICCT  =  ICCT  +  1 

DO  13  I  *  1.  NROW 

13  WRITE  (6.  1054)  (XMCII.J),  J  =  1,  NROW  ) 

1054  FORMAT  (//(IX,  3D25.16  )) 

DO  14  I  =1,  NROW 
DO  14  J  =1,  NROW 
XMA ( I »  J  =  AREA( ICCT) 

XMAIJ.Ii  =  AREA! ICCT) 

14  ICCT  =  ICCT  +  1 

DO  15  I  =  1,  NROW 

15  WRITE  (6,  1054)  (XMA(I.J),  J  =  1,  NROW  ) 

IF  (  NROW  -  1  I  16,  16,  18 

16  AMPLTD  =  XMC(l.l)/  XMA(l.l) 

El  GENS  =  DSQRT (  AMPLTD) 

WRITE  (6.  1060)  El GENS 

1060  FORMAT  ( //  3X»  15HEIGEN  VALUE  =  ,  D25.16  //) 

GO  TO  10 


0010 

0020 

3C 

40 

0050 

0060 

70 

0080 

0090 

100 

0110 

120 

130 

140 

150 

0160 

0170 

180 

0190 

0200 

0210 

0220 

0230 

0240 

0250 

0260 

0270 

0280 

0290 

0300 

0310 

0320 

0330 

U34u 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

0450 

460 

470 

480 

490 

500 

510 

520 

0530 

540 

550 

560 

570 

0580 

590 
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TABLE  11  (Continued) 


18 


20 


CALL  SMTHX  t  XHC.  XHA.  NROW.  XMB.  XU  ) 

WRITE  (6*  1070)  ((XKB(I.J).  J  *  1»  NROW  )»  I  *  1*  NROW  ) 
DO  20  I  *  1*  KROW 
DO  20  J  =  I»  NROW 

XM8( I*  J)  *  (XMBII.J)  +  XMBl J.I) )/2. 

XMB-IJ.I)  *  XMBl  I  »J) 

WRITE  (6*  1070)  ( (XMB( I.J)  »  J  *  !♦  NROK)»  I  =  1*  NROW  ) 


1070  FORMAT  (IX.  5D25.16  ) 

CALL  EIGEN  (  XMB.  NROW*  NOIT.  A.  XMD.  LIMIT.  CONV.  TELL.  NUMCYC  ) 

C  A  *  COLUMN  MATRIX  OF  EIGENVALUES 

C  XMD  -  SQUARE  MATRIX  OF  CALCULATED  EIGENVECTORS  FOR  MATRIX  PENCIL 
WRITE  (6.  1072)  TELL.  CONV.  LIMIT  •  NUKCYC 
1072  FORMAT ( /  2X.6HTELL  =.  F5.2.  3X»  20HC0NVERGENCE  FACTOR  =  .  F10.8. 

1  3X»  15HLIMITED  CYCLE  =  .  15  /  3X.  19HNUMUER  OF  CYCLE  = 

2  .  16  ) 

IF  (  TELL  )  10.  10.  30 
30  CONTINUE 

DO  AO  I  =  1.  NO IT 

40  All)  =DSORT  (  1.  /  AID)  *  4.0 

WRITE  (6.  1076)  (A(I).  I  =  1.  NOIT  ) 

1076  FORMAT  (IX.  16HEIGENVALUES  ARE  .  /  (5D25.16)  ) 

DO  44  I  =  1.  NOIT 

44  WRITE  (6.  1078)  I.  (XMD(I.L).  L  *  1»  NROW  1 
1078  FORMAT  (3X,  13.  31HTH  EIGENVECTORS  FROM  ITERATION  /(  5D25.16  )) 
NM1  =  NOIT  “  1 

DO  48  I  =  1.  NM1 
IP1  =1+1 

DO  48  J  =  IP1»  NOIT 
VAL  =  0. 

DO  46  K  *  1.  NROW 
46  VAL  =  VAL  +  XMD(I.K)  *  XMOIJ.K) 

48  WRITE  (6.  J080)  I,  J.  VAL 

1080  FORMAT  (  ?X.  14  *  25HTH  EIGENVECTORS  MULTIPY  .14.  25HTH  EIGEN  V 

1ECT0RS  EQUAL  TO  .  D25.16  > 

DO  70  I  =  1.  NOIT 

52  DO  53  J  =  1.  NROW 

53  C(J)  =  XMD ( I . J  ) 

CALL  TRAVEC  (  XU.  C»  B.  NROW  ) 

C  B  -  ORIGINAL  COLUMN  MATRIX 

BIG  =  0. 

DO  56  J  *  1»  NROW 

ABSB  =  DABS(BIJ)  J 

IF  (  BIG  -  ABSB  )  54.  56.  56 

54  BIG  =  ABSB 
56  CONTINUE 

DO  60  J  *  1.  NROW 
60  B(J)  ■  B( J)  /  BIG 

WRITE  (6.  1090)  I*  All).  IB(J).  J  =  1.  NROW  ) 

1090  FORMAT  (  2X»  12.  15HTH  EIGEN  VALUE  .  D25.16  /(  /5D25.16  )) 

IF  l  NP  )  66.  70.  66 

66  CALL  PLNODE  l  NP  ) 

70  CONTINUE 

LAST  *  LAST+  1 

100  IF  (  LAST  -  1  )  10.  300.  300 
300  CONTINUE 
STOP 
END 


0590 

0600 

610 

620 

0630 

640 

0650 

660 

0670 

0680 

0690 

0700 

0710 

0720 

730 

740 

750 

760 

0770 

0780 

0790 

800 

0810 

0820 

830 

840 

850 

860 

870 

880 

0890 

900 

0910 

0920 

930 

940 

950 

960 

970 

980 

990 

100C 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 


SUBROUTINE  XPYP  (XP.  YP.  NROW.  MODE 


1160 
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TABLE  11  (Continued) 


1 *NROW) 


DIMENSION  XP 1 21 ) »YP(21 ) 
READ!5»1000>  NROW 
1000  FORMAT  !I10> 

DO  1110  H=  l»MODE 
1110  READ(5»1100)  lCXP(I)»YPm>»  I 
1100  FORMAT  (16F5.2) 

RETURN 

END 


SUBROUTINE  DUBINT 

DOUBLE  PRECISION  XMAI21.21) »XMB!21 *21 > *XMC! 21*21 > »XI !48 >  »YI 148) 
DOUBLE  PRECISION  Wl (48) »HOR 1 21) *VER (21 ) *AREA( 462) »AREAU 1 462 1 
DOUBLE  PRECISION  AREAV( 462 ) »XWW( 462 J 
DIMENSION  XP 1 21 ) »  YPI21) 

DOUBLE  PRECISION  HXX! 21 ) *HYY ( 21 ) »HXYt 21 ) «B>21* 

DOUBLE  PRECISION  BOQ.WII.UI fVI»DU.DV.rfIJ.YPS.YMS»YUP»YUM*YVP 
DOUBLE  PRECISION  YVM.XWIJ.P 

COMMON  XMA*  XMB *  XMC»  XI*  Yl»  Wl*  HOR»  VER*  AREA*  AREAU*  AREAV* 

1  XWW.P.B.ALPHA. BETA. RATIO. NK*NR0W*XP.YP»AM1»BM1 

SMI  =  .667 

NO  =  NROW*  ( MROVI ,  +  1) 

BOO  =  1.  /  BETA 
DO  1  K=1»N0 
AREAU(K)  =  0. 

AREAVIK)  =  0. 

1  AREAIK)  *  0. 

DO  20  1*1 »NK 
WRITE  (  6*  1000)  I 
1000  FORMAT  (  3X*  3HI  =.  13  ) 

WII  =  Will) 

UI  *  0.5*(1.+XI(I)) 

VI  =  0.5*U.-XI!I>) 

DU  =  RATIO*! (l.-Ul**ALPHA>**BOQ} 

DV  =  RATIO*! !1.-VI**ALPHA)**B0Q> 

DO  14  J=1 *NK 
WIJ  =  WIIJ) 

YPS  =  0*5*11. +YI(J>> 

YMS  =  0.5*!1.-YI!J)) 

YUP  =  DU*YPS 
YUM  =  DU*YMS 
YVP  =  DV*YPS 
YVM  =  DV*YMS 

CALL  ALL  IUI.YUP*  HXX*  HYY •  HXY  ) 

IC  =  1 

DO  4  KJ*1»NR0W 
DO  4  KI=KJ*NRGW 

XWW(IC)  =  HOR(KJ)  *  HOR(KI)  -  SMI  *  <HYY!KI>  *  HXX(KJ) 

1  +  HXX!KI)  *  HYY!KJ)  -  2.  *  HXYIKI)  *  HXY(KJ)  ) 

a  IC  =  IC+1 

DO  5  KJ  =  1*  NROW 
DO  5  KI  *  KJ.NROW 
XWW !  IC )  *  VEROCJ)  *  VER (XI ) 

5  IC  *  IC+1 

CALL  ALL  (  UI*  YUM*  HXX.  HYY*  HXY  ) 

IC  =  1 

DO  6  KJ*1 *NROW 
DO  6  KI=KJ.NROW 

XWIJ  *  WIJ  *  ! XWW ( I C )  +  H0R1KI )  *  HOR(KJ)  - 

1  SMI  *  <  HYY(KJ)  *  HXX(KI)  +  HXXtKJ)  *  HYY(KI)  -  2*  *  HXY!KI)  * 
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TABLE  11  (Continued) 


2  HXY(KJ)  )  ) 

AREAU(IC)  =  AREAU(IC)+XWIJ 

6  IC  =  IC+1 

DO  7  KJ  =  1»  NROW 
DO  7  KI  *  KJ*  NROW 

XWIJ  =  WIJ  *  (XWWtIC)  +  V£R(KI)*VER(KJ)> 

AREAUUC)  =  AREAU(  IC)+XWI J 

7  IC  =  IC+1 

CALL  ALL  t  VI*  YVP*  HXX,  HW«  HXY  ) 

IC  =  1 

DO  8  KJ*l»NROW 
DO  8  KI=KJ,NROW 

XWWtIC)  «  HOR(KJ)  *  HOR(KI)  -  SMI  *  (HYY(KJ)  *  HXX(KI) 

1  +  HXX (KJ)  *  HYY(KI)  -  2*  *  HXY(KI)  *  HXY(KJ)  > 

8  IC  =  IC+1 

DO  9  KJ  =  1*  NROW 

DO  9  KI  *  KJ*  NROW 

XWW(IC)  *  VER(KI)  *  VER(KJ) 

9  IC  *  IC+1 

CALL  ALL  (  VI*  YVM*  HXX*  HYY .  HXY  ) 

IC  *  1 

DO  10  KJ*1»NR0W 
DO  10  KI=KJ»NROW 

XVJIJ  *  WIJ  *  (XWi;(IC)  +  HOR(KI)  *  HCR(KJ)  - 

1  SMI  *  (  HYY(KJ)  *  HXXIKI)  +  HXX(KJ)  *  HYY(KI)  -  2.  *  HXY(KI)  * 

2  HXY(KJ)  1  ) 

AREAV(IC)  =  ARhAV(IC)+XWlJ 

10  IC  =  IC+1 

DO  11  KJ  =  1*  NROW 
DO  11  KI  *  KJ.  NROW 

XWIJ  =  WIJ  *  (  XWWtIC)  +  VER(KI )*VER(KJ) ) 

AREAV(IC)  =  AREAV( ICJ+XWI J 

11  IC  *  IC+1 
14  CONTINUE 

DO  16  K=1*N0 

AREA(K)  =  AREA(K)+WII*(DU*AREAU(K)+DV*AREAV(K))/2. 

AREAU(K)  =  0. 

16  AREAV(K)  *  0. 

20  CONTINUE 

DO  30  K=1»N0 
30  AREA(K)  =  • 5*AREA (K) 

RETURN 

END 


SUBROUTINE  ALL  (  X*  Y*  HXX*  HYY*  HXY  ) 

DOUBLE  PRECISION  XMA( 21  * 2 1 ) *XM5 ( 2 1 » 2 1 ) *XMC (2 1 . 2 1 ) *X 1 ( 48 ) * Y I ( 48 ) 
DOUELE  PRECISION  Wl (48 ) *HOR ( 21 ) *V£R (21 ) * AREA ( 462 ) *AREAU(462 ) 
DOUBLE  PRECISION  AREAV < 462 > *XWW(462 ) 

DIMENSION  XP ( 21 ) *  YP(21) 

DIMENSION  NXP( 21 )  ♦  NYP(21) 

DOUBLE  PRECISION  HXX ( 21 ) *HYY( 21 ) *HXY( 21 ) *B  121* 

DOUBLE  PRECISION  X,Y*F*FX*FY,FXX*FYY*FXY 

DOUBLE  PRECISION  DF*XIP* YJP *G*GX.GY*GXX *6YY»6XY*DG*P.AI *AJ 
COMMON  XMA,  XMB.  XMC*  XI*  YI*  Wl *  HOR*  VER.  AREA*  AREAU*  AREAV* 

1  XWW*P.B»ALPHA*8ETA*RATI0»NK»NR0W»XP.YP»AM1»BM1 

**#*#*•»  a#***#****##***#**####******##******#***###*##***#*####### 

CALL  VECTOR  (X*Y.F»FX»FY »FXX .FYY.FXY* ALPHA *BETA*P * AMI *BM1  ) 

**********##*»***#**«**#******###•»*#**+*#*****•»*+#+#**+#*#**#**♦* 
DF  =  FXX  +  FYY 


1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1880 

1890 

19u0 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 


2180 

2190 

2200 

2210 

2220 

2236 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 
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DO  20  KK  =  1*  XROK 

2336 

NXP(KK) 

S 

XP 1  KIO 

2340 

NYPtKK) 

= 

YP(KK) 

235C 

XIP 

s 

X  **NXP(KK) 

2360 

YJP 

m 

Y  **NYP(KK) 

2370 

A! 

s 

XP(KK) 

2380 

AJ 

2 

YPIKK) 

2390 

G 

= 

XIP  *  YJP 

246j 

GX 

= 

AI  *  G  /  X 

2410 

GY 

2 

AJ  *  G  /  Y 

2420 

GXX 

= 

(AI  -  l.i  *  GX  /  X 

2430 

GYY 

2 

(AJ  -  1.)  *  GY  /  Y 

2440 

GXY 

= 

AJ  *  GX  /  Y 

2450 

DG 

= 

GXX  +  GYY 

246ci 

HOR(KK) 

2 

G  *  DF  +  F  *  DG  +  2**(FX*GX  +  FY*GY  ) 

2470 

HOR(KK) 

= 

H0R(KK)*100C0ooCUu0. 

248U 

VER(KK) 

S 

F  *  G 

2496 

VER(KK) 

= 

VER ( KK )*10C000000004 

2500 

HXX(KK) 

2 

FXX*G  +  F*GXX  +  24*FX*GX 

2510 

HXX(KK) 

2 

HXX(KK) *10000000000« 

2520 

HYY IKK) 

2 

FYY*G  +  F*GYY  +  24*FY*GY 

2530 

HYY  ( KK ) 

2 

HY  Y ( KK ) *10000000000  4 

2540 

HXY(KK) 

2 

FXY*G  +  F*GXY  +  FX*GY  +  FY*GX 

2550 

HXY(KK)  = 

CONTINUE 

RETURN 

END 

HXY(KK) *100000000004 

2560 

257C 

2580 

259v, 

SUBROUTINE  VECTOR 1X»Y*F*FX»FY»FXX»FYY*FXY * ALPHA *u£TA»  P.AMl.BMl  )  2600 
DOUBLE  PRECISION  X.Y.F.FX.FY.FXX.FYY.FXY  2610 
DOUBLE  PRECISION  XA»PYD*FR1 »FR2»DX  *DY  *P  262C 
NALPH  =  I  FIX  I ALPHA)  2630 
NBETA  =  IFIX(BETA)  2640 
XA  =  X**NALPH  2*50 
PYB  =*  P  *  Y**NBETA  2660 
FR1  =  1.  -  XA  -  PYB  267u 
FR2  =  FR1  *  FR1  2680 
F  =  FR2  2690 
DX  =  -  ALPHA  *  XA  /  X  2700 
DY  =  -  BETA  *  PYB  /  Y  2710 
FX  =  2.  *  FR1  *  DX  2720 
FY  *  2.  *  FR1  *  DY  2730 
FXY  =  2.  *  DX  *  DY  2740 
FXX  =  2.  *  FR1  *  DX  *  AMI  /X  +  2.*  DX  *  DX  2750 
FYY  •=  2.  *  FR1  *  DY  *  BM1  /  Y  +  2.  *  DY  *  DY  2760 
RETURN  2770 
END  2780 


SUBROUTINE  SMTRXI  A*  C.  N*  E.  XU  )  2790 
TO  TRANSFORM  (C-W2A)X  =  0  INTO  BX=W2X  2800 
DOUBLE  PRECISION  A(21* 21 ) *C(21»21 ) *XL(21»21) .XUI21.21) .0(21*21)  2810 
DOUBLE  PRECISION  EI21.21)  2820 
CALL  SMTRXKA*  XL*  XU*  N  )  2830 
CALL  SMTRX2  <  XL.C*  D*  N  )  2840 
CALL  SMTRX3  t  XU*  D.  E*  N  '  2350 
RETURN  2860 
END  287o 
SUBROUTINE  SMTRXI (  A*  XL.  XU.  K  )  2880 
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C  TO  FIND  L  AND  L'»  TO  STORE  IN  XL  AND  XU  2890 

DOUBLE  PRECISION  A(21»21) »XL<21»21) *XU(21»21J  2900 

DOUBLE  PRECISION  S  2910 

DO  5  I  =  I»  N  2920 

DO  5  J  =  1.  N  2930 

XU  ( I.JJ  =  0.  2940 

'5  XL  (IiJI  =  0«  '  2950 

XUC1»1)  =  DSORT ( A( 1 <1 ) )  2960 

XL(ltl)  =  XU(1*1J  2970 

DO  15  IC  =  2*  M  2980 

XU(1»  IC)  =  AU.  IC)/  XU(l.l)  2990 

15  XL( IC* 1 )  =  XU(ltlC)  3000 

DC  100  I  =  2i  N  3ol0 

IP1  "  I  +  I  3^20 

IM1  =1-1  3030 

S  =  0«  3040 

DO  20  K  =  1»  I Ml  3050 

20  S  =  S  +  XU(K.I)  *  XU(K»I )  3060 

XU(I.I)  =  DSGRT (A( I » I )  -  S  )  3070 

XL  {  I » I  *  =  XUU.I)  3080 

IF  (  I  -  N  )  23 »  100»  ICO  3090 

23  DO  30  J  =  IP1»  N  3100 

S  =  0«  3110 

DO  25  K  =  1*  IM1  3120 

25  S  =  S  +  XU(K.I)  *  XUOUJ)  3130 

XU(I»J)  =  (A(I#J)  -  $)/XUII«l>  3140 

30  X^J.I)  =  XU(I.J)  3150 

100  CONTINUE  3160 

RETURN  317C 

END  3180 

SUBROUTINE  SMTRX2  <XL»  C.  D*  ft  )  3190 

C  TRANSFORM  TC  (L)-lC  AND  STORE  IN  D  3 2 " 

DOUBLE  PRECISION  XL ( 2 1. 21 )  *C( 21  ,  21 )  *D(  21» 21 )  321- 

DOUBLE  PRECISION  S  3220 

t /u  9  t  *  n  3  2  3  C 

3  Dud)  -  CiiSiJ  / 

UO  luu  i  =  2*  n  3250 

IM1  «  I  -  1  3260 

DO  100  J  =  1»  N  3270 

S  «=  0«  3280 

DO  10  K  ■  li  IM1  329o 

10  S  =  S  +  XL( I »K)  *  DfK.J)  3300 

100  D( I »J)  *  (C(I.J)  -  S  )  /  XL  < I  *  I )  3310 

RETURN  3320 

END  3330 

SUBROUTINE  SMTRX3  <XU»  D.  E»  N  )  334o 

C  TRANSFORM  TO  <L)-1C<L')-1  AND  STORE  IN  E  3350 

DOUBLE  PRECISION  XU ( 21 *2 1 ) »D< 21 f 21 ) t E ( 21 .21 )  3360 

DOUBLE  PRECISION  S  3370 

DO  5  I  *  1*  N  3380 

5  E  ( 1*1)  =  D(  I  *  1 )  /  XU(J.l)  3390 

DO  100  J  =  2*  N  3400 

JM1  *  J  -  1  3410 

DO  ICO  I  =  1 »  N  3420 

S  =0.  3430 

DO  10  K  =  1»  JM1  3440 

10  S  =  S  +  E 1 1 » K >  *  XU(KtJ)  3450 

100  EUiJ)  ■  ( D(  I  *J )  -  S  )  /  XU(J«J)  3460 
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RETURN 

3470 

END 

3480 

SUBROUTINE  EIGEN{A,NRANK*NROOT*ANSWER*VECTOR.LIMIT»CONVER. 

TELL,  3490 

l  NUMCYC  ) 

3500 

c 

THIS  SUBROUTINE  FINDS  THE  EIGENVALUES  AND  EIGENVECTORS  bY  AN 

ITERATIVE  S  3510 

c 

using  the  method  of  reduction. 

3520 

DOUBLE  PRECISION  A(21 »21) » ANSWER! 21 ) * VECTOR! 21 *21)*Z(21)*Y(21»21)  3530 

DOUBLE  PRECISION  GREAT *TRY»R 

.DIFF  *CONVER 

3540 

DIFF  *  0. 

3550 

DO  24  1=1 »NROOT 

3560 

WRITE  (6»  200)  I 

3570 

200 

FORMAT  (2H  *  13  ) 

3580 

"  '  j=i»nrank 

3590 

1 

Y(I»J)=1. 

3600 

NUMCYC=0 

3610 

2 

NUMCYC=MUMCYC+1 

3620 

WRITE  (6*  300)  NUMCYC*  DIFF 

3630 

300 

FORMAT  (  4X.  14*  5X*  D25.16 

*  3X ) 

3640 

I F ( NUMCYC-L IM I T ) 3  *  3  •  25 

365C 

3 

DO  4  J=I*NRANK 

3660 

Z(J)=0. 

3670 

DO  4  K=I*NRANK 

3680 

4 

Z<J)=Z(J)+A(J.K)*Y( I.K) 

3690 

GREAT  =  DABStZf 1) ) 

3700 

INDEX=I 

3710 

IFI I-NRANK) 5*8*8 

3720 

5 

K=I+1 

3730 

DO  7  J=K*NRANK 

3740 

TRY  =  DABS(Z(J)) 

3750 

IF ( GREAT— TRY ) 6 *7  *7 

3760 

6 

GREAT=TRY 

3770 

INDEX=J 

3780 

7 

CONTINUE 

3790 

8 

DIFF=0. 

381-0 

GREAT=Z( INDEX) 

3810 

DO  9  J=I*MRANK 

3820 

Z ( J ) =Z ( J ) /GREAT 

3830 

9 

DIFF  =  DIFF  +  DA3S(Z(J) 

-  Y(I.JI) 

3840 

DO  10  J=I .NRAMK 

3850 

10 

Y(I.J)=Z(J) 

3860 

IF( DIFF-CONVER  .'11*11*2 

3870 

11 

ANSWER ( I ) =GREAT 

3880 

GREAT=Z( * ) 

3390 

DO  12  J=I»NRANK 

3900 

Zt J) =Zt J) /GREAT 

3910 

12 

Y( I»J)«Z( J) 

3920 

IF( I-NROOT) 13*15*15 

3930 

13 

L=I  +  1 

3940 

DO  14  J=L*NRANK. 

3950 

DO  14  K=L  *NRANK 

3960 

14 

A(J,KlsAIJ*KI-Z(JI*AII»K) 

3970 

15 

IFt 1-1)20*20*16 

3980 

16 

DO  19  J=2 » I 

3990 

L= I-J+l 

4000 

M=L+1 

4010 

R=0  • 

4020 

DO  17  K=M*NRANK 

4030 

17 

R=R+A t  L.K ) *Z ( K ) 

4040 
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R=R/ (ANSWER! X ) -ANSWER  (U )  4050 

Z(D=1*  4060 

00  18  K*M»NRANK  40?0 

18  Z(K ) =Y(L»K)+Z (K ) /R  4080 

19  CONTINUE  4050 

20  GREAT  *  DABSIZim  41C0 

IMDEX=1  4110 

DO  22  J«2»NRANK  4120 

TRY  *  DABS(Z( J) }  4130 

IF (GREAT-TRY  121*22 *22  4140 

21  GREAT=TRY  4150 

INDEX=U  4160 

22  CONTINUE  4170 

GREAT*Z( INDEX)  4180 

DO  23  J«1*NRANK  4190 

23  VECTORl I »JI*Z(JJ/GREAT  4200 

24  CONTINUE  4210 

TELL=1.  4220 

RETURN  4230 

25  TELLa“l*  4240 

RETURN  4250 

END  4260 


SUBROUTINE  TRAVEC  (XU*  X*  PHI*  NRCri  )  4270 

DOUBLE  PRECISION  XU( 21 *21 J *X ( 21 J »PHI t 21 J  4280 

DOUBLE  PRECISION  SUM  4290 

N  =  NROW  4300 

NM1  =  N  -  1  4310 

PHI ( N )  =  X ( N J  /  XU ( M * N )  4320 

DO  100  I  =  1#  NM1  4330 

J  =  N  -  I  4340 

SUN  =  0.  4350 

DO  80  K  =  Ji  NM1  4360 

KP1  s  K  +  1  4370 

80  SUM  =  SUM  +  XUIJ.  KPl)  *  PHKK.P1)  4380 

100  PHI ( J )  =  (X(J!  -  SUM  )/XU(J.J)  4390 

RETURN  4400 

END  4410 

SUBROUTINE  PLNODE  (HP)  4420 

DOUBLE  PRECISION  XMA( 21*21 ) »XMB (21 »21 > *XMC< 21 *21 ) *XI (48 ) »YI (48 )  4430 

DOUBLE  PRECISION  Wl (48) *HOR (21 J *VtR(21 > *ARLA( 462 ) *AREAU( 462 J  4440 

DOUBLE  PRECISION  AREAV<46Z) *XV:W(462)  445C 

DIMENSION  XP ( 21 ) *  VP(2l)  4460 

DOUBLE  PRECISION  B( 21 5 *VXP( 21 ) * VYP< 21 ) *R( 50 )  4470 

DOUBLE  PRECISION  V*XNP*ERRGR»STEP ,P  4480 

COMMON  XMA»  XMo,  XMC*  Xi*  YIt  V.'I*  HOR i  VlR*  AREA*  AREAu*  AREAV *  4490 

1  XWW*P*B*ALPHA.b£TA,RATlO,NK*NROu»XP.YP*AMl*BMl*  4500 

2  SWITCH.  VXP.  VYP  4510 

ERROR  *  0.0001  4520 

STEP  *  0*05  4530 

SWITCH  =  1.  4540 

XNP  *  NP  4550 

DO  50  1  «  1*  NP  4560 

AI  =1-1  4570 

V  «  AI  /  XNP  4580 

IF  (  V  I  20*  10*  20  4590 

10  DO  18  IX  =  1*  NROW  4600 
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IF  (  XP(IX)  I  16*  14*  16 

4610 

14 

VXP(IX)  *  1* 

4620 

GO  TO  18 

4630 

16 

VXPIIX)  =  0. 

■  4640 

18 

CONTINUE 

4650 

GO  TO  28 

4660 

20 

DO  26  IX  =  1*  NROW 

4670 

IF  (  XP(IX)  )  24*  22*  24 

4680 

22 

VXP(IX)  =  1. 

4690 

GO  TO  26 

4700 

24 

VXP(IX)  =  V  **  XP(IX) 

4710 

26 

CONTINUE 

4720 

28 

CALL  REGSUN  (  0.*  1.*  STEP 

.  R.  NR*  ERROR  ) 

4730 

IF  (  NR  )  50*  50*  44 

4740 

44 

WRITE  (6*  1400 )  V*  (R( J) • 

J  =  1.  NR  ) 

4750 

1400 

FORMAT  <  IX*  3HX  =  ,  F 6.3* 

2X,  3NY  =»  lVt-6.3 

/(13X*lVf6.3J ) 

4760 

50 

CONTINUE 

4770 

SWITCH  *  3* 

478C 

DO  80  I  =  1*  NP 

4790 

AI  =1-1 

4800 

V  =  AI  /  XNP 

4810 

IF  (  V  )  60.  52*  60 

4820 

52 

DO  58  IY  =  1.  NROW 

483u 

IF  «  YP(IY)  5  56*  54*  56 

4840 

54 

VYP(IY)  =  1. 

4850 

GO  TO  58 

4860 

56 

VYP(IY)  =  0. 

4870 

58 

CONTINUE 

4880 

GO  TO  68 

4890 

60 

DO  66  IY  =1.  NROW 

4900 

IF  (  YP(IY)  )  64.  62*  64 

4910 

62 

VYP(IY)  =  1. 

4920 

GO  TO  66 

4930 

64 

VYP(IY)  =  V  **  YP(IY) 

4940 

66 

CONTINUE 

4950 

68 

CALL  REGSUN  (  0.*  1.,  STEP.  R.  NR,  ERROR  ) 

4960 

IF  {  NR  )  80.  80.  74 

4970 

74 

WRITE  (6*  1600)  V.  (R( J 1  * 

J  =  1.  NR  ) 

4980 

1600 

FORMAT  (  IX.  3HY  =,  F6.3* 

2X »  3HX  =»  19F6.3 

/(13X.19F6.3) ) 

4990 

80 

CONTINUE 

5000 

RETURN 

5010 

END 

5 -'20 

A 


'« 

■i 


SUBROUTINE  REGSUN  (  A,  b»  H,  R,  N»  ERROR  )  5030 


C 

TO  FIND  ALL  ROOTS  OF  EIGENVECTOR 

5*40 

DOUBLE  PRECISION  R t  50 ) .ERROR* XL »XR  »YL.YR*Xl*H*YI 

5050 

N  =0 

5o6J 

XL  ■  A 

5070 

4 

YL  =  FUNCT I  XL ) 

5o80 

IF  (  DABS(YL)  -  0.1D-10  ) 

10*  lo*  2o 

5090 

10 

N  =  N  +  1 

5100 

R  ( N )  =  XL 

5110 

XL  =  XL  +  H 

5120 

IF  (  XL  “  B  )  4*  4*  16 

5130 

16 

RETURN 

5140 

20 

XR  =  XL  +  H 

5150 

IF  (  XR  -  B  )  22*  22.  16 

5160 

- 

22 

YR  =  FUNCT (XR) 

5170 

IF  (  DABS(YR)  -  0.1D-10  ) 

30*  30*  24 

5180 
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24 

IF 

< 

YR*YL  )  40.  30.  60 

5190 

30 

N 

«  N  +  1 

5200 

R(N) 

*  XR 

5210 

XL 

*  XR  +  H 

5220 

IF 

( 

XL  -  B  )  4.  4.  16 

5230 

40 

XI 

*  (  XR  +  XL  >  /  2. 

5240 

IF 

( 

XI  -  XL  -  ERROR  )  46.  46.  48 

5250 

46 

N 

«  N  +  1 

5260 

R(N) 

*  XI 

5270 

XL 

«  XI  +  H 

5280 

GO 

TO  4 

5290 

48 

YI 

«  FUNCT (XI ) 

5300 

IF 

( 

DABS(YI)  -  0.1D-10  )  46.  46.  50 

5310 

50 

IF 

( 

YL*YI  )  52.  46.  54 

5320 

52 

XR 

*  XI 

5330 

GO 

TO  40 

5340 

54 

XL 

*  XI 

5350 

GO 

TO  40 

5360 

60 

XL 

=  XR 

5370 

YL 

YR 

5380 

GO  TO  20 
END 


539o 

54l.j 


FUNCTION  FUNCT(O)  5410 

OOUBLE  PRECISION  XMA< 21 >21 ) »XMB< 21 .21 ) *XMC( 21 *21 ) .XI (48 ) >YI 148)  5420 

OOUBLE  PRECISION  Wl ( 48 ) .HOR ( 21 ) . VER ( 21 ) *AREA( 462 ) »AREAU 1 462 )  5430 

DOUBLE  PRECISION  AREAVt 462  )  »Xv.U(462  )  5440 

DIMENSION  XP (21 ) *  YP(21)  5450 

DIMENSION  NXP ( 21 ) *NYP ( 21 )  5460 

DOUBLE  PRECISION  3 ( 21 I « VXP l 21 ) < VYP ( 21 )  5470 

DOUBLE  PRECISION  Q»  SUM. QYP »QXP .FUNCT .P  5460 

COMMON  XMA.  XMU.  XfX.  Xl»  Yl»  Ul*  HOR*  VER.  AREA.  aREaU.  AREAV.  5490 

1  XWW.P.G. ALPHA. BETA. KATI0.NN..i.R0'..»XP>YP»Am1.l>..1.  5500 

2  SWITCH.  VXP.  VYP  5510 

DO  500  I  =  l.NROW 

NYP(I)  =  YP(I)  5520 

NXP(I)  =  XP(I)  5530 

500  CONTINUE 

IF  <  SWITCH  -  2.  )  2.  20.  20  5540 

2  SUM  -  0#  5550 

DO  10  1=1.  NROW  5560 

IF  (  YP ( I )  )  4.  3.  4  5570 

3  OYP  =  1.  5580 

GO  TO  10  5590 

4  IF  (0)  6*  5.  6  5600 

5  QYP  =  0.  5610 

GO  TO  10  5620 

6  QYP  ■  0  **NYP(l)  5630 

10  SUM  *  SUM  +  B ( I )  *  VXP(I)  *  QYP  5640 

FUNCT  *  SUM  5650 

RETURN  5660 

20  SUM  =  0.  5670 

DO  30  I  =1*  NROW  5680 

IF  (  XP( I J  )  24.  23*  24  5690 

23  QXP  =  1.  5700 

GO  TO  30  5710 

24  IF  (0!  26.  25.  26  5720 

25  OXP  =  0.  5730 

GO  TO  30  5740 


26  QXP  =  Q  **NXP(I)  5750 

30  SUM  =  SUM  +  Bill  *  QXP  *  VYP(I)  5760 

FUNCT  ■  SUM  5770 

RETURN  5780 

END  5790 
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The  principle  used  to  solve  for  the  natural  frequencies  is  the  Rayleigh-Ritz  method. 
The  plate  geometry'  is  defined  by 

F(X,  Y,  P,a,  (i)  =  (1  -Xa  -  PY&)2 


where  X  =  — , 
a 


P  =  R-P, 

<*  =  / 3  =  10, 

a  is  the  dimension  ir.  ^-direction,  and 

b  is  the  dimension  in  y-direction  for  a  rectangle  with  clamped  boundaries. 

In  the  computer  program,  the  Rayleigh-Ritz  procedure  uses  a  21-term  polynomial  in  A' 
and  Y  to  express  the  displacement  W  (Equation  (G9)).  The  integrals  of  the  Rayleigh-Ritz 
equations  are  then  solved  by  a  64-order  Gaussian  quadrature  technique.  Finally,  the  eigen¬ 
values  of  Equation  (Gl3)  are  solved  by  an  iterative  method  of  reduction. 

The  computer  program  solves  for  one  set  of  frequencies  at  a  time.  Four  sets  of  poly¬ 
nomials  completely  define  the  plate:  even-even,  odd-odd,  even-odd,  odd-even.  Manual  plotting 
of  the  nodal  points  for  the  first  quadrant  yields  the  modes  shapes  from  which  the  modal  num¬ 
bers  may  be  assigned  to  the  frequencies. 

The  eigenvalues  resulting  from  the  computer  program  are  actually  the  dimensionless 
frequencies  (note:  a  =/  2nf  in  this  program) 


o  =  a2  \Jp2 
m,n  y 


/Yh  ' 
n\gDj 


where  the  pm  represent  the  natural  frequencies.  Thus,  the  program  eigenvalues  must  be 

Eh 3 

modified  manually  to  yield  frequencies  in  hertz.  Letting  pmn  =  2 'vfmn  and  D  =  - 

12(1  -<72) 

Equation  (II)  becomes 


CO  h 
m,n 


2nd2  V  12y(l-a2) 


In  addition  to  the  eigenvalues,  the  program  computes  the  points  for  the  nodal  lines  to 
be  plotted  to  give  the  mode  shapes. 

A  sample  problem  for  eight  modes  with  32-order  Gaussian  quadrature  required  30  min¬ 
utes  on  the  IBM  7090. 
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input  Description 


The  input  data  are  in  dimensionless  form. 


Their  description  is  as  follows. 


Program  Symbol 

Theory 

Symbol 

Description 

Format 

NK  (card  1) 

N 

The  value  —  where  N  is  the  order  of 

2 

Gaussian  quadrature 

110 

Beginning  on  card  2,  start  the  XI  array  and  end  with  WI  array;  last  card  of  this  set  is 

/  NK\ 

card  i1  -  TV 

XI 

Gaussian  arguments;  NK  elements;  4  to  a 
card 

4D20.10 

WI 

Gaussian  weights;  NK  elements;  4  to  a  card 

4D20.10 

Next  8  elements  are  on  the  ^2  +  card 

ALPHA 

a 

Exponent  of  plate  geometry  equation: 

ALPHA  -  10  for  rectangle 

F5.2 

BETA 

P 

Exponent  of  plate  geometry  equation: 

BETA  =  10  for  rectangle 

F5.2 

RATIO 

R 

Aspect  ratio  i/a ,  where  b  is  dimension  in  in¬ 
direction  and  a  is  dimension  in  ^-direction 

F5.2 

MODE 

The  number  of  sets  of  modes  desired. 

If  MODE  = 

1  X,  Y  are  even  powered:  odd-odd  modes 

2  X,  Y  are  odd  powered:  even-even  modes 

3  X  even,  Y  odd:  odd-even  modes 

4  X  odd,  Y  even:  even-odd  modes 

15 

NOIT 

Number  of  eigenvalues  desired 

15 

NP 

Number  of  nodal  points  desired: 

NP=0  means  no  points 

NP  =  20  means  20  points 
for  nodal  line  plot 

15 

LIMIT 

Number  of  iterations  in  eigenvalue  solution; 
suggested  limit  is  800 

15 

CONV 

Convergence  criterion:  suggested  value 
0.00001 

F10.7 

NROW  ,(card  3  + 

Number  of  polynomials  in  X  and  Y 

no 

XP(I),  YP(I) 

/  NROW *2 

1  next -  cards 

V  16 

for  MODE  number  of  \ 
times  / 

Powers  of  terms  of  X  •  Y  polynomial; 
note  that  there  must  be  as  many  sets  as  the 
value  of  MODE  indicates  but  that  the  pre- 
gram  solves  for  only  one  set  at  a  time 

16F5.2 
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Sample  input  data  corresponding  to  the  above  description  are  shown  below: 


16 

0.*  9972638618494815600 .98561 151 15452683303*  964762255  5  8750643l>0 
0*8963211557660521200. 8493676137325699700. 794483795967942401:0 
0.6630442669302152000. 5377157572407623200.5C68999089322293900 
0*33 1868602282 1276400.2392873622521 370700 .144471 96 15827964900 
0*0070186100094700900.0162743947309056700.0253920653092620500 
0.0428356980222266800.0509980392623761700.0586840934785355400 
0.0723457941068485000.0781938957870703000.0833119242269467500 
0»091 1738736957638800. 0938443990308045600.0956387200792746560 


.934906075937739680 
.732182118740289630 
.421351276130635340 
.046307665687738310 
.034273862913021430 
.065822222776561 640 
.087652093034403310 
.096540088314727300 


lo.e 

10.0 

2 

0*0 

1.167 

4 

4 

8  2C-  80 

0  O.OCCCOOl 

0*0 

*2.0 

0.0 

4.0 

0.0 

6.0 

0.0 

8.0 

0.0 

10. 0 

0*0 

0.3 

2.0 

2.0 

2.0' 

19 

4.0 

2.C 

6.0 

2.0 

8.0 

2.0 

0.0 

4.0 

2.0 

4.0 

4.0 

4.0 

6.0 

4.0 

0.0 

6.0. 

/§ 

2.0 

6.0 

4.0 

6.0- 

0.0 

8.0 

2.0 

8.0 

e.c 

10.0 

1. 

1. 

3. 

1. 

5. 

1. 

7. 

1. 

9. 

1. 

11. 

1. 

1. 

3. 

3. 

3. 

5. 

3. 

7. 

3. 

9. 

5. 

1. 

5. 

3. 

5. 

5. 

5. 

7. 

5. 

1. 

7. 

3. 

7. 

5. 

7. 

1. 

9. 

3. 

9. 

1. 

11. 

o. 

1. 

2. 

1. 

4. 

1. 

6* 

1. 

8. 

1. 

10. 

1. 

0. 

3. 

2. 

3. 

4. 

3. 

6. 

3. 

8  . 

3. 

0. 

5. 

2. 

5. 

4. 

5. 

6  . 

5. 

C. 

7. 

2. 

7. 

4. 

7. 

0. 

9. 

2. 

9. 

o. 

11. 

1.0 

0*0 

1«0 

2.0 

1*0 

4.0 

1.0 

6.0 

1.0 

8.0 

1.0 

10.0 

3.0 

0.0 

3.0 

2.0 

3.0 

4.0 

3.0 

6.0 

3.0 

6.0 

5.0 

0.0 

5.0 

2.0 

5.0 

4.0 

5.0 

6.0 

7.0 

0.0 

7.0 

2.0 

7-0 

4.0 

9.0 

0.0 

9.0 

2.0 

11*0 

0*0 

a 

aa 

'O  l 


o  > 


ro  , 


S 


Output  Description 

The  program  yields  the  eigenvalues  and  eigenvectors,  with  nodal  points  for  the  first 
quadrant  and  many  intermediate  results.  Unless  the  user  is  particularly  interested  in  a  pro¬ 
gramming  analysis,  he  will  use  the  first  page  of  output  and  then  skip  to  the  eigenvalue 
section. 

On  the  first  page  are  some  of  the  input  data,  such  as  o  ,  /3,  RATIO,  MODE,  which  are 
labelled  accordingly.  The  index  I  is  printed  to  indicate  the  step  of  Gaussian  quadrature.  An 
underflow  message  from  the  system  may  occur;  the  program  corrects  for  small  numbers  in  the 
underflow  in  subroutine  ALL. 

The  next  several  pages  have  five  elements  to  a  row  and  are  the  following  matrices: 

1.  (7-matrix  of  Equation  (Gl4a) 

2.  A-matrix  of  Equation  (Gl4b) 

3.  B-matrix  of  Equation  (G  15b) 
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A 


A 


j 


'i 

A 


i 


i 


4 


>1 
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The  output  then  indicates  which  eigenvalue  is  being  solved  for  and  the  number  of 
iterations  needed.  The  variable  TELL  indicates  convergence:  TELL  =  1  means  convergence 
but  TELL  =  -  1  means  no  convergence.  The  convergence  limit  and  the  number  of  times  the 
iterations  are  performed  are  also  printed.  The  eigenvalues  are  printed  in  ascending  order, 
followed  by  the  eigenvectors.  The  results  of  the  orthogonality  check  are  shown. 

Finally  for  a  given  eigenvalue  the  nodal  points  for  the  first  quadrant  are  printed  out. 
Figure  28  shows,  by  way  of  a  particular  example,  how  the  mode  shapes  and  corresponding 
frequencies  are  matched.  The  eigenvalues  (called  EIGENVALUE  in  the  output  data)  obtained 
directly  as  output  from  the  computer  program  are  multiplied  by  the  frequency  factor  for 
SUNFRE  given  in  Appendix  I.  This  process  yields  the  natural  frequencies  which  are  tabu¬ 
lated  in  Table  1. 

Thus  for  a  particular  eigenvalue  (e.g.,  EIGENVALUE  =  337.0694),  a  corresponding 
natural  frequency  can  be  computed  (/  =  2179.078  for  this  case).  The  corresponding  mode 
number  can  be  determined  by  plotting  wave  shape  data  available  from  the  computer  program. 
These  data  are  plotted  in  the  first  quadrant  (Figure  28a)  and  then  projected  into  all  four 
quadrants  (Figure  28b).  From  the  latter  figure,  the  mode  number  is  evidently  (m,n)  =  (5, 2). 

YNGFRE  (see  Table  12  and  Figure  29) 

T-o  steps  are  needed  to  find  the  natural  frequencies  of  vibrations  by  the  Young 
method.  The  first,  YOUNG,  provides  preliminary  data.  The  second,  YEIGN,  computes  the 
eigenvalues  and  converts  them  to  the  natural  frequencies.  Since  the  results  of  YOUNG  could 
be  used  as  input  for  other  eigenvalue  programs,  YOUNG  was  made  more  general  than  YEIGN. 

YOUNG 

YOUNG  is  a  computer  program  which  calculates  the  members  of  the  (7-array  of  the 
eigensystem.  Equation  (Bll): 


P  *  .. 

X  2 

_  i  _ _ «  V  mn  mn'  mn 

m  *  l  n  —  l 


5  =  1  for  m  =  i  and  n  =  k 

mn 


S  =  0  for  m  4  i  or  n  4  k 

mn 

For  the  computer  program,  i  =  1,  p\  k  =  1,  g;  and  p,  q  -  10. 

The  program  YOUNG  uses  its  subroutine  YINTGR  to  compute  numerical  results  of 
Young’s  closed  form  solutions  of  the  Rayleigh-Ritz  integrals  of  a  clamped  beam.  Next 
YINTGR  constructs  the  arrays  necessary  for  the  computation  of  the  (7-matrix: 
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TABLE  12 

Program  Listing  for  YNGFRE  Computer  Program 

Table  12a  -  YOUNG 

3IBFTC  YOUN6 

DIMENSION  El 10*10) »Fll0*10)*Hll0*l0>»Kll0*10>*EP$tl0)*CI 10*10) 

REAL  K 

RFADl5*l)  M«N 
1  FORMAT! 215) 

READ(5*11)  A*B 
11  FORMAT I 2Fl?#6 ) 

PI  *  5*14159 

CALL  YINTGR|M*N*EPS*E) 

WRITE<6»5)  I  EPS  1 1 ) » I  >•  1*M) 

5  FORMAT 1 6X  *5E16*8 ) 

SI6MA  -.33 
A3«A**3 

WRITE(6*310)  A»C 
310  FORMAT 1 5X  *2F12  *6 ) 

WR1TE(6*320)  M»N 
320  FORMAT I 5X *2 15) 

NY1  «  N/2 
NY2  •  N/2  +  1 
DO  4  I»1*M 
DO  3  J»1»N 
HII*J>  ■  EII.J) 

K(1*J)  •  E!  1*  J) 

El  I »J)  «  “Ell* J) 

F( I *J)  «  El  1* J) 

3  CONTINUE 

4  CONTINUE 
XOUNT-O 

DO  400  I«1«M 
DO  300  J*1*N 
DO  200  MX«1*M 
DO  100  NY»1*N 
IFIMX.NE.I)  60  TO  9 

IFINY.EO.J)  60  TO  6 

8  C 1 MX  *NY )»SI6MA*A /B* I E I MX  » I4*F I J  *NY ) +F 1 1 »MX )*F I NY  * J ) ) 

1  +2.*I1*“SIGMA)*A/B#HII*MX)#X(J»NY1 
60  TO  T 
6  CONTINUE 

C«MX*NY)«fl/A*EPSII)»*4+A3/B3»EPSU)*»4+2*»S!6MA*A/B»EII»I)*F|J*J) 

l+2**ll«“SIGMA)*A/8»H<t*!)*X<J*J> 

CONTINUE 
100  CONTINUE 

X0UNT*K0UNT+1 

COMMENT  XOUNT  WAS  USED  FOR  ENDPUNCHING****  NOW  IT  USED  ONLY 
THE  CASF  N  IS  A  MULTIPLE  OF  2*  ******** 

IF«M-IM/3*5>>  200*230*240 

240  IF|M“(M/3*3) )  200*210*220 

220  lFIM-IM/2*2> >  200*222*230 

210  WRITE! 6*20)  I  CIMX»NY>»  NY  «1*N) 

WRITEI 8*20)  I  CIMX»NY)»  NY  «1*N) 

GO  TO  200 

222  WRITEI6»22)  I  C(MX*NY)»  NY  »1*NY1>*  XOUNT 
WRITEI 6*22)  «  CIMX»NY)»  NY  «NY2*N>  »  XOUNT 
WRITE! 8*22)  I  ClMX*NY)*  NY  *l*NYl ) *  XOUNT 
WRITEI 8*22)  I  CIMX»NY) »  NY  «NY2»N)»  XOUNT 
GO  TO  200 

250  WRITEI6.24)  I  CIMX»NY)*  NY  »1»NY1) 

WR!7EI6*24)  I  CIMX»NY)»  NY  «NY2»N> 

WRITEI8*24)  1  CIMX»NY) »  NY  «1»NY1> 
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TABLE  12*  (Continued) 


M*ITE(8a24)  (  C(MX.NY) •  NY  *riY2«N) 

200  CONTINUE 
300  CONTINUE 
400  CONTINUE 
ENDFILE  • 

20  F0NMATC3E16a8) 

22  F0(MATf4E16a8»12XaI4) 

24  F0m*ATfSE16a8) 

230  STOP 

END 

*If»FTC  YINTE» 

StftWOUTINE  YlNTGRtH.N'FPS.Al 

DIMENSION  ALPdOJ.EPSdOJ.AdOalOI 

PI  ar  3 a  14139 

ALPC11  -  0.95250226 

ALP!2)a>  la000TT731 

ALP13J  *  0.99996643 

ALP  141  .at  la0000014S 

ALPI51  at  0.99999994 

ALPIftI  «  laO 

EPSdl  at  4a73 004080 

EPSI2I  •  Ta 89320460 

EPSI3)  *  10a99 56078 ft 

EPS 14)  *  14.13716930 

EPSI51  *  17.27879960 

EPSI6)  *  20.4239522 

0O  10  J  •  7aH 

ALPIJt  «  1.0 

AJ  -  J 

10  EPS! J)  *  f (2aO*AJ  ♦  laC)*PI)/2.0 
00  25  K  *  1*M 
00  35  L  *  1*N 
XL  «  X.  *  L 
1FU.NF.LJ  CO  TO  40 

AOC.L1  •  ALPfK)*EPS(K)*<ALPfKJ#EPSU)-2.ft) 

60  TO  35 

40  ACICaL)  <.  CC4.0*EPS(U**2*FPStK)**2)*IALFlL)*FPSCL) 

1  -ALPIK)*EPSIKI> 

2  *C1.4i-l.)**nCL  )1)  /  (EPSIKHK4  -  EPS<L)P*4) 

35  CONTINUE 

25  CONTINUE 

WRITE(6«50)  ( (A (KM. KM) »KM  ■  l.MJ.KN  »1»N) 

50  F0ft«AT(2Xf5E16aS> 

RETURN 

ENO 
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Table  12b  -  YEIGN 


PROGRAM  YEIGNt INPUT tOUTPUT * TAPE5= INPUT *TAPE6=OUTPUT ) 

DIMENSION  B( 64*64) »RTR(64) >RTI (64) *U( 64*64) *1X1(64) *1X2(64) » 

A  X1(64)*X2(64)*X3(64) »X4(64i »XX1(64»64) >XX2(64»64) » 

B  XX3(64*64) »EVLRAD(64) *EVTRAD(64f 64) *X5(64> »X6(64) » 

C  X7C64)  *X8(64)  *X9(64)  *X10t64)  *XIKb~: ) 

DOUBLE  PRECISION  BDP(64»64) »RTR:mP(64) .RTIIKP!©'*) *XIMP(64.64) » 

A  DPX1I64) *DPX2($4« 

COMMENT  AS  OF  11/20/70  LIM  MUST  BE  A  MULTIPLE  OF  3*4*  OR  5 
READ(5»110)  LlMtLUP 
110  FORMAT (21 10) 

N  «  LIM  **  2 
READ (5*115)  CONST 
115  FORMAT CE16.8) 

DCONS  =  DBLE< CONST) 

WRITE(6*3) 

3  FORMAT (1H1) 

WRITE! 6*1 )LlM»N 

1  FORMAT (21 10) 

IF(M0D(LIM*3)*EQ*0)  GO  TO  410 
IF ( MOD { LIM»5 i *EQ*0 )  GO  TO  420 
READ (5 *91) ( (B( IA*JA) *JA=1*N) *IA=1*N) 

91  FORMAT (4E16* 8) 

GO  TO  99 

410  READ(5*94)  ( (B( IA*JA) »JA=1»N) *IA=1»N) 

94  FORMAT (3E16«8) 

GO  TO  99 

420  READ (5*430! ( (B( IA*JA) *JA=1 *N) * IA*1 *N) 

430  FORMAT (5E16. 8) 

99  WRITE(6*4)((B( IA»JA)»JA=1»N) *IA=1»N) 

DO  10  1=1 »N 
DO  10  J«1»N 

10  BDP(I*J)=B(I»J) 

4  FORMAT! 1X.6E18.8) 

CALL  VARAH1(B*N»RTR»RTI»U»64»IX1*IX2*X1*X2»X3*X4.XX1*XX2»XX3) 
WRITE(6*3) 

WRITE! 6*5) ( I»RTR( I l »RTI ( I ) » I*1*N) 

5  FORMAT! I5.2E17.8) 

WRITE(6.3) 

DO  9  J=1»N 

WRITE(6*6).!»(U(I»J)»I=l»N) 

6  FORMAT (// 15/ (6E20.8)) 

9  CONTINUE 

DO  11  K  *  1»LUP 

CALL  VARAH2(8DP.N.2.0«*(-95)»RTR»RTI»U*RTRIMP»RTIIMP»EVLRAD#XIMP» 

1  EVTRAD..TRUE**  64* IX1»X1»X2»X3»X4»X5»X6»X7»X8»X9» 

2  X10»X11*DPX1»DPX2»XX1»XX2*XX3) » 

3  RETURNS (97) 

DO  12  I*1»N 
RTR( I )=RTRIMP( I ) 

RTI ( I )*RTIIMP{ I ) 

DO  12  J=1»N 

12  U( I »J)=XIMP(I»J) 

11  CONTINUE 
WRITE( 6*3) 

92  DO  14  I *1 *N 

IF  (RTRIMP(I).GE.l.O  D-12)  GO  TO  13 
DPX1(I)«“1.0  DO 
GO  TO  14 

13  DPXKI)  ■=  DCONS  *  DSQRT  ( RTR  IMP  ( I ) ) 

14  CONTINUE 
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TABLE  12b  (Continued) 


WRITE!6»250) 

250  FORMAT I 1H1»5X»*THE  FOLLOWING  IS  INTENDED  AS  A  GUIDE  IN  INTERPRETIN 
1G  THE  OUTPUT.*/  5X»*  THE  SUBSCRIPT  PRINTED  WITH  THE  EIGENVALUES  AN 
2D  FREQUENCIES  ON  THE  LAST  PAGE*/5X»*IS  THE  SUBSCRIPT  OF  ABS  LAMBDA 
3(  )  IN  THE  MAIN  SECTION  OF  OUTPUT —  EACH  EIGENVALUE  IS  PRINTED**/ 
45X » *FOLLOWED  IMMEDIATELY  BY  ITS  EIGENVECTOR. THE  SECOND  SUBSCRIPT 
5  OF  THE  EIGENVECTOR  COMPONENTS  AGREE*/5X»*WITH  THE  SUBSCRIPT  OF 
6 LAMBDA*) 

WRITE!6*240) 

240  FORMAT C5X>*HHcN  READING  EIGENVECTORS* LOOK  FOR  THAT  COMPONENT*/ 
l*WHOSE  VALUE  =1.0  .THE  FIRST  SUBSCRIPT  OF  THIS  COMPONENT*/ 

2  5X»*INDICATES  THE  MODE  NUMBER  OF  THE  FREQUENCY.*/ 

3  5X.*INTERPRETATI0N  SCHEME  BELOW  WITH  M»N  8EING  THE  MODE  NUMBER*/ 

4  6X.*JA**12X**M**7X»*N*) 

KOUNT  =  1 

DO  210  KM  =  1»LIM 
DO  202  KN  =  l.LIM 
WRITE! 6*310)  KOUNT »KM»KN 
31u  FORMAT (5X»I4*10X*I4*5X*I4) 

KOUNT  *  KOUNT  +  1 
202  CONTINUE 
210  CONTINUE 

WRITE(6*260> 

260  FORMAT (5X**THUS  BY  LOOKING  AT  THE  EIGENVECTOR  OF  EACH  LAMBDA*/5X* 
1*USER  MAY  ASSIGN  MODAL  NUMBERS  TO  THE  FREQUENCIES  BELOW*) 
WRITE(6.120) 

120  FORMAT (6X»  *£IGENVALUES  AND  CORRESPONDING  FREQUENCIES  *  ) 
WRITE(6»15)  ( I »RTRIMP ( I ) *DPX1 1 1 )  » I  *  l.N) 

15  FORMAT! I6.D25. 16*5X»D25. 16) 

STOP 

97  WRITE! 6*98) 

98  FORMAT !5X»*  PROGRAM  ABORTS  UNNATURALLY  *  ) 

RTRIMP! I )*RTRl I ) 

RTI IMP! I )*RTI ! I ) 

GO  TO  92 
END 


Figure  29  —  Flow  Chart  for  YNGFRE,  Computer  Program  for  Computing 
Natural  Frequencies  of  a  Plate  by  Young  Method 


/  READ  IN  / 


f  CALL  YINTGR  \ 
COMPUTE 

^INTEGRAL  VALUER 

REARRANGE  A- 
ARRAYS  FROM 
YINTGR  INTO 

H,  K,  E,  F  ARRAYS 

DO  400  |. 
DO  300  J 

=  1,  M 
-  1#  N 

<D 


DO  200  MX  =  1,  M 
DO  100  NY  =  1,  N 

COMPUTE  C  iyJ 

MX,  NY 


1 - 100  CONTINUE 

WRITE  AND  PUNCH 

CUJ 
MX,  NY 

- 200  CONTINUE 

- 300  CONTINUE 


400  CONTINUE 


0 


YINTGR 

INITIALIZE 

ALP  AND 

EPS  ARRAYS 

SOLVE  CLOSED 
FORM  INTEGRALS, 
STORE  ANSWERS 

IN  A-ARRAY 

Q  RETURN  ^ 


Figure  29  a  -  YOUNG 
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V. 


. _ I _ 

/  READ  LiM:  NO.  OF  TERMS 

/  in  one  direction  const: 

/FREQUENCY  CONSTANT  B-ARRAY: 
/  GASRAY  FROM  YOUNG 


CALL  VARAH1 

COMPUTE  INITIAL 
APPROXIMATE 
EIGENSYSTEM 


<CALLVARAH2  \ 
REFINE  AND  BOUND  \ 
APPROXIMATE  / 
v  EIGENSYSTEM  / 


MULTIPLY  EIGEN¬ 
VALUES  BY  CONST 
TO  GET  NATURAL 
FREQUENCIES 


WRITE  OUT: 
f SUBSCRIPT  SCHEME 
fFOR  MODE  NUMBERS;/ 
EIGENVALUES; 
FREQUENCIES 


Figure  29b  -  YEIGN 


T 


c„' m  -  c  f  B-r F*.  *  rj  *  m - c)  t  nlm ain 


format  or  n=fk 


Same  as 
(B13) 


....  I  .  a3  .  a  a  Same  as 

e0*' a--*!*  —  <i +2,  J-  «a  ♦  *(1  ■ -  c)  j  »a (B14) 


for  *  =  i  and  *  =  i 


Finally  the  main  program  compotes  the  C-matrix.  These  data  are  punched  out  on 
cards  for  use  in  a  program  for  solving  the  eigensystem. 

Only  two  cards  are  needed  for  YOUNG: 


Card 

Symbol 

Description 

Format 

1 

M 

Number  of  terms  in  z-direction,  M  -  10 

215 

N 

Number  of  terms  in  y-direction,  N  -  10; 

If  output  of  YOUNG  is  to  be  used  with 
YEIGN,  M=N 

2 

A 

Length  in  z-direction 

2F12.6 

B 

Length  in  y-direction 

The  printed  output  consists  of  the  array  of  integral  values  E  (/,  J),  five  elements  to  a 
row.  Then  comes  the  £PS- array  (values  of  €{),  again  five  elements  to  a  row.  A,  B,  M,  N  are 
printed  next.  Finally  the  array  Ny  is  both  {Hinted  and  punched  on  cards.  There  are 
N/2  elements  per  card,  (or  N/ 3  if  N  is  a  multiple  of  three)  with  the  order  cycling  first  through 
NY  =  1,/V,  then  MX  =  l,Af,  next  J  =  1,N,  and  finally  l  -  1,M. 

For  f?|'  g  ,  YOUNG  required  2  minutes  on  the  IBM  7090. 

YEIGN  Step 

YEIGN  is  a  computer  program  for  the  CDC  6600  which  uses  the  eigensystem  programs 
VARAH1  and  VARAH2.  The  latter  two  NSRDC  programs  are  FORTRAN  IV  adaptations  of 
algorithms  of  J.  M.  Varah.29 

VARAH1  computes  an  initial  approximate  eigensystem.  The  eigenvalues  are  computed 
using  the  QR  method  of  Francis30  after  the  system  is  reduced  to  Hessenberg  form.*  The 
eigenvectors  are  found  by  the  inverse  iteration  method  of  Wielandt.*  Finally  VARAH2  refines 


.  * 


i 


I 

i 


! 


♦See  Reference  33. 
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I 

I 


and  bounds  the  approximate  eigensystem  as  suggested  by  Wilkinson.31* 32  For  farther  infor¬ 
mation  about  both  the  mathematical  processes  and  the  programs,  complete  with  listings,  see 
Reference  33. 

Because  the  CDC  6600  has  a  60-bit  word,  the  high  degree  of  accuracy  needed  in  the 
inverse  iteration  might  not  be  achieved  on  smaller  word  computers.  Also,  the  largest  problem 
tested  was  a  64  x  64  matrix,  which  took  6.85  minutes. 

The  problem  to  be  solved  is  Equation  (Bll).  However,  the  double  summation  is  treated 
as  a  single  summation  for  use  in  YEIGEN.  The  problem  becomes 

|  (B  (IA,JA)  -  \I)AjA  =0,  JA  =  1, A' 

i  A  s  1 

where  £V  =  (LIM)2  (LIM  is  the  number  of  terms  p  of  Equation  (Bll);  p  must  equal  q  for 
YEIGN); 

/  is  the  identity  matrix  to  which  the  Kronecker  delta  reduces; 

A  is  the  single  dimensional  matrix  replacing  AJBn; 

B  is  the  matrix  of  two  dimensions  replacing  the  5-matrix; 

JA  is  the  subscript  replacing  at  and  n,  cycling  through  n  first,  then  m;  and 
I  A  is  the  subscript  replacing  i  and  k,  cycling  through  k  first,  then  i. 

An  example  of  the  transition  from  c£a  to  B(lAtJA)  is  shown  below,  with  LIM  =  3; 


C\\  =  B(  1,  1) 

5*1  =  6(2,  1) 

511  =  6(5,  1) 

5f!  =  6  (9,  1) 

511=6(1,2) 

511  =  6(2,2) 

# 

• 

• 

z 

C\UB(1,2) 

z 

511  =  5(5,  9) 

5||  =  6(9,  9) 

C\\  =  5(1,  4) 

511  =  6(2,  9) 

5jf  =  6(6,  9) 

5“  =  5(1,  5) 

6lf  =  6(3,  1) 

z 

C\\-B{  1,6) 

z 

5fi  =  6(7,1) 

53*1=5(1,  7) 

5l3  =  5(3,  9) 

; 

511  =  5(1,8) 

6?i  =  6(4,  1) 

C\\  =  6(8,  1) 

511=6(1,9) 

5fl  =  6(4,  9) 

511  =  5(8,  1) 

A(JA)  associates  with  m,n  in  a  similar  manner.  The  vector  A  does  have  two  sub¬ 
scripts  for  computer  storage  purposes;  however,  the  printed  output  of  the  eigenvectors  has 
two  subscripts  with  the  first  of  these  referring  to  JA.  The  eigenvector  yields  the  frequency 
modal  number  (m,n)  from  the  </.4-value  of  the  eigenvector  component  whose  amplitude  is  equal 
to  1.0.  The  subscripis  JA  are  related  to  their  respective  ( m ,  n)  values  in  the  final  section  of 
the  printout. 
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YEIGN  produces  many  pages  of  output.  The  user  should  look  first  at  the  last  few 
pages  of  the  output  for  the  eigenvalues  and  corresponding  natural  frequencies  and  for  the 
eigenvector  snbscript  scheme.  Then  the  user  shor’d  go  to  the  main  body  of  the  output  to 
locate  each  eigenvalue,  followed  immediately  by  its  eigenvector.  Now,  from  the  component 
with  the  value  of  1.0,  he  can  assign  the  frequency  a  modal  number,  as  directed  above. 

A  sample  output  for  each  eigenvalue  of  YEIGN  is  given  in  Table  13.  The  eigenvalue 
and  vector  components  are  given  with  their  error  bounds.  In  the  given  case,  the  frequency 
has  modal  number  (3,  4). 

The  data  cards  needed  for  YEIGN  are  as  follows: 


Card 

Symbol 

Description 

Format 

1 

LD< 

Limit  on  summation  of  Equation  (Bll) 

Note:  N  =  (LIM)2  is  number  of  eigenvalues 

2110 

LUP 

Number  of  iterations  for  refining  eigen* 
system.  For  engineering  purposes  LUP 
=  1  yields  adequate  frequencies 

2 

CONST 

A  I  E 

Value  of  -  |  / - 

2  a  a2  yi2y(l-er2) 

E16.8 

FREQUENCY  =  CONST  *  y/  EIGENVALUE 

3 

B(L4,JA) 

C-array  of  Equations  (Bl3)  and  (B14),  with 

JA  changing  most  rapidly;  that  is 

(JA  =  1,  N)  for  each  IA  value,  (IA  =  1,  N) 

4E16.S 

TABLE  13 


Staple  Output  Data  for  Each  Eigenvalue  of  YEIGK 


A8S 

(LAMBDA  (4)— (.704 

*BSf 

*1 

1* 

4) 

-  f  0 

*8Sf 

At 

2» 

41 

—  { 

*BSt 

At 

3* 

41 

-  f  0 

40SC 

AC 

4* 

4] 

-  c  - 

AoSt 

*f 

5» 

4) 

-  t  0 

495  ( 

Xt 

6* 

4) 

•  (  • 

AoS  ( 

xt 

7» 

4) 

-  t  0 

A9SI 

At 

8. 

41 

•  (  • 

AOS( 

Af 

9t 

4) 

-  f  0 

AoSt 

At 

10* 

4) 

-  t  9 

ABSt 

Xt 

11» 

4) 

-  t  0 

AbS( 

At 

12* 

4) 

-  t  0 

AoS  ( 

Xt 

13» 

41 

-  (  0 

A9S( 

X| 

1*» 

4) 

-  t  0 

AOS  f 

At 

159 

4) 

-  f  0 

AoSf 

xt 

16* 

4) 

-  t  0 

AOS  ( 

Xt 

17. 

4) 

-  I.  0 

AoS  I 

xt 

18. 

4* 

•  (  - 

AoS  t 

Xt 

19. 

41 

•  t  0 

AOS  ( 

*t 

20* 

4) 

«»  ( 

AUSt 

At 

21* 

4) 

•  (  0 

AOS  ( 

A* 

22. 

4) 

•  ( 

AUSt 

*1 

23* 

41 

•  C  0 

AOSt 

Xt 

2*. 

4) 

•  ( 

AOSf 

At 

25* 

41 

-  t  0 

AoS  { 

xt 

26. 

41 

-  t  0 

AOS( 

xt 

27. 

41 

-to 

AOS  ( 

*1 

28* 

*1 

-  (  • 

AbS( 

At 

29* 

41 

•  f  0 

AoS  ( 

*t 

30* 

41 

•  f  0 

AbS{ 

xt 

31. 

41 

•  f  0 

AOSl 

Xt 

32* 

41 

-  t  0 

AOS  ( 

xc 

33. 

41 

-  t  0 

ABSt 

Xt 

34*  " 

4) 

AOSf 

xt 

35. 

41 

-  t  0 

AOS  ( 

xt 

36. 

4) 

•  ' 

AMS  f 

xt 

37. 

41 

•  (  0 

AOS  C 

xt 

3U. 

41 

•  { 

AOSt 

xt 

39. 

4) 

-  1  0 

AOS  f 

Xf 

*0. 

41 

•  (  • 

AOSt 

xt 

41. 

41 

•  t  0 

AUSt 

Xt 

*2* 

41 

-  t  0 

ABSt 

Xt 

43. 

41 

-  t  0 

AOSt 

Xt 

*4* 

41 

-  (  0 

AOSf 

Xt 

45* 

41 

-  t  0 

AOSt 

Xt 

46* 

4) 

-  t  0 

AUSt 

Xt 

47. 

4) 

-  (  0 

AUSt 

-XC 

*8. 

41 

•  (  0 

AOSt 

Xt 

49. 

41 

-  (  o 

AOSf 

Xt 

60. 

41 

•  (  • 

AOSt 

Xt 

51* 

41 

-  (  0 

AoSf 

X( 

52* 

4) 

«  ( 

AUSt 

Xt 

53. 

4) 

-  (  0 

AOSt 

*( 

5*. 

41 

•  (  • 

ABSC 

X( 

55. 

4) 

-  (  0 

AOSt 

*( 

56* 

4) 

-  {  • 

».OS  ( 

X( 

57* 

4) 

-  (  o 

ABSt 

X{ 

5B» 

4) 

-  (  0 

AOSt 

Xj 

59. 

4) 

-  (  o 

AbS( 

X{ 

6  0» 

4) 

(  0 

AOSt 

X( 

61* 

4) 

-  t  0 

AbS'f 

“XT 

62* 

4) 

-  (  0 

AOSt 

X( 

63. 

41 

-  (  0 

ABSt 

X't 

64* 

4) 

-  (  0 

13340 1370142D?05)  ) 

.LE. 

.170619608-16 

1 

) 

.LE. 

.9?Q24537fc-20 

5995n7?*0300S593D-02) 

) 

.Lt. 

.778R49B0F-20 

1 

) 

.It. 

.90760Ai«bt-20 

7218*696378 J6?29D-01 1 

) 

.LE. 

,7691*3l3F-20 

1 

) 

.Lt. 

.23766565E-1V 

66939400787818310-021 

) 

.LE. 

.51119409E-20 

1 

) 

.LE. 

.6056U449E-20 

1724742043731H4O-O2) 

) 

.LE. 

. 174 9972 7E-23 

) 

) 

.LE. 

.39R509R9E-20 

) 

) 

.LE. 

.90) 19974E-21 

) 

). 

.LE. 

.422783R2E— 20 

) 

) 

.LE. 

•14423BB3E-20 

> 

) 

.LE. 

.R69A?96l£-20 

) 

) 

.LE. 

.10R142R7F.-20 

> 

) 

.le. 

.I003*l*7£-20 

) 

) 

•  LE. 

•293S3291E-21 

) 

) 

•  LE.  . 

.70895791E-20 

68745376609784760-01 ) 

) 

.LE. 

.S2992614E-20 

) 

) 

,.LE.  .. 

.1 16031 65E-19 

) 

.LE. 

•116S3756E-20 

.*29*5351707250020-01 » 
) 


.13912685146662720-01) 

_ i 

) 
) 
I 

...  , 

_  _  ) 

I 

I55is963982i92<402)' 

.S222242782532458D-C1  )' 
) 

.60344840431606750-0*) 

) 


) 

.15462A06H6B01onlO-01> 

> 

S543B6006337130-03) 
) 

■.519000**649220580-03) 

) 
) 
) 
) 

_ > 

) 

_ _ _ ) 

) 


•LE.  .  .SRR9166SE-23 
.LE.  ,2*lS*93B£-20 
.yLE»„  ..24*6*n*E“2e. 

.LE.  . 1005693*c— 20 
_.LE._  .9795531 0Er2». 
.LE.  .7735*1 37E-20 
_.LE.  . ,10R4*S66E-19 
.LE.  . 176271 56E-20 
•LE.  _ .30239793E-29 
•LE.  .67036259E-2X 

..LE. _ ,77401 143E-21 

,LE.  .29003447E-21 

..LE. _ . 3409571 8E-20 

.LE.  .97511100E-20 
„«LE.  _  ,.52*9*0*2E-20 
•L£»  .*62301 A1E-2P 

.LE,  . .30222952E-20 
•LE.  ,1717*83HE-?0 
•LE.  .11997536E-20 
.LE.  .86786923E-21 
•LE.  152900 llF-20 
•LE.  .62553825E-21 
•LE,._,177B6«90F-20 
.Lt.  ,5B73*800E-2l 
,LE._  .12221567E-20 
.LE.  .32236593E-21 
.LE.  _.50079662E-21 
.LE.  .12632940E-21 
.LE. _  .93936974E-21 
.LE.  .70007703E-21 
.LE.  .14363602E-20 
.LE.  .979500B3E-2) 
,LE.  . .869133B4E-21 
.LE.  .64457120E-21 
•LE.  .31352B43E-21 
.LE.  .24050771E-21 
.LE.  .25239636F-21 
.LE.  .21729957E-21 
.LE.  .392BC60HE-21 
.LE.  .10B45599E-21 
•LE.  .34350941F-21 
.LE.  .72577677E-22 
_.L£.  ,263*0*96E-21 

.LE.  .61120439E-22 


In  this  table,  the  eigenvalue  represents  the  frequency  with  modal  number  (3,  4). 
Notice  that  the  vector  component  AB3  (X(20,  4»  has  bounded  value  of  1.0. 
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Cloassen-Tfcoroe  Manual  Method  of  Computation 

Claassen  and  Thorne10  give  an  exact  analysis  of  the  problem  of  sinusoidal  free  vibra¬ 
tions  of  a  thin  rectangular  isotropic  plate.  For  comparison  with  the  results  of  the  present 
report,  the  frequency  parameter  K1  was  modified  Manually  to  frequency  f  using  the  formulas 

shown  below.  The  results  are  shown  in  Table  1. 
a  < 

For  —  =  i  - 1,  the  corresponding  value  Kx  is  obtained  from  a  table  in  Reference  10. 

Then:* 


nh  I  ~ 

f-*t  —  V - 

2c2  f3o  (l-o 


1  ,  &hv  [ 

For  —  >  1,  k —  <  1,  and  K'.  =  KJ&  so  that/=  K.  -  i/- 

t  h  1  0-2  y3 


2<*2  l3P»(l-'flr2) 


Sample  Problem 
Given: 


0.C313 

a  =  2  ft;  £  =  2.33  ft,  h  (half  thickness)  =  —— —  ft. 


Then: 


E  =  4175  x  106  lh/ft2,  pw=  466.56  lh/ft2,  <r  =  0.33, 
1  -  a2  =  0.8911,  g  =  32.2  fl/sec2 

it  =  4-  =  0-858 
0 


The  corresponding  value  of  K1  is  obtained  from  Table  H  of  Reference  10  by  interpolation  of 
values  of  Kx  (designated  K  in  the  reference)  corresponding  to  k  =  0.84  and  k  =  0.86  given  in 
the  table.  The  result  for  the  1, 1  mode*  is  Kx  =  3. 184789.  Then 


'u 


hn 


=  (3.184789)  (63.8047)  =  203.204 


*The  table  and  therefore  interpolation  of  tabulated  values  yield  different  values  of  for  different  modes, 
i.e.,  K  j  is  unique  for  a  particular  mode.  * 
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